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Abstract 



Lepton flavour violation and neutrino masses are a signal for new Physics beyond the 
Standard Model and are deeply related. The minimal extension of the Standard Model to 
make it include neutrino masses is not satisfactory from a conceptual point of view, since it 
requires a severe fine-tuning of Yukawa couplings. Seesaw models provide a consistent and 
natural mechanism to generate neutrino masses and require Physics beyond the Standard 
Model. In this thesis, the connections between models for neutrino masses and processes 
that violate lepton flavour are explored. The reconstruction of high energy parameters 
from neutrino data is partially possible within the framework of these seesaw models and 
it is enhanced by the knowledge of phenomena outside the neutrino sector, as it is the case 
of lepton number violation processes. Grand Uniflcation SUSY models offer a consistent 
theoretical framework for Type I, Type II and Type III seesaw models. On the other 
hand. Type II seesaw has the attractive feature of producing a lepton asymmetry through 
triplet decays which can be converted into a baryon asymmetry with leptogenesis. Thus 
being, we studied a model to explain the baryonic density of the Universe constrained by 
the phenomenological neutrino data. 



Keywords: 

Neutrino Masses; Lepton Flavour Violation: Baryogenesis: Leptogenesis; Majoron; Seesaw 
Models. 



Acknowledgments 



This thesis is a teamwork. A team made of people to which I had close contact with 
or I just know from books and scientific articles. From this second group, I select three 
particular sets of authors to whom I would like to pay tribute. The first is Steven Weinberg 
and his "Quantum Field Theory" in three volumes from where I took inspiration and 
learned much of the fundamentals of the Physics I think I know today. It was particularly 
useful his discussion on anomalies in [1] and on supersymmetry in [2]. The second set 
includes Carlo Giunti and Chung W. Kim, authors of "Fundamentals of Neutrino Physics 
and Astrophysics" [3] that was invaluable in setting the details straight and also as a 
source of experimental and bibliographical data. The third is comprises the authors of 
"Theory and Phenomenology of Sparticles" [4], which I believe found just the right balance 
between theoretical formalism and phenomenological implications of SUSY. 

From the first group of people, I begin with my Ph.D. supervisor, J. C. Romao, who, 
despite the multiple requests he has due to his academic responsibilities, was always 
available to discuss scientific matters with me and to answer all questions (even the 
unexpected) I put to him, thus making this thesis also a consequence of his perseverance. 
Next, the AHEP group in Valencia and Werner Porod in Wiirzburg also had a major 
role in the research that is presented here: the published paper [5], which is a joint work 
with J. C. Romao, A. Villanova del Moral, M. Hirsch and W. Porod, had its genesis on 
an idea by Martin and Werner, if I remember correctly, like it did also the article [6] 
published with S. Kaneko, J. C. Romao, M. Hirsch and W. Porod. In these two papers 
the publicly available Fortran code SPheno [7], which is mostly a one man's job, played 
a crucial role. I decided to reproduce here the contents of these published works because 
I find hard, if not impossible, to select the parts to which I gave a contribution even 
if a modest one and also as not to adulterate the scientific content and language style. 
The director of AHEP, Jose W. F. Valle, was also invaluable in the work presented in 
Chapter 4 as this was mainly his idea and also in an invitation to visit Valencia for 2 
weeks, where I deeply appreciated the scientific and human environment I found there. 
My gratitude also to Filipe Joaquim, now at CERN, who pointed out an inconsistency 
in the work "Leptogenesis in an A4 model" and that is described in detail in Chapter 
3. I also did benefit from discussions with Mariam Tortola, Albert Villanova del Moral, 
Avelino Vicente and Stefano Morisi. 



iv 



I also thank Gustavo Castelo Branco, director of CFTP, for accepting me at tliis re- 
search center where I liad exceptional working conditions and a very friendly atmosphere. 

Last but not least I wish to thank to my Mother, for obvious reasons and others not 
so obvious, as without her this thesis would not have been possible and also to Maria 
dos Anjos Carapau for a very professional work in reviewing the text and pointing out 
its inconsistencies in the English language. Of course, I'm the sole responsible for any 
spelling or grammar error that remains, or for some lack of proficiency. It's not an easy 
task to write in a language other than the native one. 

This thesis was done under the grant SFRH/BD/29642/2006 from "Fundagao para a 
Ciencia e Tecnologia" . 



Joao Esteves 



Contents 



Abstract i 

1 Supersymmetry and Neutrino Masses 3 

1.1 The Minimal Supersymmetric Standard Model 3 

1.1.1 Effective Potential and Higgs Bosons Masses 7 

1.1.2 Neutralino and Chargino Masses 15 

1.1.3 Lepton and Slepton Masses 17 

1.1.4 Renormalization Group Equations 18 

1.2 Neutrino Masses and Mixing 20 

1.2.1 Neutrino oscillations 20 

1.2.2 Lepton flavour violation 25 

1.2.3 Seesaw Models 27 

1.3 Solar, atmospheric, reactor and accelerator neutrinos 33 

1.3.1 Solar neutrinos 33 

1.3.2 Atmospheric neutrinos 39 

1.3.3 Reactor neutrinos 42 

1.3.4 Accelerator neutrinos 44 

1.3.5 Neutrino parameters 45 

2 Baryon Asymmetry, Leptogenesis and Dark Matter 49 

2.1 Anomalies, Baryonic and Leptonic Currents 49 

2.2 Instantons and Sphalerons 52 

2.3 Baryogenesis and Leptogenesis 56 

2.4 Cosmology and Dark Matter 58 

2.4.1 The Standard Cosmological Model 58 

2.4.2 Boltzmann equation 64 

2.4.3 Dark Matter 72 

2.4.4 The Majoron as a Dark Matter candidate 76 

3 Leptogenesis in an A4 Model 81 

3.1 The Model 81 

3.2 Triplet Decays and Leptogenesis 84 



vi 



CONTENTS 



3.3 Higgs Potential Minimization 89 

4 yl4-based neutrino masses with Majoron decaying dark matter 93 

4.1 Introduction 93 

4.2 The Model 94 

4.3 A4 Invariant Higgs Potential 95 

4.4 Neutrino parameter analysis 97 

4.5 Majoron Dark Matter 99 

4.5.1 Decay into neutrinos 100 

4.5.2 Decay into photons 100 

4.5.3 Numerical results 105 

4.6 Conclusions 106 

5 Seesaw Models and Lepton Flavour Violation 109 

5.1 Theory setup 110 

5.2 Numerical results 114 

5.2.1 LFV stau decays 115 

5.2.2 Total production cross section of X2 118 

5.2.3 Total production of X2 tiines BR to /i-r lepton pair 120 

5.3 Conclusions 122 

6 Dark Matter in Seesaw Type II Model 125 

6.1 Theory setup: mSugra and SU{5) motivated type-11 seesaw 127 

6.2 Numerical results 131 

6.3 Conclusions 139 

7 Lepton Flavour Violation and Dark Matter in Seesaw Type III Model 141 

7.1 Introduction 141 

7.2 Supersymmetric seesaw type-Ill 141 

7.3 Effects of the heavy particles on the MSSM spectrum 142 

7.4 Lepton flavour violation in the slepton sector 146 

7.5 Numerical results 147 

7.5.1 Lepton flavour violation 148 

7.5.2 Dark Matter 151 

7.6 Conclusions 155 

8 Conclusions 159 



List of Figures 



1.1 Evolution of the SM and MSSM gauge couplings with energy scale. Taken 
from [4] 19 

1.2 / — > l' ^ within a minimal extension of the Standard Model with massive 
neutrinos 26 

1.3 / — > l' ^ within a minimal extension of the Standard Model with massive 
neutrinos 26 

1.4 / — > I'-f within the MSSM 27 

1.5 Dimension 5 operator responsible for neutrino masses 29 

1.6 Type-I seesaw mechanism 30 

1.7 Type-II seesaw mechanism 31 

1.8 Neutrinoless double beta decay 32 

1.9 The pp chain of stellar thermonuclear reactions. Taken from [3] 34 

1.10 The CNO cycle of stellar thermonuclear reactions. Taken from [3] 34 

1.11 Energy spectra of neutrino fluxes from the pp and CNO chains as predicted 

by the SSM. Taken from [8] 35 

1.12 Neutrinos from cosmic rays. Taken from [3] 39 

1.13 Symmetry in neutrino flux distribution. Taken from [3] 40 

1.14 Asymmetry in atmospheric muon neutrinos as computed with Monte Carlo 
methods. Taken from [9] 41 

1.15 Ratio of measured to expected flux of different reactor experiments. The 
dotted curve corresponds to the best-fit values Am|oL = 5.5 x 10~^ eV^ 

and sin^ 26'sol = 0.83. Taken from [3] and [10] 42 

1.16 Exclusion curves for the CHOOZ experiment. Also shown is the allowed 
region for — > Ue from the Kamiokande experiment. Taken from [3] and 

[11] 43 

1.17 (a) Allowed region in the tan^ Am^ plane obtained in the KamLAND 
experiment. The lines show the regions allowed by solar neutrino data, (b) 
Combined analysis of two-neutrino oscillations with KamLAND and solar 
neutrino data. Taken from [3] and [12] 44 



viii 



LIST OF FIGURES 



1.18 Number of events as a function of the reconstructed neutrino energy. Points 
witii error bars correspond to data, tlie solid curve corresponds to best fit 
with oscillations and the dashed line refers to no oscillations. Taken from 



[3] and [13] 45 

1.19 Determination of the leading solar oscillation parameters from artificial 
and natural neutrino sources. Curves for data from solar, KamLAND and 
global analysis are shown. Taken from [14] 46 

1.20 Determination of the leading atmospheric oscillation parameters from ar- 
tificial and natural neutrino sources. Curves for data from atmospheric, 
MINOS and global analysis are shown. The dot, star and diamons indicate 
the best fit points of atmospheric, MINOS and global data, respectively. 
Taken from [14] 47 



1.21 Left: constraints on sin^ 6*13 from different parts of the global data. Right: 
allowed regions in the {612 — On) plane at 90 and 99.73% CL as well as 
the 99.73% region for the combined analysis. The dot, star and diamond 
indicate the best fit points of solar, KamLAND and combined data, respec- 



tively. Taken from [14] 48 

2.1 One loop gauge 3- vertex that contributes to the triangle anomaly 50 

2.2 Anomaly in the baryonic current 51 

2.3 Higgs Funnel 75 

3.1 One Loop diagrams contributing to asymmetry from the decay 86 

3.2 log(eA) contours as a function of the triplet mass and the triplet-lepton 
coupling A 87 



3.3 Contour plots for log^g^s as a function of Ma and A^, after the effect of 
washout has been taken in account by solving the Boltzmann equations. 

For comparison we note that (at 95% CL) we have: —9.35 < log^Q r/^ < —9.18. 89 

3.4 Contour plots for log^^rjE as a function of A^ and A^, after the effect of 
washout has been taken in account by solving the Boltzmann equations. 

For comparison we note that (at 95% CL) we have: —9.35 < log^Q r/^ < —9.18. 90 

4.1 Correlation between the neutrinoless double beta decay amplitude parame- 
ter \mee\ and the atmospheric mixing parameter. Experimental sensitivities 



are also given for comparison 98 

4.2 CP violating phase 6 and CP-invariant J in terms of the reactor mixing 
parameter. The 3 cr-excluded range for sin^ Oij is given for comparison. . . 98 

4.3 Top: One loop diagrams for the decay J — >■ 77. Bottom: Effective J77 
vertex 101 

4.4 Radiative majoron decay into photons 102 



LIST OF FIGURES 



ix 



4.5 Left paneL Tj^,y as function of tlie Majoron mass respecting Eq. (4.19) 
for tiA =1 eV (turquoise), 100 eV (daric green), lOiceV (magenta), IMeV 
(grey), lOMeV (dark blue) and 100 MeV (blaclc). Right panel: Tj^^ as 
function of the Majoron mass for the same values of the triplet vev as in 
the left panel. The upper orange shaded region is the excluded region from 
X-ray observations taken from ReL[15] 106 

5.1 Br(f2 fi + Xi) (Isft panel) and Br('r2 — >■ e + Xi) (right panel), in the 
mo, Ml/2 plane for our standard choice of parameters: fi > 0, tan /3 = 10 

and Aq = GeV, for type-I seesaw, imposing BT{fi — )> e + 7) < 1.2 • 10"^"^. . 116 

5.2 Br(f2 fi + Xi) (l^ft panel) and Br('f2 — >■ e + Xi) (right panel), in the 
mo,Mi/2 plane for our standard choice of parameters: fi > 0, tan /3 = 10 

and Aq = GeV, for type-I seesaw, imposing BT{fi — )> e + 7) < 10"^"^. . . . 116 

5.3 Br(f2 fi + Xi) (l^ft panel) and Br('f2 — >■ e + Xi) (right panel), in the 
mo, Ml f 2 plane for Ai = 0.02 and A2 = 0.5 and our standard choice of 
parameters: /i > 0, tan /3 = 10 and Aq = Q GeV, for type-II seesaw, 
imposing Br(^ e + 7) < 1.2 • 10"^^ 117 

5.4 Br(f2 fi + Xi) (Isft panel) and Br('r2 e + Xi) (right panel), in the 
mo, Ml/2 plane, for Ai = 0.02 and A2 = 0.5 and our standard choice of 
parameters: i_i > 0, tan /3 = 10 and Aq = GeV, for type-II seesaw, 
imposing Br(^ — > e -|- 7) < 10"^^ 118 

5.5 Br(f2 fi + Xi) (Isft panel) and Br('f2 e + Xi) (right panel), in the 
mo, Ml/2 plane for standard choice of parameters: /i > 0, tan/3 = 10 but 
different Aq = —300 GeV, for type-I seesaw, imposing Br(/i — )> e -|- 7) < 

1.2 • 10-^1 119 

5.6 Br(f2 fi + Xi) (Isft panel) and Br('r2 — > e -|- Xi) (right panel), in the 
mo, Ml/2 plane for Ai = 0.02 and A2 = 0.5 and standard choice of parame- 
ters: > 0, tan /3 = 10 but different Aq = —300 GeV, for type-II seesaw, 
imposing Br(^ e -h 7) < 1.2 • lO'^^ 119 

5.7 Br(f2 ^ + Xi) (Isft panel) and Br('r2 e + Xi) (right panel), in the 
mo, Ml/2 plane for standard choice of parameters: > 0, Aq = but 
different tan j3 = 30, for type-I seesaw, imposing — >■ e -|- 7) < 1.2 • 10"^"^. 120 

5.8 Br(f2 fJ'+Xi) (left panel) and Br(f2 e-hXi) (right panel), for Ai = 0.02 
and A2 = 0.5, in the mo,Mi/2 plane for standard choice of parameters: 
H > 0, Aq = 0, but different tan /3 = 30, for type-II seesaw, imposing 
Br(// ^ e + 7) < 1.2 • lO-^^ 120 



X 



LIST OF FIGURES 



5.9 Production cross section (at leading order) of X2 versus M1/2 for varying 
mo, and for our standard choice of parameters; /i > 0, tan/3 = 10 and 
Aq = GeV, in type-I seesaw (left panel) and type-II seesaw (right panel) 

for Ai = 0.02 and A2 = 0.5 121 

5.10 Production cross section (at leading order) of X2 thiies BR of X2 goii'^g to 
H-T lepton pair versus M1/2 for rriQ = 100 GeV (red), 200 GeV (green), 
300 GeV (blue) and 500 GeV (magenta), and for our standard choice of 
parameters: fi > 0, tan/3 = 10 and ^0 = GeV, for type-I (left panel) 
and for type-II seesaw (right panel) with Ai = 0.02 and A2 = 0.5, imposing 
Br(// ^ e + 7) < 1.2 • lO-^^ 121 

5.11 Production cross section (at leading order) of X2 times BR of X2 going to 
/x-T lepton pair versus A/1/2, for our standard choice of parameters: fi > 0, 
tan/3 = 10 and ^0 = GeV, for type-II seesaw, imposing Br(/i — >■ e-l-7) < 
1.2 • 10~^^, for a fixed value of rriQ = 100 GeV and different values of 

A2 = 0.1 (green), 0.5 (red), 0.9 (blue) 122 

6.1 Analytically calculated running of scalar (to the left) and gaugino mass 
parameters (to the right), leading order only. The mass parameters are 
calculated as a function of Mi^ for the mSugra parameters rriQ = 70 GeV 
and Ml/2 = 250 GeV. For M15 ~ 2 x 10^^ GeV the mSugra values are 
recovered. Smaller M15 lead to smaller soft masses in all cases. Note that 
the running is different for the different mass parameters with gaugino 
masses running faster than slepton mass parameters 130 

6.2 Top: Contours of equal dark matter density {Q^oh'^) in the (mo,Afi/2) 
plane for the "standard choice" tan ^ = 10, Aq = and > 0, for mSugra 
(left panel) and type-II seesaw with Mt = W^^ GeV (right panel). The 
lines are constant Vt^oh? with Q,^oh? = 0.1, 0.2, 0.5, 1, 2. Bottom: Range of 
parameters allowed by the DM constraint at 3 cr c.l. To the left: mSugra; 

to the right: Mt = 10^"^ GeV. For a discussion see text 132 

6.3 Allowed region for dark matter density (0.081 < Vl^oh? < 0.129) in the 
(mo,Mi/2) plane for the "standard choice" tan/3 = 10, Aq = and /i > 0, 
for five values from Mt, Mt = 10^'' GeV (red), to Mt = 10^^ GeV (cyan), 
to the left. To the right: Variation of the mass difference m^^ — m^^ (top 
lines) and of Qh^ (bottom lines), as a function of Mt for four different 
values of mo: (cyan), 50 (magenta), 100 (blue) and 150 GeV (green) for 
one fixed value of M1/2 = 800 GeV. The yellow region corresponds to the 
experimentally allowed DM region 133 



LIST OF FIGURES 



xi 



6.4 Limits for mSugra with tan j3 = 10, and > for = —300 GeV (left 
panel) and Aq = —500 GeV (right panel). The blue regions are allowed by 

the DM constraint, for the explanation of the bounds see fig. (6.2) and text. 134 

6.5 Logarithmically scaled zoom into the focus point region. In red the allowed 
region for 0.081 < flh^ < 0.129 and in cyan the allowed region due the 
variation of mtop = 171.2 ±2.1 GeV. The left panel is for mSugra case and 
the right panel for Mt = 10^^^ GeV. The other parameters are taken at our 
"standard" values 135 

6.6 Allowed region for dark matter density in the (mo, M1/2) plane for Aq = 0, 
// > and tan/3 = 45, for (from top to bottom) Mt = 5 x 10^"^ GeV (red), 

Mt = 10"(green) and Mt = lO^^ GeV (blue) 136 

6.7 Allowed region for the dark matter density in the (mo, M1/2) plane for 
Ao = 0, fi > and tan/3 = 45, for Mt = 5 x 10" GeV and (to the 
left) for three values of rritop = 169.1GeV (blue), rritop = 171.2 GeV (red) 
and rritop = 173.3 GeV (green). To the right: The same, but varying m?,. 

f^bot = 4.13 GeV (blue), rrihot = 4.2 GeV (red) and mbot = 4.37 GeV (green). 137 

6.8 Allowed region for dark matter density in the (mo, Afi/2) plane for the 
"standard choice" of mSugra parameters for Mt = 10^"* GeV. To the left: 
For one fixed value of A2 = 0.5 the allowed range for negligibly small 
neutrino Yukawa couplings (red) and Yt fitted to correctly explain solar 
and atmospheric neutrino data (blue lines). To the right: the DM allowed 
range of parameters for 3 different values of A2, A2 = 0.5 (red), A2 = 0.75 
(green) and A2 = 1 (blue). Note the logarithmic scale 137 

6.9 Allowed region for dark matter density in the (mo, M1/2) plane for our 
"standard choice" of mSugra parameters and for two values of M^: Mt = 
5 X 10" (left panel) and for Mt = 10^^ (right panel). Superimposed are 

the contour lines for the Br{fi — >■ 67) 139 

7.1 Mass parameters at Q = 1 TeV versus the seesaw scale for fixed high scale 
parameters mo = M1/2 = 1 TeV, Aq = 0, tan /3 = 10 and /i > 0. The full 
lines correspond to seesaw type-I, the dashed ones to type-II and the dash- 
dotted ones to type-Ill. In all cases a degenerate spectrum of the seesaw 
particles has been assumed 143 

7.2 Example of spectra at Q = 1 TeV versus the seesaw scale for fixed high 
scale parameters mo = M1/2 = 1 TeV, tan /3 = 10 and fi > 0. On left panel 
Mh,m^o,m^+ while on the right panel we have M^,m^o,m-+. The line 
codes are as in Fig. 7.1 144 



xii 



LIST OF FIGURES 



7.3 Four different "invariant" combinations of soft masses versus the mass of 
the 15-plet or 24-pIet, Afis = Af24- The plot assumes that the Yukawa 
couplings are negligibly small. The calculation is at 1-loop order in the 
leading-log approximation. The lines running faster up towards smaller M 

are for type-Ill seesaw, the values for type-II seesaw are shown for comparison. 145 

7.4 The limits of the invariants in seesaw type-Ill models. Left: (m|— m|.)/Aff , 
right (mg — m%)/Mf. The blue lines are for SPSS, the light blue one for 
mo = 500 GeV and M1/2 = 1 TeV, and the red one for mo = M1/2 = 
1 TeV; full (dashed) lines are 2-loop (1-loop) results. The black line is the 
analytical approximation, for comparison 146 

7.5 Values of the gauge coupling at Mqut = 2 x 10^^ GeV as a function of 
the seesaw scale, black lines seesaw type-II and green lines seesaw type-Ill 
with three 24-plets with degenerate mass spectrum; full (dashed) lines are 
2-loop (1-loop) results. For the calculation of the electroweak threshold the 
spectrum corresponds to mo = M1/2 = 1 TeV, ^0 = 0, tan /3 = 10 and /i > 0.147 

7.6 Br{^L — >■ 67) as a function of the seesaw scale for seesaw type-I (red line), 
seesaw type-II (blue line) and seesaw type-Ill (magenta line). In case of 
type-I and type-Ill a degenerate spectrum has been assumed. On the left 
panel mo = mi/2 = 300 (GeV), on the right panel mo = mi/2 = 1000 
(GeV). In both cases we take tan/3 = 10, Aq=Q and /i > 149 

7.7 5r(/i 67) versus for mo = M1/2 = 1000 GeV, tan/3 = 10, A = 
GeV and /i > 0, for seesaw type-I (solid lines) and seesaw type-Ill (dashed 
lines), for Mseesaw = 10"^ GeV. The curves shown are for 2 values of the 
Dirac phase: (5 = (red) and 5 = ti (blue), both for normal hierarchy. . . . 150 

7.8 Brij — >■ 67) versus s^g (left) and Br{T — )> /i7) versus (right) for mo = 
M112 = 1000 GeV, tan/3 = 10, v4o = GeV and > 0, for seesaw type-I 
(solid lines) and seesaw type-Ill (dashed lines), for Afseesaw = 10^^ GeV. 

The curves shown are for 5 = (red) and 5 = tt (blue) for normal hierarchy. 151 

7.9 Branching ratios for li — > Ij^ (solid lines) and li — > 3/j (dashed lines) versus 
the seesaw scale for tan/3 = 10, /i > 0, Aq = GeV, M1/2 = mo = 1000 
GeV. On the left panel we scan on Mr^ with Mr^ = Mr^ = 2 x 10^^ GeV 
while on the right panel we scan on Mr.^ with Mr^ = Mr^ = 2 x 10^^ GeV. 
The color code is red for /i — > 67 or — >■ 3e, blue for r — )> /i7 or r — > 3/i 

and green for r — > 67 or r — > 3e 152 



LIST OF FIGURES 



xiii 



7.10 Difference between tlie masses and tlie liglitest stau and tfie iightest neu- 
tralino (upper row) as well as tlie corresponding Qh^ (lower row) as a 
function of the seesaw scale. The left (right) plots are for seesaw type-II 
(III). A degenerate seesaw spectrum has been assumed in case of seesaw 
type-Ill. Ml/2 = 800 GeV, Aq = 0, tan/3 = 10 and /x > 0. The lines 
correspond to full blue line rriQ = 0, red dashed line niQ = 50 GeV, green 
dashed dotted line mo = 100 GeV, black dashed line rriQ = 150 GeV and 
orange full line mo = 200 GeV. The gray band shows the preferred range 
according to eq. (7.11) 153 

7.11 Mass of the lightest neutralino (left plot), its higgsino content (middle plot) 
and the corresponding fi/i^ (right plot) as a function of mo for a seesaw 
type-II model with Mr = 10^"^ GeV, rritop = 171.2 GeV, Aq = 0, tan /3 = 10 
and /.i > 0. The lines correspond to full blue line M1/2 = 195 GeV, red 
dashed hue M1/2 = 200 GeV, green dashed dotted hne M1/2 = 205 GeV, 
black dashed line M1/2 = 210 GeV and orange full line M1/2 = 215 GeV. 

The gray band shows the range eq. (7.11) 154 

7.12 Mass of the lightest neutralino (left plot), its higgsino content (middle plot) 
and the corresponding Qh"^ (right plot) as a function of mo for a seesaw type- 
Ill model with a degenerate seesaw scale Mw = 10^^ GeV, rritop = 171.2 
GeV, Aq = 0, tan /3 = 10 and // > 0. The lines correspond to full blue line 
Ml/2 = 400 GeV, red dashed line M1/2 = 405 GeV, green dashed dotted 
line Ml/2 = 410 GeV, black dashed line M1/2 = 415 GeV and orange full 

line Ml/2 = 420 GeV. The gray band shows the range eq. (7.11) 154 

7.13 Dark matter allowed region (in blue) for mSUGRA (left panel) and for 
type-I seesaw (right panel). The parameters are tan/? = 10, ^0 = 0, /i > 
and Mt = 10^^ GeV for rritop = 171.2 GeV. Also shown (in yellow) are the 
regions excluded by LEP (small values of M1/2), and by LSP constraint 
(small values of mo). Also shown are the Higgs boson mass curves for 



7.14 Like in Fig. 7.13 but for seesaw type-II (left panel) and type-Ill (right panel). 156 

7.15 Like in Fig. 7.14 but for Aq = —300. Seesaw type-II (left panel) and type- 



7.16 Comparison between using 1-loop (red) or 2-loop (blue) RGBs on the dark 
matter allowed region for type-II (left panel) and type-Ill (right panel). 
The parameters are: Aq = 0, ji > and Mseesaw = 10^"^ GeV, rritop = 171.2 



Mh = 110 GeV (in red) and for M^ = 114.4 GeV (in magenta) 



155 



III (right panel) 



156 



GeV and tan /? = 52 for type-II and tan /3 = 49 for type-Ill 



157 



List of Tables 



1.1 Matter superfields in the MSSM 3 

1.2 Higgs and gauge superfields in the MSSM 4 

1.3 Standard Model Predictions (BP2000): solar neutrino fluxes and neutrino 
capture rates, with la uncertainties from all sources (combined quadrati- 
cally). Taken from [16] 34 

1.4 Solar Neutrino Rates (units are in SNU for Chlorine, Gallex+GNO and 
SAGE; for *B and hep are respectively 10^cm~^s~"^ and 10'^cm~^s~"'^): The- 
ory versus Experiment. Taken from [16] with data published up to the year 
2000 - see the references therein 35 

2.1 Standard Model quantum numbers 50 

3.1 Characters of A4 82 

3.2 Lepton multiplet structure of the model 82 

4.1 Lepton muhiplet structure (Q = T3 + Y/2) 94 

7.1 Examples masses in GeV for M1/2 = 1000 GeV, tan/3 = 10, ^0 = GeV 
and fi > 0, for seesaw type-Ill for a degenerate seesaw spectrum with 
M24 = 10^^ GeV 149 



Chapter 1 



Supersymmetry and Neutrino Masses 



1.1 The Minimal Supersymmetric Standard Model 

Here we give a brief review of the Minimal Supersymmetric Standard Model (MSSM), 
following mainly [4] and [2]. The left chiral matter and Higgs superfields on one hand, 
and vector gauge superfields on the other, of the MSSM are shown on tables 1.1 and 
1.2, respectively. The Lagrangian density will be a sum of a term £susy that is fully 
supersymmetric and a term £soft that breaks supersymmetry softly: 

-^MSSM = ^SUSY + ^SOFT- (1-1) 

Moreover, Csusy is by itself a sum of a pure gauge part Cg, a matter part C]\i and a 
Higgs part Ch- 

^SVSY = + + ^H- (1-2) 

The pure gauge Lagrangian contains the spinorial gauge field strengths associated with 
SU{3), SU{2) and U{1) that define the kinetic and self interaction terms for the gauge 
bosons and gauginos: 

^9 = \j + "^w ■ Wwe. + W^Wy^ + h.c, (1.3) 

Table 1.1: Matter superfields in the MSSM 

Lepton Doublets gC/(3) ® SU{2) ® U{1) Quark doublets gj7(3) SU{2) ® U(l) 

L2=(y (1.2,-1) o^=te) (3.2,i) 

Anti-lepton singlets 5(7(3) ® SU{2) ® (7(1) Anti-quarks singlets S(7(3) ® SU{2) ® (7(1) 

(1,1,2) (3,1,-|),(3,1,|) 
Ej, (1,1,2) U_2,D2 (3,1,-|),(3,1,|) 

Er (1,1,2) (73, (3,1,-|),(3,1,|) 



4 



Supersymmetry and Neutrino Masses 



Table 1.2: Higgs and gauge superfields in the MSSM 

Vector Gauge Superfields Left Chiral Higgs Superfields 

Hypercharge 
Weak Isospin 
Vg" Color 




The gauge field strengths are obtained from tire vector superfields through the relations 

1 

4 
1 



Wa = -^-VVe-^V^e"", (1.4) 



W = --VVe^V'^e-'' , (1.5) 

with V = 2gV°'T'^ and Wa = 2gW^T°', where T"* are matrices from the adjoint representa- 
tion of the gauge group with coupling constant g. The vector superfields V are expanded 
as functions of the Grassmann variables 0, 0, in the Wess-Zumino gauge, 

v^zix, G, 0) = ea'^eA^ix) + Ma"(x) + mx'ix) + ^eeeoD'^ix), (i.e) 

which give the expansions for the field strengths W"', expressed in the chiral coordinates 

W: = + D^{y)e^ - {a'^^e)^F;M + ^eearv^\^f^{y). (1.7) 

As for the matter Lagrangian density, it is a D-term, which means that it can be expressed 
as a 4-integration in the Grassmann variables: 

_^f)tgfeV7A»+siy^y) ^. _^ Qtgfev°A"+32V"'-c?+5iy^y)g, 

with a the Pauli Matrices, A" the Gell-Mann matrices, as usual, and an implicit sum over 
i = 1, 2, 3. It gives the kinetic terms for the matter fields (SM particles and superpartners) 
and also the couplings between these and the gauge fields (gauge bosons and gauginos). 
Assuming the expansion of the left chiral superfield $ in terms of the chiral coordinate y 
and the Grassmann variable 6 

$(y, e) = <^(x) - i0a^'e^^(f){x) - ^96998'' d^(j){x) 

+ V29^{x) + ^9edf,^a^'9 + 99F{x), (1.9) 
v2 
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with (f) a scalar field, ^ a spinor field and F an auxiliary field, the component Lagrangian 
can be obtained. Finally, the Higgs Lagrangian establishes, in the D-term, the kinetic 
parts for the 2 Higgs superfields as also the couplings between them and the gauge su- 
perfields, and in the F-term, the couplings of the Higgs superfields and the matter ones, 
through the superpotential Wmssm- 

Ch = J d'e [/fte(s^^--+^i^"^)i/p + Wmssm<^'(^) + ^^IssM^^f^)] . (1-10) 
Wmssm = ea,{^iH'^Hl - ff^H^L'^E, - f^^H^Q'.D, - f^^QtH'^Uj), (LH) 

with an implicit sum over p = 1,2. The dimension four Yukawa couplings from the Stan- 
dard Model are obtained from the F-terms of the products of three left chiral superfields. 
That is why in the previous expression there are 2 delta functions in the Grassmann vari- 
ables. Also, the need for at least two Higgs superfields can be understood from the fact 
that the superpotential must have a definite chirality (holomorphy) and so it is impossible 
to define the couplings from a superfield and its complex conjugate in the same functional. 
This is also related with the cancelation of anomalies. In the previous Lagrangian Cmssm, 
there are auxiliary fields that do not have a kinematic term. In the general expansion of 
a left chiral superfield, there is an F-term that is associated with a 2-power of Grassmann 
variables, and in the general expansion of a vector superfield, there is a D-term that is 
associated with a 4-power in these variables. The functions of x that are the coefficients 
of these terms do not have a kinetic part, and so can be eliminated with the Lagrange 
equations. For a generic superpotential of the type 

W = + ^.hj^.<^, + ^^fc$,:$,<I>,., (1.12) 

where the <I>i are generic left-chiral superfields and a sum over repeated indexes is assumed, 
the equations of motion for the coefficient F- of the F-term of give 

Fi = -fi - fij<t>j - ^fijk4>j(pk, (1-13) 

with (pi the scalar components of <I>i and where it is assumed that the Lagrangian contains 
both W and W''. This is the same as taking the derivative of W with respect to and 
projecting it onto the scalar component: 



(1.14) 

6»=0 



In the same way, we have 

F=^ . (1.15) 
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Also, as can be seen by eqs. (1.3), (1.6) and (1.7), the gauge part of the Lagrangian 
contains a field, where a is a gauge index, that can be eliminated, giving the relation 



(1.16) 



assuming that the chiral superfields belong to the fundamental representation of the 
gauge group, labeled by the multi-index i. In the end, we will get a scalar potential 
y{(l)iA*j) given by 

V{^,A*i) = F}F, + ]^D'^D\ (1.17) 

The Higgs sector is specially important in what follows, so lets concentrate on the 
part of the scalar potential that is responsible for EW symmetry breaking and the Higgs 
bosons masses. Its origin is in the D-terms (1.16), associated with SU{2) and U{1): 



si 



h\hi 



h\h2 



(1.18) 



where hi are the scalar components of the superfields i/j, i = 1,2. With the relation 



(1.19) 



the previous expression for Vn can be written as 



Vn = 4 



+ 



gt + 92 



(1.20) 



Also, as we can see by expression (1.10), there is an interaction term between Hi and H2 
that originates the F-terms 



dW 



dhi 
dW 



dho 



/i/ii 



(1.21) 
(1.22) 



and gives the contribution to the scalar potential 



{h\hi + 4^2) 



(1.23) 



We must add also the part from the supersymmetry breaking Lagrangian Csoft 
^SOFT = rn\h\hi + m\ 



This Lagrangian may be obtained from first principles in a theory of supergravity, where 
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there are superfields in a hidden "world" that communicate with the MSSM superfields 
only via gravitation. We shall not go into the details here, but merely state that such a 
Lagrangian is soft, in the sense that it does not introduce quadratic divergences. Also, 
the explicit and soft supersymmetry breaking for the effective Lagrangian of the MSSM 
is the most well motivated process of breaking supersymmetry, since breaking at tree 
level with an F-term or a D-term is phenomenologically excluded, and it is known that, 
if supersymmetry is not broken at tree level, then it cannot be broken to all orders in 
perturbation theory^. 



1.1.1 Effective Potential and Higgs Bosons Masses 



The total tree level Higgs scalar potential is then 



V =Vd + V^, + Fsoft 

.2 , „2 

[h\h2 



.52 

' 2 



h\hi 



h\h2 



(1.25) 



+ {ml + 1^1^) {l-i\hi) + [ml + 1^1^) (4^2) 
+eabBnRe{hlhl). 

The condition for this potential to be bounded below can most be seen by making the 
charged components of hi and /12 to vanish and the neutral components to be equal. After 
doing this, one gets the relation 



+ m\ + m\ > Bfi, 



(1.26) 



assuming, as usual, that Bfi is positive. Also, to look for a minimum of the scalar 
potential with electromagnetic gauge invariance, we must set the charged components in 
V also equal to 0, so that V becomes 



V 



_ri+92 



+ {ml + \fif)\h'i\ 



'n| I "-2 



[1.27) 



With the vacuum expectation values of that break SU{2) (g) ?7(1), the fields can be 
expanded around it^ 

h'',=Vi + (f), (1.28) 



^Details can be found in [2]. 

^With this convention we have \/vl+v^ = 174 GeV. 
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which gives, up to second order in 0^, 

yo ^9U_9i (|^^|2 _ |^^j2^, [2Re(^,*^^ (1.29) 

+ 3l±^ [2Re(t;t</)i - v^cl^^)] + (m? + [2Re(i.t(/)i + 

+ {^1 + l/f^l') [2Re(?;2>2 + 1-^21')] (1.30) 
—B^Re{vi4>2 + W201 + 0i<?!*2) + constant. 

The minimization of this expression in order to the fields (pi forces the linear terms to 
vanish: 

{ml + \fif) V, + ^^i^ {vl - vl) V, - ^B^iV2 = 0, (1.31) 
{ml + V2 - {vf - vl) V, - \B^,v, = 0. (1.32) 

where, without loss of generahty, we have assumed that the vacuum expectation values 
are real. It is convenient to define the parameters 

tan/3=— , (1.33) 
m\=2\i_if + ml + ml, (1.34) 

and also to consider the Z boson mass 

ml = l{gl + gl){vf + vl). (1.35) 

Multiplying (1.31) and (1.32) by V2 and Vi, respectively, and summing and subtracting, 
gives 

Bfi=m\sm2^, (1.36) 
ml — ml = — {m\ + m|) cos 2(5. (1.37) 

Also, the following relations will be useful: 

ml + liJ-l"^ = ^m\-^{m\ + ml)cos2(3 (1.38) 
ml + = ^m\ + ^{m\ + ml)cos2(3 (1.39) 

From (1.34) and (1.37) we see that the square of the /i parameter can be fixed by these 
equations: 

lup = ^ (mntan^/3 - m?) ml. (1.40) 

l-tan2/3^ 2 M i; 2 ^ ^ ^ 
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There is an important feature in this equation. We see that for tan /3 > 1 ml must be 
negative. This behavior is typical from spontaneous symmetry breaking where usually it 
is necessary to put by hand a negative mass of the scalar field. However in the MSSM is 
possible to drive ml negative through the running of RGEs from some unification scale. 
This fact is known as Radiative Electroweak Symmetry Breaking. This equation will be 
important in the discussion of neutralino masses. 

We have three equations'^ for five quantities, mi, m2, B, fi and A, so only two of them 
are independent, modulus the signs. These can be chosen to be mi and m2- In the next 
section it will be clarified why, in fact, only one parameter mo is taken as arbitrary, so 
that mi = m2 = mo at some high scale. 

With the vacuum expectation values that are the solutions of (1.31) and (1.32), the 
quadratic scalar potential (1.30) is 

= lm|cos2/3 [|(/)ip - |(/)2p] +m|[Re(cos/3(^i-sin/3(/)2)]^ (1.41) 
+ + l'^2|') - \{m\ + m|) cos 2/3 [l^ip - \<j)2\^] 

— m\ sin 2/3Re {4>i(t)2) + constant. 

From this, it is possible to obtain the mass matrices of the neutral Higgs bosons. First, 
for the imaginary components of (f)i, we have 

^2^1 M(l-cos2/3) imisin2/? \ 

1^ \m\sm2(3 \m\{l + cos2P)) ' ^' ' 

This has a zero eigenvalue, which is identified with one of the Goldstone bosons from 
electroweak symmetry breaking. The other eigenvalue is equal to the trace, which is just 
m\. This will be the CP-odd neutral Higgs boson A. For the real parts of the (/)j, we get 
for the elements of the mass matrix 

(A^Re)ii =\^\{^ - COS 2/3) + ^m|(l + cos 2/3) (1.43) 
(M2Ji2 ={ML)21 = -^(mi + m|)sin2/3 (1.44) 
{ML)22 =\m\{l + cos 2/3) + ^m|(l - cos 2/3). (1.45) 



'^The Z boson mass is, of course, given by its phenomenological value and tan/3 is assumed to be an 
independent parameter of the theory, only barely constrained. 
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The eigenvalues equation gives 



m 



H 



m\ + m| + \ {m\ + ?^|)^ — 4,m\m\ cos^ 2/3 



+ m| — a/ (m^ + m\Y — 4,m\m\ cos^ 2/3 



(1.46) 
(1.47) 



With the natural assumption uia > mz,we can expand the previous expressions in powers 
of m%/m\ to conclude that 



ml < m|. 



(1.48) 
(1.49) 



The second inequality is obviously excluded by LEP data, but it should be emphasized 
that these are tree level results, and radiative corrections can change this picture in a 
significant way, as we shall see. 

Finally, the masses of the charged scalars can be obtained by taking in (1.25) the 
neutral scalars equal to their vev and the charged components to be and /i^. Thus, 
the quadratic part of the potential is 



V^ = § HKT + v,ht\' + '-^{vf - vl) {\h-\^ - \ht\') 



+ (mf + \^i\ )\h^ \ + (^2 + \l-i\ )\h2 \ + Re{Bnh^ 
With formulas (2.11), (1.38), (1.39) and the W boson mass 



(1.50) 



mw = 252(^^1+^2) 



V'^ can be written as 



which, in matrix form, is 



(1.51) 



= (ml + m|) (^sm^f3\h^\^ + cosl3^\h+\'^ + ^ sin2/3 {h^h+ + {h^)*ih+)*)^ (1.52) 



V^ = {mi + m'z) 



sin^/? i sin 2/3 
i sin 2/3 cos^/3 



{ht)\ 



(1.53) 



which gives a mass matrix with a zero eigenvalue (the other Goldstone boson) and the 
other eigenvalue equal to the trace, + m\. This is the squared mass of the charged 
scalar. 

To see how quantum corrections affect the mass of the Higgs neutral scalars, we must 
compute the effective potential [17], [1]. Lets start to consider a simple (f)'^ model with 
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Lagrangian 



1 



L = ^{d^cpf + ^m^"^ - -gcf)^ + c.t., 



(1.54) 



1 

where ^ is a real scalar field, m and g are the physical constants and c.t. denote counter- 
terms necessary to renormalize the theory. Assuming that the scalar field acquires a vev 
expanding (j) at v, with 

d) — > (j) + v. 



(1.55) 



we obtain for the physical action 

/[(/) + v] = -V4 ^-^m^w + ^gv^^ + (^i^v - j d^x(, 



2^2 



(1.56) 



Here, V4 is the four- volume that we can regularize in some way and /i^ is the field depen- 
dent mass 

- ' ' (1.57) 



2 2,^2 

H = -m + -gv . 



The effective action F and the effective potential K//, at v, are related by 

m = -v,v,ff{v). 



(1.58) 



Also, the exponential of effective action is given by the functional integral of the expo- 
nential of the action, with the integration restricted to one particle irreducible diagrams 
[1]. This allows us to conclude that the zero-loop term is just the constant in (1.56) 



iT^^){v) = -tV4-^m\' + ^v'' 



(1.59) 



and the one-loop term is 

exp ('ir(^)(t))) = J Ud(f){x)exp U j d^x {d^d{x)f - ^i^(t){:i 



(1.60) 



This can be put in a Gaussian form with an integration by parts, so that one needs to 
compute 

/■,K\ 

(1.61) 



in order to have 



InDet ( — 



,rW[^;] = -1 InDet (— 

2 V TT 



(1.62) 
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where K is the Klein-Gordon operator: 



K. 



xy 



d d 



It is useful to consider the identity 



(1.63) 



DetA = exp(lnTrA) 



(1.64) 



and to Fourier transform the operator K 



(2^)2 
.2 I ,,2 



e-^p-^—-L-e">-yK. 



(2vr)2 

(V + /^')<5'(p-g), 



a;, J/ 



(1.65) 



so that 



iV^'-)[v] = -i I d^p 



V4 



ln(^ 

TT 



d"'pln [ — (— + /i^ 



(1.66) 



where the four- volume factor arrives from the delta function in K{p,p), as shown by its 
Fourier representation. This integral is obviously ill-defined, but it can be made convergent 
with dimensional regularization [18]. Also, performing a Wick rotation, we have 



d^p 
(2^ 



In - (-/ + I? 



TV 



d'^p 



,lnU(-/ + M^ 



(27r) 
d'^PE 

d f d'^pE 



In -(p| + ^2 



-ITT 



a=0 



'da J {27iY \pI + h^ 

.d_ / (-i7r)" r(a - f ) 1 
''~d^ (^(47r)^/2 r(a) (/i2)"-'i/2 

.. r(-f) 1 

\AtiY/^ (^2)-d/2' 



a=0 



(1.67) 



noting that r'(a— d/2) is finite and T[a) — > l/a as a — > 0. Then, the effective potential 
up to one-loop is 



(1,68) 



The one-loop term is divergent in the limit e = A — d — > 0, but the counter-terms allow 
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us to eliminate the divergences [17]. Using the well known properties of the T function 

r(z + l) = 2r(^), (1.69) 
r(e) = ^-7 + 0(e), (1.70) 

where 7 is the Euler-Mascheroni constant, we see that 

^j-j) / 2W/2 ^ I - ih n^)d/2 Q 

^- - 7 + ln(47r) - ln(/.2) + (1.72) 



2(47r 



In the MS renormalization scheme, the divergent part as the ln(47r) and the 7 terms are 
subtracted using the counter-terms. In the end, we get the following expression for the 
effective potential up to one-loop: 



where an arbitrary mass scale Q was introduced for the expression to be dimensionally 
correct. This mass scale can be understood in the framework of the Renormalization 
Group equations. 

When one has, in addition, Dirac fermions that couple to the scalars through the typ- 
ical Yukawa term Xtptpcf), one has to consider the Gaussian of some operator on fermionic 
variables x- 

DetA = exp ( / dxdxxAx] ■ (1-74) 



On the other hand, it is known that the square of the Dirac operator is the Klein-Gordon 
operator, with 4 multiplicity, so we can reproduce the computation performed above, with 
some extra-factors of 2 and 1/2. Putting all together, we get for the scalar potential in a 
model with a Dirac fermion and a real scalar boson, 



T. 1 2 2 , 5 4 1 4 A /^^ 3 

Vpff = — m V H V -a in — ■ 

2 24 Mn^^ V Q 2 



16^2 ^ \ Q2 2 

with the Dirac mass rrin = Af . 



+ T7^<fln^-|L (1.75) 



In a supersymmetric theory, we have the same degrees of freedom for super-partners, 
so if there is a Dirac fermion with 4 degrees of freedom, 2 complex scalar bosons must 
correspond. This means that we must multiply the one-loop scalar part by 4, assuming a 
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full boson mass degeneracy: 

" -Tb{< f - f ) - ('■' - f ) ) ■ ™' 

where and ruf are the boson and the fermion masses. Remarkably, if they are equal, 
the two terms cancel and there is no one-loop correction. This is one more evidence for 
the need of breaking supersymmetry. We have seen that the tree level mass value for 
the lightest neutral Higgs boson is too low to be correct. This means that higher order 
corrections are essential to bring it to acceptable values, and for that there must be an 
unbalance between fermion and boson masses. 



Higher order corrections to the Higgs boson masses are dominated by loops with the 
top and the stop, so that a color factor of 3 must be included. On the other hand, the 
top mass is given by the Yukawa coupling from the last term in (1.10) and the stop mass 
comes from the SUSY soft breaking Lagrangian with a mass parameter in and also from 
the F-terms in the superpotential. Neglecting generation mixing and assuming Left-Right 
mass degeneracy, the field dependent top and stop masses are 

= uzmw (1.77) 

ml = m' + if!^,m\'. (1.78) 

The one-loop effective potential will give a contribution to the total scalar potential equal 
to 

= ^ ((^^ + u^,rm (iog^^±i|i^ - ^ 

-(/3^3)>^r(log^^^^-0). (1.79) 
This will change the minimization equation (1.32), introducing the extra term 

^^■f^'^' Ll(ln^-l]-m^,(ln%-l]], (1.80) 



87r2 I M Q2 ; ' \ Q 

which gives an extra contribution to (1.39) and, as a consequence, also to (1-41) and (1.45). 
However, this contribution partially cancels the second derivative of (1.79), leaving us with 
the positive correction 

where 
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Then the matrix element M^^)22 in (1.45) will have this additional contribution which 
changes the eigenvalues to 

ml,H =1 + m| + 

/ \ sm (3 ^ 

±i + rnlf sin^ 2/3 + (^(m| - m\) cos 2/3 + \ . (1.83) 



1.1.2 Neutralino and Chargino Masses 

The gauginos and higgsinos mass eigenstates are known as neutralinos and charginos. 
They are a superposition of the neutral gauginos and neutral higgsinos in the case of 
neutralinos and the charged ones for charginos. To obtain the relevant mass matrices, 
one has to consider the Higgs Lagrangian (1.10) and the general expansion (1.9) applied 
to the Higgs superfields Hi and H2. This generates terms that generically are of the type 

-V2g{T%,\^^,<t>*+h.c., (1.84) 

where A°, (pi and are a gaugino, a Higgs and a Higgsino fields respectively. When the 
neutral components of the Higgs fields acquire vevs, a mixed Higgsino/gaugino mass term 
appears. Thus, for the charged fermions, 

r^^gg = --^ {vi\+hl + V2\~h\ + h.c.) - (a/sA+A" + iih\h\ + h.c.) (1.85) 

where the gaugino states are A^ = l/v^(Ai =F ^'^2). Defining two component fermionic 
fields as 

^-^(^:), (1.86) 

r = f V (1.87) 



Hi 



the chargino mass Lagrangian is 



Cl,^ss = -{rY M,ij+ (1.88) 



with the chargino mass matrix 



\ v^Mm/cos/3 /X / 
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Upon diagonalization of this matrix, we obtain the mass eigenstates which are hnear 
combinations of gauginos and higgsinos. 



As for the neutrahnos, the mass Lagrangian is 



nn 

'-MASS 



--^Aa {vih\ - V2hlJ + -^Aq (vih\ - 
M2A3A3 - MiAoAo + fLhlhl + h.c. 



(1.90) 



Defining a four component fermionic field as 



/Ao\ 

A3 
hi 

\ ^2 / 



the neutrahno mass Lagrangian can be written in the compact way 

1 



nn, 
'-MASS 



with the neutrahno mass matrix 



Ml 




-Mzcpsw MzSpSw \ 

M2 MzCpCw -MzSfsCw 

MzCfsSw MzCfsCw -/x 

y MzS/^Sw -MzSpCw -fi ^ / 



(1.91) 



(1.92) 



(1.93) 



with sw = sm9w,cw = cos6'h',s^ = sin/3 and = cos/3. Upon diagonahzation, we 
obtain four mass eigenstates that are linear combinations of the neutral gauginos and 
higgsinos. 



We see that the neutrahno and chargino masses are strongly dependent on the soft 
breaking mass parameters Mi and M2 and also on the fi parameter. Since this last one 
is fixed by the electroweak breaking conditions (1.36) and (1.37) and assuming universal 
boundary conditions at the GUT scale, as it is the case in mSUGRA, in strict sense 
the neutrahnos and charginos masses depend on the scalar rriQ and fermionic M1/2 soft 
breaking mass parameters and also, of course, on tan /3. This allows us to obtain some 
limit cases for the neutrahnos. Since the /_i parameter belongs to the Higgsino part of the 
neutrahno mass matrix, we see that for fi <C Mi, M2, the lightest neutrahno state will be 
mainly a Higgsino and for the opposite condition it will be mainly a gaugino. This can 
happen in some regions of the mSUGRA parameter space, as we shall see further on. 
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1.1.3 Lepton and Slepton Masses 



To obtain the slepton masses we need to consider the relevant part of the soft SUSY 
breaking Lagrangian: 



(1.94) 



There are also F and D terms that arrive from (1.11) and (1.8) through (1.14) and (1.16), 
respectively, and that contribute to the scalar potential (1.17). The F terms are 



(1.95) 



and the D terms are 

= \g! {\hi\' - \h2\') E (l^'iP - 2|e,i?,P) + ^gl (h\ah + hlah,) V^^ak^. (1.96) 



When the neutral Higgs fields acquire vevs, the previous expressions lead to the following 
mass terms 



-C\j = V* {{Mf),, + M| cos 2(55,j) v, 

+ e*i ((Mf ),, - Ml cos 2/3(1/2 - sin^ 0^)5,, + m^J,,) e.^ 
+ ~<R [i^Dzj - ^1 cos 2/3 sin^ e^Sij + m^^Jy) e^R 
- i^ki^ij* + mei5ij/itan/3)ej7? + h. 



.c. 



Defining a six-component slepton field f 



h 
Ir 



[1.97) 



(1.98) 



where each f^ R carries a generation index i, we can write the sfermion mass Lagrangian 
in a compact way, 

-4/= Ef^Mff' (1-99) 



with Mf a 2 X 2 block matrix, each block being a 3 x 3 matrix: 



M- = ^ ^ 

\ fRL fRR 



(1.100) 
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Explicitly, for sneutrinos, 



, , / M? + l/2M|cos2/3 1 0\ , , 

M,'=i ' / ^ 1 (1.101) 



and for selectrons, smuons and staus, 



M, 



/ M?-M|(l/2-sin2 6'vF)cos2/3 1 + me^me -me(A^* + /itan/3 1) 

" \^ -(A^^ + fi* tan (3 l)mj M? - M| sin^ Bw cos 2/3 1 + me We 

(1.102) 

with (me),:j = nie^Sij. In the case of a significative mixing, either by large values of 
or of the /i parameter, the mass of the lightest slepton can be reduced considerably. 



1.1.4 Renormalization Group Equations 

It is well known that the 1-loop Renormalization Group Equation (RGE) for the gauge 
coupling in Yang-Mills theories is [17] 



dt 




C2{G) + ^ 5^ %{R) + i 5^ r„(i?) j g\ (1.103) 



where C2{G) is the quadratic Casimir factor of the gauge group G {N for SU{N)), Ti{R) 
refers to the representation constant to which the i chiral fermions belong and a refers to 
the complex scalars in the loop, with the same meaning for Ta{R). In a supersymmetric 
theory, for each chiral fermion in some representation, there is a complex scalar in the 
same representation. In this way for the MSSM we must add the second and the third 
terms in the previous expression for each i and with a = i. Also, for each gauge boson 
in the adjoint representation there will be a gaugino in the same representation and this 
means that we must add — 11/3C2(G) with 2/3C2(G). Then for SUSY Yang-Mills theories 
we have 

To get the explicit evolution of the three gauge couplings in the MSSM we must apply 
this formula to its gauge, Higgs and matter content as shown in table 1.1. For the U{1)y 
gauge group we must add for the matter and Higgs fields which gives 11 so that, 

for gi, ^ 

'\gl (1.105) 



dgl^n_ 

dt 8tt^' 



For SU{2)^C2{G) = 2 and the gauge contribution is —6. The matter and Higgs fields 
contribute with 14/2 as there are 14 SU{2) doublets in the MSSM (counting color) and 



1.1 The Minimal Supersymmetric Standard Model 



19 



T,;(G) = 1/2. So 



dgl 

dt 



792- 



(1.106) 



Finally for SU{3) the gauge part contributes with —9 and there are 12 color triplets and 
anti-triplets that contribute to '^^T.i{R) with 6. Then 



dt 



87r2 



9t- 



Writing generically 



dt 8tt' 

we can integrate this equations to give at 1-loop 



9^ 



9iM 



(l-^Ai?K^)ln(^/Mo))' 



(1.107) 



(1.108) 



(1.109) 



with t = ln(/i/yUo). The evolution is shown on figure 1.1 with the usual SU{5) normal- 
ization for the Ui gauge coupling, gf — )> 3/5gl. Note the unification of the three gauge 
couplings on an energy scale around 10^^ GeV. For the sake of completeness we present 



1.5 - 



1.0 



-to 



0.5 



0.0 




5/3 g^=gf 



10+ 



10 



12 



10 



16 
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Figure 1.1: Evolution of the SM and MSSM gauge couplings with energy scale. Taken from 
[4]. 



here other RGBs in the MSSM. We list only the ones relative to the third generation. 
Details can be found on [4]. 



Yukawa Couplings 



dft _ 


ft 


dt 


167r2 


dh_ 


fb 


dt 


167r2 


dlr_ 


fr 


dt 


167r2 



r -[2 , r2 ^^2 

Gft + fb - y ^3 



352 



13 



Qfb + ft + fr- y 53 - 3^2 - 



^fl+^fl 



^9l 



7 2 

9^1 



(1.110) 
(1.111) 
(1.112) 
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/i parameter 

^ i^f^ + sfS + fr-^dl-gl) (1.113) 



dt 167r2 

Trilinear soft SUSY breaking parameters 



dA' 


1 


~dt 


~ 8^ 


dA^ 


1 


dt 


~ 8^ 


dA"- 


1 


dt 


~ 8^ 



+ flA' - ^{glM, - 3glM, - fgjM, ) (1.114) 



1 fi 7 

r2 Ab . f2 At , f2 AT ^'^ „2 n ,r Q^2»^ '2, 



Bilinear soft SUSY breaking parameter 

dB 1 



f^Sr - 3gi\Mi\ - gt\M,\' - -gfSt (1.120) 



(1.115) 
(1.116) 



8^2 (-3/'^* - ^fSA' - f^A^ + 3glM2 + glM^) (1.117) 



Soft SUSY breaking Higgs boson masses 

'^'^1 1 ( o f2o I -p2 o o „2 1 11 .r2 I ,.2 1 , r 1 2 1 



¥ = 8^ (3/f 5* - 35IIMII - ^,?|Mip + \glS'^ (1.119) 



Soft SUSY breaking third generation slepton masses 

dm? 1 / 1 

'3 _ ^ f2Q Q„2|,^2| ,.2|i,r|2 -^^2c2 



dt 87r2 V ' 2 

diTi— 1 

^ = ^ (2/,25, - 4gl\M,\' + glSl) (1.121) 
Soft SUSY breaking gaugino masses 

^-^--•^•^ 

1.2 Neutrino Masses and Mixing 
1.2.1 Neutrino oscillations 

Lepton flavour changing currents are absent in the Standard Model. This happens because 
neutrinos are massless and it is always possible to choose a physical basis where the 
leptonic Yukawa couplings are diagonal. However, the evidence of neutrino oscillations is 
an indirect proof for neutrino masses. In fact, lets assume that a neutrino of a specific 
flavour Ui is created at t = and propagates as a superposition of mass eigenstates Ua- 
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We would like to know the probability of finding the flavour j at a time t. The flavour 
and the mass states are related through a change of basis performed by a unitary matrix 
U: 



> = ^ f^m W) , i = e, i^T. (1.123) 



As a simplifying assumption, we take a beam of definite momentum p, related to the 
energy of each neutrino with a definite mass through the familiar relation 

= VW+^a- (1-124) 
For quasi relativistic neutrinos, we can expand this in powers of m^/[p|^, obtaining 

^»-b1(l + ^)- (1-125) 

The time evolution is given by the usual exponential factor and the fiavour state at a time 
t is 

Kt)) = ^e-'^"'U,aWa), (1.126) 

a 

SO that the ampHtude for finding a flavour j at time t is given by 

(t) = 5^ e-^»*/7;,C/„ (1.127) 

a 

and the conversion probabifity is 

= Yl \U]aU^aU*,U,h\ cos{{E, - E,)t + (|),h^J), (1-128) 

ab 

with 

<f>abij = Aig{U;,U,,U*,U,b). (1.129) 
Using now equation (1.125), we obtain the result 

P.,^.^{X) = \U*,U,aU*,U,b\ C0S(^ + (l^aUj), (1-130) 

ab 

where t has been replaced by x and the oscillation length Lab is defined by the formula 

Lab = ^ (1.131) 

'<b 

4irE 
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with E some common energy scale and m^^ the square masses difference 

rnlb = ml-ml (1.133) 

It is now obvious that, in the limit of zero masses, the conversion probability is zero. 
Also, neutrino masses cannot be fully degenerate, because in that case we also get a null 
oscillation probability. 

It is possible to give some variations of formula (1.130) that are useful in some appli- 
cations. From (1.128) we can write 

P,^^,^,(t) = U*JJ,aUlU,,eMKE, - E,)t), (1.134) 

ab 

From the unitarity relation 

Y.U*JJ.,aU*,U,, = 5,, (1.135) 

ab 

and the following fact 

2 UlU,aU*,U,, = U*JJ,am,, + Y U*JJ^aU*,U,, + Y U*,U,aU*,U,, (1.136) 

ab a^h b>a b<a 

we know that^ 

Y \U^a?\U,a? = " 2Re ( f/Jat/^f/^f/, J (1.137) 

a \b>a J 



and then 



P,,^,^,(t) = Re [y U*JJ^aU*,U,')j COS - Im U*aU,aU*,U,^ sin 



2E 



^ Y P'M^ia? + 2Re \ YUlcP^-^*bUob\ COS ( 
a \b>a / ^ 



2E 



~. b>a 



2E 



^ - 21^^ U*JJ^aU*,U,)j [l - cos ) 
2Im [Y U*JJ,,UIUA sin (1.138) 



*In the last term just change dummy indices (a, 6) — > (&, a) 
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or else 



P^^^^^it) = 5,, - 4Re 5^ U*aU;aUlU,, sin' 

\b>a / ^ 

- 2Im U*^U,M:,uJ sin 



,b>a 



(1.139) 



If we want to compare this result with the one from oscillations for anti-neutrinos, we 
must start by taking the complex conjugate of formula (1.123): 

a 

Then we can reproduce the same computations but with Uia replaced by [7*^ and vice- 
versa. In the end we get for anti-neutrinos 



V b>a 

-2 



iE 



+ [Y. U*aU.aUUJ,ij Sin . (1.141) 

This means that a way to test CP violation is to measure the asymmetry 

= 4Im (^g UlU,M:,U,}j sin . (1.142) 

We see that for CP violation to occur the imaginary part must be non-zero. It can be 
shown [3] that the quantities 

s,fc.,,J = Im {U*JJ,aUlU,b) (1.143) 
are all equal up to a sign Sab-.ji and J, given by 

J = lm{U,,U,2U;^U:,), (1.144) 



is known as the Jarlskog invariant, because it is invariant by a re-phasing of the neu- 
trino fields or equivalently by a re-parametrization of the mixing matrix. A convenient 
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parametrization of this matrix is given by 

(1.145) 



U 



( C12C13 S12C13 Sige"'"^" \ 

-S12C23 - Ci2S23Sl3e**'' C12C23 - Si2S23.Sl3e**'' S23C13 
\ S12S23 - Ci2C23Si3e'*i3 -C12S23 - Sl2C23Sl3e'*'^ C23C13 / 



where Cab = cos 6ab, Sab = sin 6ab, with the mixing angles varying in the range < < 7r/2 
and the Dirac CP-phase varying in < 5r3 < 27r. In this case the Jarlskog invariant has 
the expression 

J = - sin 2^12 sin 2^^23 cos ^13 sin 2^^l3 sin 5i3. (1.146) 



In the approximation of only two neutrino flavours the mixing matrix takes the simpler 
form 

(cos 61 sin 6 \ 
(1.147) 
— sin 9 cos 6 I 

and we can give also a simpler formula for the transition probability from (1.134): 

P..^., = \ sin' 2^ (l - cos ^) , ^ ^ J (1.148) 

or 

P,^.^, = \ sin^ 2esm^ , z ^ j. (1.149) 

The survival probability is obtained in the obvious way: 

P,.^,. = 1 - P,,^,^ (1.150) 
= l-isin^2esin^(^), z/j. (1.151) 

Noting that the average in x of the cos function is zero, we get 

(P.,^.,) = ^sin2 20 (1.152) 

These results can be useful if some experiment is only sensitive to the mass differences 
between two states, as it is already known that Aml^j <C l\m\^^. 

In reactor oscillation experiments neutrino energies are of the order of 1 MeV and it 
is convenient to write the transition probability in these specific units: 

. , /1.27m? (eV^)a:-(m)\ 
P,,^., = sin^ 9 sin^ ' (1.153) 



^(MeV) 

In real experiments it is not possible to determine E and x with 100% accuracy, so that 
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it is necessary to average the probability formula as 



with 



{P..^..) = \ sin^ (l - (cos )) , ^ ^ J (1.154) 



This average can be computed analyticahy if the density (^) is gaussian with an average 
value (-1) and a standard deviation cJx/e- In other cases it must be computed numerically. 

From formula (1.154) we get an upper bound for sm^ 6 in terms of the maximum 
oscillation probability and the cos average: 

2^111 ax 

sin' 9 < "'^"'^ , . (1.156) 



1- (cos(^)) 



This allows to make exclusion plots for Am' = m'j, vs sin' 2^, like the one in figure 1.16 
if the maximum probability is known. 

Two limiting cases are worth mentioning. The first one is 

1^i^^)^^^^l^^2e<P--^,^ (1.157) 
which gives the most stringent bound on sin' 29. The other is 

^ (|) » vr ^ sin' 29 < 2P-- ^, (1.158) 
In this case the argument of the cos function oscillates rapidly which gives a zero average. 



1.2.2 Lepton flavour violation 

In principle we could extend minimally the Standard Model, including an SU{2) singlet 
for the right handed neutrino and generate neutrino masses through the familiar Higgs 
mechanism. Even though this process wouldn't be completely satisfactory because it 
would require a fine tuning of Yukawa couplings, we would have lepton flavour violating 
(LFV) processes like the one in fig. 1.2. However, these processes would be highly 
suppressed due to the GIM mechanism, in virtue of the smallness of neutrino masses. To 
see this, note that the amplitude for the process in fig 1.2 is proportional to 



(1.159) 



26 



Supersymmetry and Neutrino Masses 




Figure 1.2: I — > I 7 within a minimal extension of tlie Standard Model with massive neutrinos. 




Figure 1.3: / — > I 7 within a minimal extension of the Standard Model with massive neutrinos. 



with 1^ and rrib the momentum and the mass of the loop neutrino and U the matrix that 
diagonalizes neutrino masses, considering a leptonic weak basis where the charged leptons 
mass matrix is already diagonal. Because U is unitary, we clearly see that, in the limit 
rrifc = or fully degenerate masses, we get no lepton flavour violation. Now, we focus on 
the first factor of the previous formula and, assuming that neutrino masses deviate little 
from a common value m, expand it in powers of Am^/(/c^ — m^), with Ami — "^6 ~ "^^5 
to get an amplitude proportional to 



b 



This illustrates the fact that LFV processes are encoded in the off-diagonal elements of 
the mass matrices. Usually, one adopts the formalism of mass insertions, where these mass 
matrix elements are considered as interaction vertices that flip the flavour and propagators 
are functions of some common mass scale. In this way, we could describe the process with 
the diagram of fig. 1.3. 

In the MSSM the situation is rather different. Fig. 1.4 shows two typical contributions 
for the LFV process / — > I 7, mediated by neutralinos and charged sleptons and charginos 
and sneutrinos, respectively. 




Figure 1.4: I — > I'-f within the MSSM. 



1.2.3 Seesaw Models 

The Dirac Lagrangian 

Cd = ii'^^d^il) — rwipxlj (1.161) 

with the decomposition of the Dirac spinor ijj in left and right chiral states '0L and iIjr 
can be written as 

= i^LYd^%l)L + i^RYd^%l)R - mipLi^R - rmljRil}^. (1.162) 

The Lagrange equations for -0^ and -0/? give a set of two coupled equations 

il^'d^^L = mi)R (1.163) 
t^d^ijR = m^jL (1.164) 

which decouple in the limit m = 0, giving origin to what are known as Weyl spinors. 
However, even in the massive case it is not mandatory to have a Dirac spinor with four 
degrees of freedom. With the Majorana condition 

^ljR = Ci^^, (1.165) 

where C is the charge conjugation operator, it is possible to show that the equations 
(1.163) and (1.164) are equivalent. The decomposition of the Dirac spinor 

i^ = i^L + i^R (1.166) 
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translates into 

which means that a Majorana fermion is self conjugate; 



;i.i67) 



(1.168) 



Thus only neutral particles can be of Majorana type, as it could be the case of neutrinos. 
The Dirac mass term 

/:™ = -m^Vfl + H.c. (1.169) 



in the case of a Majorana particle gets transformed into the more economical way 

= (Ml + . (1.170) 



The 1/2 factor accounts for the fact that xpl and xpL are not independent. Noting the 
relation 

^=-^IC-\ (1.171) 

and the properties of the charge conjugation operator, the Majorana mass term is usually 
written as 



i^lCijL + H.c. 



The full Lagrangian for a Majorana particle is then 
1 



M 



(1.172) 



;i.l73) 



which can be written in a more compact way: 



(1.174) 



where -0 is the Majorana fermion (1.167): 



;i.i75) 



We have chosen to express -0^ in terms of -0^ and obtain the Lagrangian as a funtion of 
but we could equally well have chosen the opposite procedure and obtain a Lagrangian 
with respect to This is more suited for models that include Majorana right handed 
neutrinos that are singlets under the gauge group of the Standard Model. This is one 
particular possibility for the implementation of the effective dimension five operator that 
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Figure 1.5: Dimension 5 operator responsible for neutrino masses. 



violates lepton number and generates neutrino masses (fig. 1.5), 

O = jLHLH, (1.176) 

long ago proposed by Weinberg [19]. For these theories, after electroweak symmetry 
breaking, a Dirac mass matrix rriD that couples uj^ and Uji is present with generation 
indexes implicit, as well as a Majorana mass matrix M for Ufi. In the static limit for Uji 
the Lagrange equation gives 

UR = -M-^mlCW^, (1.177) 
which replaced in the Dirac equation for 

i^d^T^L - rriDUR = (1.178) 

gives 

i^df^i^L - M^jCuf = 0, (1.179) 
that it is the equation for Majorana neutrinos with a mass matrix Mj given by 

= -niDM-^m^^. (1.180) 

This is the celebrated seesaw mechanism, in this case type I seesaw (fig. 1.6). There are 
other possibilities for generating the operator (1.176). In Type-II seesaw (fig. 1.7) an 
SU{2) triplet T mediates through a t-channel the LH scattering which, after electroweak 
symmetry breaking, generates neutrino masses. The relevant couplings are 

-Cii = XiL^C (irsT) L + Aji/^ (iraT^) H + H.c. (1.181) 



For light neutrino masses of the order of 0.05 eV, the scale of the new physics A in (1.176) 
can easily be estimated. In Type-I, mjj is of the order of electroweak symmetry breaking, 
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so 

A ^ ^ ilL X 10^ ~ 10^^ GeV. (1.182) 
0.05 ^ ^ 

Of course, this depends on the magnitude of the Yukawa couplings in the interaction 

terms between the Higgs doublet and the left leptons, so there is some freedom in the A 

scale. On the other hand the right neutrinos are singlets under the Standard Model gauge 

group, which means that there is no symmetry that prevents their mass from being high. 

So this upper scale provides a nice justification for the smallness of light neutrino masses, 

which otherwise would require a severe fine tuning of Yukawa couplings if it were to be 

generated by the usual Higgs mechanism. 

The estimate for Type-II seesaw models shows that 

M| ~ X 10^ GeV. (1.183) 

Supersymmetric models require at least 2 triplets because of the holomorphy of the su- 
perpotential and the cancelation of anomalies. In that case A2 is proportional to one of 
the triplets mass My, leaving us with roughly the same estimate as in Type-I. But in 
non supersymmetric models Mt can be as low as 10* GeV (note however that A2 has the 
dimension of mass). 

The seesaw models predict that light neutrinos are of Majorana type. The neutrinoless 
double beta decay 2/3oi> (fig. 1.8), if observed, would be the final answer to the Dirac 
or Majorana nature of neutrinos. In fact, this process is equivalent to neutrinos being 
Majorana particles, as it was first proved in [20]. Even if 2/3q^ is an effect of some new 
physics beyond the Standard Model, like supersymmetry, Left-Right symmetric models, 
or other, the same physics would imply a Majorana mass term for neutrinos. 

The counting of independent parameters depends on the seesaw type one is considering. 
For Type-I, from equation (1.180), we know that Mf is the mass matrix of Majorana 
neutrinos so that it can be diagonalized via an orthogonal transformation with a unitary 
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s ^ 



Figure 1.7: Type-II seesaw mechanism. 



matrix U 



= diag(mi, ms, mg) = U^M'^U 



(1.184) 



with m.;, i 



1,2,3 in general complex. An miitary nxn matrix has "^"""^^ independent 
modulus and "^"^"^^^ independent phases. This means that U will have 3 mixing angles 
and 6 phases. Five of these phases can be absorbed by the charged lepton and neutrino 
fields redefinitions. Not all phases can be absorbed because there is an overall phase that 
is related with the lepton number and that leaves the leptonic charged current invariant. 
On the other hand to change to the physical masses we must multiply Dj by a diagonal 
matrix of phases. This means that 



y = t/ • diag(e-**i/2^ e-**^/2^ 1) 



(1.185) 



where we have factored out one phase and U is of CKM form 



U 



( C12C13 

-S12C23 - Ci2S23Sl3e 



\ S12S23 - Ci2C23Sl3e"^" 

Then we have 



S12C13 Si3e 



-id 



C12C23 - Si2S23Sise S23C13 
-C12S23 - Sl2C23Sl3e'*" C23C13 J 



diag(mi, m2, m3) = V^MjV, nii > 0. 



;i.l86) 



(1.187) 



Now we must repeat the procedure for the right-handed neutrinos. Then the overall 
counting is, for the light neutrino sector, 3 positive eigenvalues for the mass matrix, plus 
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d 



Figure 1.8: Neutrinoless double beta decay 

three mixing angles, plus 3 CP phases (1 Dirac and 2 Majorana), which gives a total of 
9 independent parameters. In the end we have 18 independent parameters for Type-I 
seesaw. This can be put in a more formal way. In [21] it is shown that the leptonic 
Yukawa couplings Y^, responsible for the Dirac mass term m.o are given by 

n = D^,RD^V^ (1.188) 

where D ^ is the diagonal mass matrix with positive eigenvalues for the right-handed 
neutrinos, is the same but for the light neutrinos, i? is a generic 3x3 complex or- 

thogonal matrix and V is given by (1.185). Of course, this is the most general situation; in 
many applications it is possible to reduce drastically the number of parameters, imposing 
certain assumptions. For instance we can assume that the R matrix is the identity, which 
amounts to say that the right-handed neutrinos are already in the mass base. This was 
done in [5, 6], as it is mentioned in Chapters 5 and 6. In [22] we made an even more 
drastic assumption, that the right-handed neutrinos are fully degenerated because of the 
symmetry group, which reduces the parameters in the right-handed sector from 9 to just 
1. 

With respect to Type-II the number of independent parameters is smaller. We have 
also the 3 positive eigenvalues from the light neutrinos mass matrix, plus the 3 mixing 
angles and the 3 phases. Then we have the complex coupling A2 and the triplet mass 
Mt- But two phases can still be removed by field redefinitions. This means that the 
Yukawa couplings Ai in (1.181) will be a function of 9 parameters from the light neutrinos 
as in Type-I plus two real parameters, which gives a total of 11 real parameters. Both 
Type I and Type II seesaw mechanisms have too much degrees of freedom to be fully 
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reconstructed purely from neutrino data. In fact at the moment we only have rigorous 
data from the neutrino square mass differences (3 parameters), 2 mixing angles and an 
upper bound on the third mixing angle. For now, no information about the CP phases 
is available. This means that some data outside the neutrino sector (like lepton flavour 
violation constrains) will always be necessary to make the models more predictive. 

1.3 Solar, atmospheric, reactor and accelerator neutrinos. 

In this section we present a brief survey on neutrino properties including the most relevant 
experiments that lead to the discovery of neutrino oscillations. This follows closely [3], 
where full references can be found for the experimental results presented here. We end 
with a review of neutrino parameters obtained from a global analysis of the data from the 
major neutrino oscillation experiments. 

1.3.1 Solar neutrinos 

The sun is a powerful source of neutrinos through the various thermonuclear reactions 
that occur in its core. It is estimated that the overall neutrino flux from the sun is 

$® = 6.54 X 10^° cm-2 s~\ (1.189) 

There are two major groups of thermonuclear reactions responsible for the emission of 
neutrinos and radiation: the pp-chain (fig. 1.9) and the CNO cycle (fig. 1.10). The net 
result of these reactions is the conversion of four protons and two electrons into a '^He 
nucleus plus two electron neutrinos 

4p + 2e- — >^He + 2zye + Q, (1.190) 

where Q is equal to 26.731 MeV and is released as radiation or neutrinos kinetic energy. 

Solar neutrinos were first observed during the 60's in the Homestake experiment, a 
radiochemical apparatus based on inverse /3 decay: 

z/e + ^^Cl — >^^Ar + e- (1.191) 

with a neutrino threshold energy = 0.814 MeV. Therefore this experiment was only 
able to detect intermediate and high-energy neutrinos (see fig. 1.11). 

The Homestake experiment, which was proposed in 1964 and built in the period 1965- 
1967, released its first data in 1968, declaring a solar neutrino flux less then 3 SNU^, well 
below the rate predicted by the Standard Solar Models (SSM) (see tables 1.3 and 1.4). 
Since the experiment was based on the chemical extraction of the radioactive isotope "^^Ar 

^1 SNU = 10-36 events atom-is-^ 
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Figure 1.9: The pp chain of stellar thermonuclear reactions. Taken from [3] 



(«0) 



"^O + p^ "F + 7 



l5N + p^"C + 'He (cn) i=C + p^"N + 7 



"N + p ^ '»0 + 7 



Figure 1.10: The CNO cycle of stellar thermonuclear reactions. Taken from [3] 

Table 1.3: Standard Model Predictions (BP2000): solar neutrino fluxes and neutrino capture 
rates, with la uncertainties from all sources (combined quadratically) . Taken from [16]. 



Source 



Flux 



5.95 (l.OO^-i) 



CI Ga Li 
(SNU) (SNU) (SNU) 



PP 

pep 

hep 

^Be 

8B 

13N 
150 



X 10-^ (i.oo;h1^) 

10-7 
10-1 (1.001°}° 



1.40 
9.3 > 
4.77 

5.05 X 10-4 (l.OOlEj ?^ 
5.48 X 10-2 [imtri]) 

4.80 X 10-2 {imtoil) 

5.63 X 10-4 (l.00l°-|) 



0.0 


69.7 


0.0 


0.22 


2.8 


9.2 


0.04 


0.1 


0.1 


1.15 


34.2 


9.1 


5.76 


12.1 


19.7 


0.09 


3.4 


2.3 


0.33 


5.5 


11.8 


0.0 


0.1 


0.1 



7.6l[ii 128f? 52:3^ 



Total 



which is produced in very small quantities, the uncertainties with the data were large and 
dominated by statistical fluctuations. However the accumulation of data over 23 years of 
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Table 1.4: Solar Neutrino Rates (units are in SNU for Chlorine, Gallex+GNO and SAGE; for 
and hep are respectively 10^cm~^s~^ and 10"^cm~'^s~^): Theory versus Experiment. Taken 
from [16] with data published up to the year 2000 - see the references therein. 



Experiment 



BP2000 



nrs 
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hep-Super-Kamiokande 
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5.05 
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2.80 [1.00 ±0.14] 
2.40 [1.00±1H^] 
11.3(1 ±0.8) 
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Figure 1.11: Energy spectra of neutrino fluxes from the pp and CNO chains as predicted by 
the SSM. Taken from [8] 



operation allowed to reduce the statistical uncertainty to an average rate of 6%. The solar 
neutrino flux measured at Homestake is about 1/3 of the predicted by the SSM, with a 
discrepancy of more than 3a. 

In the 90's other experiments (Gallex/GNO and Sage) appeared based on the process 

Ue + ^^Ga — >^^Ge + e- (1.192) 

that has an energy threshold E'^ = 0.233 MeV. This allowed to detect solar neutrinos 
from all sources (see fig. 1.11). The GALLium Experiment (GALLEX) started operating 
in May 1991, in Gran Sasso laboratory, Italy and it used a detector with 101 tons of liquid 
gallium chloride, GaCla-HCl solution. This experiment was followed by Gallium Neutrino 
Observatory (GNO) which operated from May 1998 until April 2003. The total collected 
data of the two experiments from May 1991 until April 2003 resulted in an average solar 
neutrino capture rate [23]: 

^Ganex/GNO ^ gg 3 ^ ^ ^ ^^_^g3^ 

which is about one half of the predicted by the SSM (see table 1.3). 
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Another Gallium experiment is SAGE (Soviet American Gallium Experiment) that 
started collecting data in 1990. It is located in the Baksan Neutrino Observatory of the 
Russian Academy of Sciences, in the northern Caucasus mountains. The average neutrino 
capture rate obtained from the data collected in the period of January 1990 to December 
2001 is [24] 

i?^age^ = 70.8tli SNU (1.194) 

and it is in agreement with GALLEX/GNO experiments. It is about one half of the 
predicted by the SSM. 

Water Cerenkov detectors provide another way of detecting indirectly solar neutrinos. 
If a charged particle travels with velocity d > 1/n in a medium with a refraction index n 
then it emits a cone of light that can be detected by arrays of photomultipliers, allowing a 
precise determination of the arrival time which is to say of the interaction point. Examples 
of water Cerenkov detectors are Kamiokande, Super-Kamiokande and SNO. 

Kamiokande (Kamioka nucleon detector) was built initially to search for proton decay. 
It started operating in 1983 and was upgraded in 1986 to Kamiokande-II in order to be 
able to observe *B solar neutrinos, which generate events in the detector with an energy 
around 10 MeV. In 1990 a new upgrade allowed to reduce the electrons threshold energy 
down to 7.0 MeV by the replacement of 100 photomultipliers, which was the start of 
Kamiokande-III. 

The solar neutrino flux is measured by the elastic reaction 

Ua + e~ — > Va + e~ (1.195) 

that is sensitive mostly to electron neutrinos because the cross section is about six times 
larger then for muon and tau neutrinos. The average neutrino flux measured in 
Kamiokande from January 1987 to February 1995 is [25] 

$Kam ^ 2.80 ± 0.38 X 10^ cm-^s-^ (1.196) 

It shows a discrepancy with the SSM of about one-half (see table 1.4). 

In April 1996 a new experiment in the Kamioka mine started operating, about 500 m 
from the cavity of the Kamiokande detector that contains now the Kamland experiment. 
This was Super-Kamiokande, still a water Cerenkov experiment, that went through its first 
phase from April 1996 until July 2001. It collected solar neutrino data with an electron 
recoil threshold energy of Ef^ = 6.5 MeV for the first 280 days and E^^ = 5.0 MeV in the 
remaining 1216 days. The corresponding neutrino threshold energies are E^^ = 6.2 MeV 
and Ef = 4^.1 MeV and this means that Super-Kamiokande was also only sensitive to ^B 
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solar neutrinos. The measured flux is [26] 

$||^ = 2.35 ± 0.08 X 10^ cm-2s-\ (1.197) 

which once again correspond to one-half of the SSM flux. It was possible to compare day 
and night fluxes which could be a sign for matter effects, since night neutrinos must cross 
the Earth before entering upwards in the detector. The results are [26] 

$SK,day ^ ^2.32 ± 0.03);°:°? X 10' cm-h-' (1.198) 
^SK.night ^ ^2.37 ± 0.03)1^8 X 10' cm-^s-i (1.199) 



with an asymmetry 



^SK.day _ ^SK,night 
8b 8b 



Aay.night " ^ / sK,day , ^SK,night\ ^^'^^^^ 
2 y^SB + '^'Sb ) 

= - 0.021 ± o.o2o;EJ:E^}^ 

that is consistent with a null result. 

A variation of the water Cerenkov technique was pursued by the Sudbury Neutrino 
Observatory (SNO) in Ontario, Canada. In this apparatus neutrino interact with heavy 
water (D2O), which allows to detect neutrinos with three different processes: 

CC: Ue + d — >p + p + e~ 
NC: Ua + d — > p + n + Ua 
ES: Ua + e~ — > + e~ 

The charged current reaction has a neutrino threshold energy E't'^'^c _ ]^ 442 MeV but 
because of the high background at such low energies the effective neutrino threshold energy 
is about 7 MeV, which means that the charged current is only sensitive to neutrinos. 
The neutral current reaction allows to detect the three neutrino flavours and it revealed 
to be extremely useful in solving the observed deficit of solar neutrino fiuxes, known as 
the Solar Neutrino Problem. It has a neutrino threshold energy _£;*h,NC _ 2.224 MeV and 
so it is also only sensitive to flux. The elastic scattering is the same process that is 
used in H2O Cerenkov detectors and it has a neutrino threshold energy E'^'^'NC = 5.7 MeV 
and once more it is only sensitive to the channel. 
The SNO experiment is divided in three phases: 

D2O Phase: the flnal neutron in NC is detected through the reaction 

n + d — ^ + 7 (6.25 MeV) (1.201) 
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This phase operated during 306.4 hve days, from 2 November 1999 until 28 May 
2001. 

NaCl Phase: in the so-called salt phase 2 tons of NaCl were added to the heavy 
water in order to detect the neutron from NC, through the reaction 

n + 3^01 — > ^^Cl + several 7 (8.57 MeV) (1.202) 

This improved significantly the efficiency since the photon distribution from this 
process is very different from the one in Cerenkov radiation, which allows to distin- 
guish between NC and CC reactions with good accuracy. 

Third Phase: in this phase, a grid with three hundred "^He proportional counter 
tubes was introduced in the heavy water tank. This isotope has a large cross-section 
for the capture of thermal neutrons, producing an energetic pair proton-triton that 
triggers an electric pulse. This process, that started in January 2005, improved 
further the sensibility of the NC detection. 

In the SNO salt phase were observed 2176 ±78 CC, 2010 ±85 NC and 279 ±26 ES events. 
This corresponds to the following neutrino fluxes [27]: 

(1.68 ± 0.06;E]:E]^) X 10^ cm-^s"^ (1.203) 
{AM±0.2ltoil) X 10^ cm-^s"^ (1.204) 
(2.35 ± 0.22 ± 0.15) x 10^ cm'^s"^ (1.205) 

We see immediately that these fluxes are not compatible with each other. More specif- 
ically, the fact that the recorded flux for NC (that is sensitive to the three neutrino 
flavours) is larger than the CC one (only sensitive to electron neutrinos) indicates that 
some electron neutrinos were converted to muon or tau neutrinos in their travel from the 
core of the sun to its surface. From the experimental data the fluxes of muon and tau 
neutrinos can be estimated: 

$SNO,NC ^ (3 26 ± 0.25;°^) X lO'^ cm-^s-i (1.206) 
$SNO,ES ^ (4 3g ^ i.52;oj^) X 10*^ cm-^s-i (1.207) 

which are in good agreement with each other. Thus, the neutral current channel used in 
SNO was able to suggest a solution to the Solar Neutrino Problem, providing the first 
evidence that neutrinos oscillate and as so have mass. 



^cc 



SNO 
NC 

SNO 
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1.3.2 Atmospheric neutrinos 



Cosmic rays that hit the Earth are mainly made of protons, with some small component 
of heavy nuclei. These primary cosmic rays interact in the atmosphere and produce all the 
hadrons and mesons, in particular pions, that decay into leptons and respective neutrinos. 
These lepton are mainly muons because the ratio of the decay widths is 

Re/, = ^JI^^^ (1.208) 

_ ml{l- ml/ ml) 
mli'^- ml/ ml) 
ci: 1.28 X lO"''. 



The muons will decay afterwards into electrons and neutrinos (see fig. 1.12). Thus, it is 
estimated that the neutrino flavours ratio is of the order of 2: 

S> -I- S)- 

~ 2 (1.209) 

in the absence of neutrino oscillations and at low and moderate energies, roughly below 5 
GeV. At higher energies the production and decay of Kaons will contribute to the emission 
of neutrinos. Also, the flavour ratio cannot be measured directly in real experiments 
because what is observed is the Cerenkov light produced by charged leptons. Since electron 
and muon neutrinos have different cross-sections for the interaction with matter and the 
efficiencies and detection criteria are different for e-like and mu-like events, it is preferred 
to report the experimental data in terms of e-like and mu-like events and to reveal the 
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\ 




Figure 1.13: Symmetry in neutrino flux distribution. Taken from [3] 



anomaly through a ratio of ratios: 




(^^i-like/^e-like), 
(^^i-like/^e-like) 



data 



MC 



(1.210) 



where N^, Nfj_ are the number of events, the numerator is the measured ration and the 
denominator is the one calculated with Monte Carlo methods. 

For high energetic cosmic rays, roughly above 1 GeV, the muons may hit the ground 
before they decay and this can suppress the ratio between the fluxes of muon and electron 
neutrinos. On the other hand, low energy cosmic rays may be severely constrained by the 
geomagnetic field. Below some energies it is expected an overabundance of cosmic rays 
in the geomagnetic poles and a suppression in the equator. This asymmetry is expected 
to vanish above few GeV: in this case, if a neutrino is observed in some underground 
laboratory A (fig. 1.13) with an azimuthal angle 9^^^ it will hit another laboratory B 
with an azimuthal angle vr — Of^. So, for high energies we have 




(1.211) 



Since the flux distribution is assumed to be independent of the location, 



(1.212) 



we conclude that it has a symmetry 




(1.213) 



This means that it is expected an equality between upwards and downwards neutrino 
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Figure 1.14: Asymmetry in atmospheric muon neutrinos as computed with Monte Carlo 
methods. Taken from [9] 



fluxes. Now, upwards neutrinos observed at some underground facility have crossed the 
Earth and are exposed to matter effects. This means that an asymmetry in the downward 
and upward fluxes is an indication that interaction with matter has induced neutrino 
oscillations. This can be quantifled by the formula 

^r'°™=(^^) , (1-214) 

where U and D are respectively the up and down neutrino fluxes integrated in the ranges 
0.2 < cos6'2 < 1 and —1 < cos6'2 < —0.2 with a = e^fi. Figure 1.14 shows Monte 
Carlo calculations of the up-down asymmetry in the absence of neutrino oscillations for 
two locations, Kamioka and Soudan. Note the cut-off on the asymmetry for energies 
bigger than 3 GeV: the geomagnetic field is no longer strong enough to deviate the cosmic 
rays which results in a equality between up and down fluxes. The different signs for the 
asymmetry are due to the different geographic locations. Kamioka is near the geomagnetic 
equator and Soudan is near the geomagnetic pole. 

In 1988 Kamiokande published data from sub-GeV events that indicated clearly an 
anomaly in atmospheric neutrino fluxes. In terms of ratio of ratios the final results of the 
Kamiokande sub-GeV and multi-GeV events are [3] 

= Qmt°om ± 0-05 (1.215) 
^multi-GeV ^ 0.571^? ± 0.07 (1.216) 

These anomalies are well explained by neutrino oscilations, even though Kamiokande was 
unable to distinguish between /i o r and /i f-> e transitions. However the results by the 
CHOOZ experiment in 1997 excluded fi ^ e transitions. 
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Figure 1.15: Ratio of measured to expected Ve flux of different reactor experiments. The 
dotted curve corresponds to the best-fit values ArrigQ^ = 5.5 x 10~^ eV^ and sin^ 20goL ~ 0.83. 
Taken from [3] and [10]. 



1.3.3 Reactor neutrinos 

The fission process in nuclear reactors produces electron anti-neutrinos in abundance from 
/3-decay in neutron-ricli nucleus. However the energy of these anti-neutrinos is small, in 
the order of a few MeV and this excludes the possibility of detecting fi and r neutrinos 
from charged current interactions. So in these experiments only the disappearance of 
z7g can be investigated. Figure 1.15 shows the ratio of expected to observed flux for 
various experiments. The experiments with L ~ 10 — 100 m are short-baseline (SBL), the 
ones with L ~ 1 Km are long-baseline (LBL) and finally the KamLAND experiment with 
L ~ 200 km is very- long-baseline (VLBL). It is clear that only KamLAND found evidence 
of neutrino oscillations. However the absence of a signal for the other experiments leads 
to exclusions curves in the Am|QL , sin^ 2^sol plane and to upper bounds for AttiIql . 
Recall that the sensitive to Am^Q^ is determined by the value oi x/E for which we have 

Am^x 

(1-217) 

It is now established that AmgQL is of order 10~^ eV^ and for an energy of few MeV this 
requires a source-detector distance around 100 km in order to get an evidence of neutrino 
oscillations, as it was the case with the KamLAND experiment. 

CHOOZ Experiment: The CHOOZ experiment was located near the CHOOZ power 
plant in France which is composed by two nuclear reactors. The distance between the 
detector and the two reactors was 998 m and 1115 m which, by (1.217), means that the 
experiment was sensible to Am^ of the order of xlO""^ eV^. The CHOOZ experiment 
gathered data from April 1997 to July 1998 and the ratio of measured to expected Ue 
events was [3] 

i? = 1.01 ±0.028 ±0.027 (1.218) 
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Figure 1.16: Exclusion curves for tlie CHOOZ experiment. Also shown is the allowed region 
for —7- Ve from the Kamiokande experiment. Taken from [3] and [11]. 



The exclusion curve is shown in figure 1.16 together with the allowed range from Kamiokande 
for — )> Ue transitions as the solution for the atmospheric neutrino deficit. The data im- 
plies that [3] 

sin2 26l<0.1, for Am^ > 2 X 10"^ eV^ (1.219) 
Am^ < 7 X 10"^ eV, for sin^ 20 = 1. (1.220) 

and excludes the Kamiokande solution for the atmospheric neutrino anomaly. 

KamLAND Experiment: this experiment [3] was designed to measure Ve fluxes 
from nuclear reactors of 53 power plants in Japan. The detector is located in the Kamioka 
mine where the Kamiokande experiment was previously installed. It consists of 1200 m"^ 
of liquid scintillator that occupies a spherical balloon with a diameter of 13 m. Neutrinos 
are detected through the inverse j3 decay 

Ue+P — >e+ + n (1.221) 

where the delayed coincidence from the light signals emitted by the positron and the 
neutron capture allows an efficient reduction of the background. The ratio of measured 
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Figure 1.17: (a) Allowed region in the tan^6',Am^ plane obtained in the KamLAND exper- 
iment. The lines show the regions allowed by solar neutrino data, (b) Combined analysis of 
two-neutrino oscillations with KamLAND and solar neutrino data. Taken from [3] and [12]. 



to expected Ue events from March 2002 to January 2004 is [3] 

= 0.658 ± 0.044 ± 0.047 (1.222) 

which deviates from unity by 5a. The best fit of the data is obtained with 

Am^ = 7.9t°oi X 10-5 eV^ (1.223) 

Figure 1.17 shows the data in tlie Am^,tan6' plane and it is clear that there is a big 
uncertainty in the mixing angle. But a global analysis with the solar data allows a much 
better resolution for 9. 



1.3.4 Accelerator neutrinos 

These experiments detect muon neutrinos from the decays of muon and kaons produced 
by the collision of a proton beam with a target. We will discuss only the K2K experiment. 
K2K Experiment: the accelerator is located in tlie KEK laboratory in Japan and the 
neutrino flux is detected at the Superkamiokande experiment. The two facilities dist 
around 250 km. The neutrino beam is an almost pure beam with roughly 1% of muon 
anti-neutrinos and 1% of electron neutrinos. The events for K2K at Superkamiokande are 
selected using GPS synchronization with the proton beam at KEK and the initial neutrino 
flux is monitored by a near detector similar to the one of Superkamiokande. K2K has 
been divided in two phases. The first one, K2K-I, lasted from June 1999 to July 2001 
and the second one, K2K-II, started in January 2003 and ended in February 2004. In 
total the experiment observed 107 fully contained yU— like events with an expected number 
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Figure 1.18: Number of events as a function of tlie reconstructed neutrino energy. Points witli 
error bars correspond to data, the solid curve corresponds to best fit with oscillations and the 
dashed line refers to no oscillations. Taken from [3] and [13]. 

without oscillations of 15ll}o- Figure 1.18 shows the energy distribution of K2K events 
plotted with the best fit curve with oscillations (solid) and the curve without oscillations 
(dashed). The best fit values are [3] 




sin^ 29 = 1.0 Am^ = 2.8 x 10"^ eV^. (1.224) 
1.3.5 Neutrino parameters 

In [28] a global analysis of neutrino oscillations was done, taking into account the latest 
data at the time of KamLAND and K2K, as well as state-of-the-art solar and atmospheric 
neutrino fluxes. In [14] the discussion was updated taking also into account the data 
released in 2008 by the MINOS collaboration, the Sudbury Neutrino Observatory (SNO), 
KamLAND and Borexino. We present here the major results of these references. In figure 
1.19 the curves for solar, Kamland and global analysis are shown. For this result new 
data from the SNO experiment was taken into account, in particular the data from CO 
and NC fiuxes, that quote a lower value for the flux ratio, (pcc/'pNC = 0.301 ±0.033 which 
leads to a stronger upper bound on sin^0i2. Also included are direct measurements of 
the solar neutrino rate performed by the Borexino collaboration. For the global study 
a recent re-analysis of the Gallex data was also taken into account. More details can be 
found on [14]. At la the best fit points are [14] 

sin^ 9^2 = 0.304l°:°?2 , Am^i = T.GSlj]:! x 10"^ eV^ . (1.225) 

For the atmospheric neutrino parameters, long-baseline accelerator and atmospheric neu- 
trino measurements from Superkamiokande were combined. The results are shown in 
figure 1.20. In [14] were found the following best fit points for atmospheric parameters, 
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Figure 1.19: Determination of tlie leading solar oscillation parameters from artificial and 
natural neutrino sources. Curves for data from solar, KamLAND and global analysis are shown. 
Taken from [14]. 



at la: 

sin^ 023 = O.bOtZl > \Aml,\ = 2.401°;}? x IQ-^eV^ . (1.226) 

We note that the data is consistent with maximal mixing 623 = it /A and both normal 
^ mi ~ m2 or inverted <C mi ~ 1712 hierarchy. The resolution of this last 
ambiguity lies in a direct determination of neutrino masses, as oscillation experiments are 
only sensitive to the squares of mass differences. As for ^13 we see from (1.142), (1.144) 
and (1.146) that a non-zero value for this mixing angle is a necessary condition to observe 
CP violation. Figure 1.21 on the left summarizes the information on ^13 from the present 
data, where the different experiments are complementary and all contribute to give a 
more rigorous bound. At 90% confidence level {3a) [14] reports the following limits: 

{ 0.060 (0.089) (solar +KamL AND) 
0.027 (0.058) (CH00Z+atm+K2K+MIN0S) (1.227) 
0.035 (0.056) (global data) 

The interplay of solar and KamLAND is shown on the right of figure 1.21 where the 
slightly smaller ration CC/NC for SNO mentioned above gives a hint for a non-zero value 
of 6'i3, even though is still compatible with a null value for this mixing angle. Details can 
be found on [14] . 

In summary, we see that neutrino data is consistent with the Tri-Bi-Maximal mixing 
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Figure 1.20: Determination of tlie leading atmospheric oscillation parameters from artificial 
and natural neutrino sources. Curves for data from atmospheric, MINOS and global analysis 
are shown. The dot, star and diamons indicate the best fit points of atmospheric, MINOS and 
global data, respectively. Taken from [14]. 



simplifying hypothesis [29] that leads to the HPS mixing matrix U (see 1.145): 



U 



2 
3 
1 

Ve 
1 



4. 



1 

V3 

1 
V3 

1 
V3 



\ 



1 

"73 / 



(1.228) 



In Chapters 5 and 6 we will assume that the leptonic mixing matrix has this form. 
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sin sin^e^2 

Figure 1.21: Left: constraints on sin'^ ^13 from different parts of the global data. Right: 
allowed regions in the {9i2 — O13) plane at 90 and 99.73% CL as well as the 99.73% region for the 
combined analysis. The dot, star and diamond indicate the best fit points of solar, KamLAND 
and combined data, respectively. Taken from [14]. 



Chapter 2 

Baryon Asymmetry, Leptogenesis and Dark 
Matter 

2.1 Anomalies, Baryonic and Leptonic Currents 

Anomalies play a central role in the Standard Model, where its absence or total cancelation 
is essential to the renormalizability of the theory. The anomaly associated with the chiral 
gauge transformations is related with the non invariance of the fermionic measure in the 
path integral formulation of the theory. This was first pointed out by Fujikawa [30] and 
it is in connection with a deep result in Mathematics, the Atiyah-Singer theorem [31], as 
it is shown in [1]. 

The explicit computation of the divergence of the gauge 3- vertex function with fermion 
loops (fig. 2.1) shows that the anomaly is given by 

[5.4la„on. = -^2D^P'i^'''"'FLFl. (2.1) 

where F^^ is the non-Abelian field strenght associated with the gauge bosons A'^ that 
belong to the adjoint representation Tq, of the gauge group, Dap^ is the totally symmetric 
quantity 

D^p, = h:r{{T^,Tp}T,) (2.2) 
and is the fermionic classically conserved current 

J^ = Xl'T^X- (2.3) 

Here x is a fermionic field that unifies all left-handed fermions and anti-fermions that 
transform nontrivially under the representations of the gauge group. 

There is always some freedom in choosing which of the 3 fermionic currents on the 
triangle diagram is anomalous, even though it is not possible to remove simultaneously 
the anomaly. In this way, we could as well have written the anomaly (2.1) for or J^. 
This is particularly useful when one has a global symmetry that it is not gauged, coupled 
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Figure 2.1: One loop gauge 3- vertex that contributes to the triangle anomaly. 
Table 2.1: Standard Model quantum numbers 



Particle 



ULi 

+" 

,0 



gluons 



9 



B 



L 


5f/(3) 


SU{2) 


U{1) 





3 


2 


1 

3 





3 


1 


_4 





3 


1 


3 


1 


1 


2 


-1 


-1 


1 


1 


+2 





1 


3 








1 


2 


1 





3 


1 






to 2 gauge fields. In this case, one can choose the anomaly to be in the non-gauged 
current. That is what happens with the baryon number current, 



Jb = BXeYXb 



(2.4) 



where B is the baryon number, which is conserved at the classical level but gets anomalous 
with the triangle diagram (fig. 2.2). To compute the anomaly, we must calculate the 
coefficient (2.2) for all the possible couplings with the gauge fields of the Standard Model. 
It is easy to see that only pairs of fields in the same gauge group contribute, since the 
traces of the SU{?>) and SU{2) generators are zero. The relevant quantum numbers are 
shown in table 2.1 (see also sec. 2.3). The normalization condition adopted here for 
SU{N) generators is 

Nr 

Tr {T^T^) = ^6^f,, (2.5) 



2.1 Anomalies, Baryonic and Leptonic Currents 
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Figure 2.2: Anomaly in th.e baryonic current. 

where Nq is a constant that depends on the representation of the group G = SU (N). For 
the baryonic current coupled to 2 SU (3) fields there is no anomaly, because 

= 0. (2.6) 

3,3 

In the case of SU{2), the coefficient Dais^y is 

= ^^2, (2.7) 

2 

for one generation, where a color factor of 3 was taken into account and N2 is the factor 
Nq above for SU{2), and for t/(l), the same coefficient is 

^BY'^ = -2. (2.8) 

This shows that, in the Standard Model, the baryonic current has a total anomaly 

with Ng the number of generations. 

In exactly the same way one can compute the anomaly for the leptonic current and 
find that it is equal, assuming the same number of generations for quarks and leptons, as 
it is required by the cancelation of gauge anomalies in the Standard Model. So, even if the 
leptonic and baryonic currents are anomalous, their difference Jg — is conserved, which 
makes B — L a, good quantum number. This may have serious cosmological implications 
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in the way the asymmetry between matter and anti-matter is obtained, as we will see. 
We finish noting that the differences between lepton flavour currents are also conserved, 
so that Li — Lj are also good quantum numbers. 



2.2 Instantons and Sphalerons 

The variation of macroscopic charge Q associated with a given classically conserved 
Noether current is usually given by 



Q = J d-^xd,.r = 0. (2.10) 
In the case of the anomalous baryonic current, the variation of baryon number is 

B = ^.e-''" J d'x {f,^,F,,^ - ^Fl^F^^)i , (2.11) 



so that 



Bf = Bo + ^e"^^'' / d'x I^F,^,F,,^ - ^F,%Fi^}l . (2.12) 



Now, it can be seen that the current 

Ki = eujk^Tr (^AF^k + jAAjA,^ , (2.13) 
where we have changed to Euclidean space through a Wick rotation, gives 

dm^^m — '^^ijklFotijFakli (2'14) 

SO that the variation of baryon number can be given as a surface integral of a total baryon 
current density. Assuming that the field strengths and F^ij vanish at infinity, the U{1) 
gauge field does not contribute and then 

with a three dimensional sphere with a radius approaching infinity. Naively, we would 
be tempted to consider that this integral is zero for gauge fields vanishing sufficiently fast. 
However, there are configurations for A" that do not vanish at infinity, but give there 
zero field strength. These are known as Instantons, after Belavin, Polyakov, Schwarz and 
Tyupkin [32] and have asymptotically a pure gauge configuration: 

iA, g{x)-^dig{x), (2.16) 
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where Ai = Aft"' and g{x) is a group element of SU{2). Belavin et al. [32] found a 
solution: ^ 

^^'(^) = (nTi^^) (2.17) 

where, explicitly, 

= (£i±|i:l) , (2.18) 

and t are the usual SU{2) generators. This allows us to define the quantity 

= / mW,-Tr (2.19) 

in such a way that the variation of baryon number is 

AB = ^X. (2.20) 

The quantity I is well known by mathematicians and it is called the Maurer-Cartan form. 
It classifies homotopy classes of maps^ of 5^ into a Lie group. A computation of I can 
be found in [1] , where it is shown that 

X = 247r^ (2.21) 

because g(x) belongs to the homotopy class of the identity. Since the homotopy group 
7^3(0) of a compact connected simple Lie group G that classifies homotopy classes of 
maps from into G is^ Z, a general homotopy class is identified by an integer u called 
the winding number and can be obtained by superimposing several solutions that have 
winding number equal to 1 or -1. So, in general, 

I = 247rV. (2.22) 

This shows that Instantons change the baryon number by an integral value: 

AB = NgU. (2.23) 

In virtue of the inequality 



< - / ( ^ -e,,MFo,M ] d^x, (2.24) 



^An homotopy class of maps on a Lie group is a set of maps that can be continuously deformed into 
each other. 

^5''^ is homeomorphic to SU{2), so this is the same as considering maps from SU{2) into G. Moreover, 
it is possible to restrict these maps to a SU{2) subgroup of G 
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the Euclidean Yang-Mills action 



1 



will satisfy the inequality 



4g. 



^ — A „2 I FaijFaijd X, (2.25) 
"2 



(2.26) 



Because the field strengths from (2.17) are self-dual 



■^i - 2''^jkiFki, (2.27) 
equation (2.26) is in fact an equality for Instantons, which allows us to write 

5 = Stt^ |z/[ /52- (2.28) 

In the Standard Model, g2 = e/sinOw and, for e^/Air ~ 1/129 and sin^ 9w — 0.23 
evaluated at Mz^ equation (2.28) gives a suppression factor of exp(— 186), which indicates 
that Instanton processes are unlikely to be observed nowadays. However, the potential 
barrier could be surpassed if thermal fluctuations are important, as it could be the case 
in the early Universe. 

Until now the discussion has been somehow oversimplified because we are considering 
a pure gauge theory. The inclusion of matter and Higgs fields can change the picture 
in a significant way. This issue was first addressed in two papers by Klinkhammer [33] 
and Manton [33] [34] , where the name Sphaleron appeared to describe a saddle point for 
the energy on the Weinberg-Salam model configuration space. Here, we present a short 
summary, following [35]. In the limit of vanishing mixing angle sin 9w 0, the Sphaleron 
has the following form, in the Aq = gauge: 

A =fi92vr)U^d,iU°^)-' (2.29) 
$ =^h{g,vr)U^ ^ j , (2.30) 

with U°° = g{xo = 0, x), where g{x) is given by (2.18), and the functions / and h have 
the asymptotic behavior 

no ^ / , (2-31) 
[1, e^oo 

hiO ^ , (2.32) 

1, e ^ oo 
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The Euclidean energy functional 



E = J d^x -^F^jF^^ + (A$)^(A^) + A (^^'^ - i 



(2.33) 



can be reexpressed as 



E = — / de 

92 Jo 



+l[/(i-/)jni^(|) + 



(2.34) 



The Sphaleron solution (2.29) and (2.30) interpolates between a vacuum state for ^ — )> cx3 
and the maximum of the Higgs potential ($° = 0) for ^ — )■ 0. The vacuum solution for 
the gauge fields is the same as the Instanton with winding number 1, and superimposing 
the solution n times we get a solution with winding number n. So, the Sphaleron 
is a saddle point in field configuration space, the lowest barrier between two topological 
distinct vacua and transitions from one vacuum state to another will induce violations 
of baryon and lepton number. Using dimensional analysis, one estimates roughly the 
Sphaleron energy: 



An 



(2.35) 
(2.36) 



where / is a typical length scale for the Sphaleron solution. Minimizing E{Aj) + E{^), 
one gets 

(2.37) 



1 



/,„ ~ ~ 10"^'^ cm 

g2V 



and 



Ap ~ ~ 10 TeV, 

92 



(2.38) 



which sets the order of magnitude for the Sphaleron potential barrier. A more accurate 
result was found in [33] by means of variational methods: 



E^ 



sp 



92 92 



(2.39) 



where i? is a function that depends weakly on X/g2- -S(O) ~ 1.52 and i?(oo) ~ 2.72. The 
effects of a mixing angle 9u, different from zero are small. 



In [36] the rate of Sphaleron processes was computed: 

2Mw{T) 



sph 



TC exp 



T92{T) 



B{\{T)lg,{T)) 



(2.40) 
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where C is a constant. This should be compared with the expansion rate of the Universe 



where Ngjf ~ 100 are the effective relativistic degrees of freedom of the Standard Model 
and rupi is the Planck Mass. One sees that the Sphaleron rate exceeds H for T > T* (see 
the comments in [36]), where 



This will depend on the value of A(T'), but for \{T) ~ §2, one finds that B = 2.1 and 
T* ~ 0.6Tc, where is the temperature for the electroweak symmetry breaking transition. 
Since Sphalerons are only active below T^, we see that there is a temperature range where 
they are in thermal equilibrium. On the other hand, the trivial relation 



and the fact that B — L is a, good quantum number, suggest that any previous baryon 
asymmetry generated by any process that gives B — L = will be erased by Sphalerons 
when in thermal equilibrium, since they violate B + L. However, if some process creates a 
non-zero value of i? — L at a high scale, for instance, by means of the decaying of a heavy 
particle into leptons or anti-leptons, then it could survive through the period of thermal 
equilibrium, even though the values of B and L could change. This was first proposed by 
Fukugita and Yanagida [37] and it is the basis of Leptogenesis. 

2.3 Baryogenesis and Leptogenesis 

In order to estimate the baryon number that survives the period of thermal equilibrium, 
one must express it as a function of the conserved quantum numbers. This was done in 
[38], but the argument presented here is from [39]. 

From a set of conserved quantities Qa^ like isospin, hypercharge or B — L and their 
quantum numbers qai for each particle specie i, we can express the chemical potentials as 



(2.41) 



T* = {2Mw{T)/g2{T)ln{m,i/T*))B{X/g2). 



(2.42) 



B = -{B + L) + -{B-L) 



(2.43) 




(2.44) 



a 



with some linear coefficients fia- The relativistic particle density is 




(2.45) 



(2.46) 
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where gi are the effective degrees of freedom for the particle specie i. The antiparticle 
density will be given by the same formula, but with fii replaced by —fii, so the difference 
is 

n,-W, = 8,,.(^)'ff (2.47) 

assuming that the imbalance between particles and antiparticles is small, that is, \ <C 1. 
The integral has the value for bosons and n^/Q for fermions. This means that 

ni-ni = —Qiiiu (2.48) 
6 

where 'gi are the effective degrees of freedom, with an additional factor of 2 for bosons. 
So, from equation (2.44) we have 

ni-ni = —Qi'^qaiiia- (2.49) 

a 

The density of the conserved quantity Qa is 

ria = ^^qai{ni - rii) = —"^MabfJ^b, (2.50) 

i b 

where M is the positive definite matrix 

Mab = ^MaiQbi- (2.51) 

i 

Since M has an inverse, we can invert equation (2.50) and use this in equation (2.49) to 
obtain 

ni-ni = ^giqatM'j^rib. (2.52) 

ah 

From this, we can obtain the baryon number density as 

Ub = ^Biirii - Ui). (2.53) 

i 

It is a simple exercise to compute the matrix elements Mab from table 2.1: 

Mb-l,b-l = (2.54) 
Mb-l,y = (2.55) 
My,y = ^ + N,, (2.56) 



58 



Baryon Asymmetry, Leptogenesis and Dark Matter 



where is the number of scalar doublets, from which we can obtain M^j, : 

MsUr = (2-58) 



13iV 
~3D 



M-\ = (2.59) 



with the determinant 



D = ^ + (2.60) 



22N' UNgNd 

This allows to compute the baryon number density from equation (2.53): 

riB = Y,9^B^ {{B - L),M-\^_^ + Y,My,b-l) Ub-l (2.61) 

i 

( 8iVg + 47V, ^ 

= U2iV, + 13ivJ "---- ^2-^2) 

taking note that the hypercharge density is zero, from gauge invariance. For Ng = 3 and 
Nd = 1 this gives Hb = {28/79)nB-L- In any case, Ub is of the same order of magnitude 
as riB-Li as it would be expected from the naive relation (2.43). 



2.4 Cosmology and Dark Matter 

In this section we review briefly the Standard Cosmological Model, based on the Firedmann- 
Robertson- Walker metric, then we discuss in detail the Boltzmann equation that describes 
the out of equilibrium processes in the early Universe and we finish with the applications 
of this equation to the cosmological bound on neutrino masses and to the mass estimate 
of WIMPs as the major constituents of Dark Matter. Major references are [3, 40-42]. 



2.4.1 The Standard Cosmological Model 

Modern Cosmology is based on the General Theory of Gravitation, first proposed by 
Albert Einstein in 1915. The dynamics of space-time are described by Einstein's equations 

R'"' - ^Rg'"' = SttGnT'"' + Agf"" (2.63) 

where g'^'^ is the space-time metric, R^" is the Ricci tensor, R the scalar curvature, T^'" 
the energy-momentum tensor and A the cosmological constant. The Ricci tensor and the 
scalar curvature are obtained from the full curvature tensor 
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with the Christoffel symbols being given by the metric 



More specifically we have 



and 



R^"" = gp^R^P'"' (2.66) 



R = gp.R^". (2.67) 
The free falling motion is determined by the geodesies equation 

with the proper time given by 

dr^ = g^fjdx^dx^. (2.69) 

The energy momentum tensor describes the matter content of the model. For a perfect 
fluid we have 

T'"' = {p + p) u^v!" - pg^"" (2.70) 
with u the proper velocity, p the energy density and p the pressure. 



Observations of the large scale structure of the Universe suggest that it is spatially 
homogeneous and isotropic at distances roughly bigger than 10 Mpc. A solution for 
Einstein equations for such an Universe is given by the Friedmann-Robertson- Walker 
metric 

dT^ = de - R{tf I + r\de^ + sin' edcp')) (2.71) 
y 1 — kr^ J 

which describes a perfect fluid with the energy momentum tensor given by (2.70) in a 
co-moving frame. The scale factor describes the expansion of the Universe and it has the 
dimension of length; k is the spacial curvature normalized to 

-1 open Universe 
k = { flat Universe (2.72) 
1 closed Universe 

The fact that R{t) depends on time implies an apparent cosmological doppler effect known 
as redshift, even though its origin is in the expansion of spacetime. In fact, since photons 
are obviously light-like (dr = 0), from (2.71) we have for a particular wave front 

*° dt _ n dr 

W)~Jr 7r^' ^ ' 



60 



Baryon Asymmetry, Leptogenesis and Dark Matter 



where (t,r) are relative to the emission and (to,^o) are relative to the reception. If 
we consider the next wave front differing from the previous in space and time by the 
wavelength and the period we have that 



so that 



which means that 



dt _ r° dr _ dt 

t+T 

_dt dt 

to 



R{t) J, R{t) 



(2.74) 



(2.75) 



An A 



R{to) R{t) 
Aq -R(io 



where z is the redshift: 



A R{t) 

= l + z, (2.76) 

z = (2.77) 
A 

The redshift can be directly measured observing the electromagnetic spectrum from stars, 
nebulae and galaxies. Observational data shows irrefutably that z > 1 and as so the 
Universe is expanding. 

When the FRW metric is used in Einstein's equations with the perfect fluid energy- 
momentum tensor, we obtain Friedmann equations 

= (2.78) 

with the Hubble parameter, given by 

H = ^. (2.79) 



R 

This definition has its origin in the expansion 

1 _ R{t) 



1 + 2 i?(t, 



(2.80) 



0) 



l + f||y(t-io) + 0((t-to)') (2.81) 



which can be seen to give the Hubble law for z <C 1, 



z = HodL 



(2.82) 
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where is the luminosity distance. 

From the Friedmann equation we see that for a flat Universe {k = 0) the density is 
equal to the critical density 

Pc = (2.83) 

Its present value is [3] 

° - (10.5369 ± 0.0016) keV cm"^ (2.84) 



87rG 



N 



where the present value of the Hubble parameter is Hq = lOO/i km s ^Mpc ^ . The Particle 
Data Group [43] quotes h = 0.72(3). 

The energy density is usually expressed in terms of a relative density 

VL = —, (2.85) 

Pc 

which allows to write Friedmann equation as 

- 1 ~. (2.86) 



Matter, radiation and vacuum have different evolutions. Matter is understood as 
non-relativistic particles with vanishing pressure. Radiation obeys the equation of state 

PR = \pR (2.87) 

and vacuum has a negative pressure which is a consequence of the energy conservation 
equation. As the Universe evolves matter density is suppressed by a factor 

Pm oc R-^ oc (1 + zf (2.88) 

because the volume expands by a factor . As for the energy density, it is suppressed 
by an additional factor of R~^ that takes into account the redshift 

PR oc R-'^ oc (1 + zf. (2.89) 

Finally the vacuum energy density is constant. These different evolution rates suggest 
that the Universe has three distinct epochs: radiation dominated, its early epoch just 
after inflation; matter dominated, an intermediate stage; vacuum dominated, its latest 
era. It will be instructive to obtain an estimate of the Universe age, taking into account 
these different contributions by radiation, matter and vacuum. Friedmann equation (2.86) 
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allows to obtain the curvature constant as a function of current parameters 

k = {no- 1) H^Rl (2.90) 
and replacing back in (2.86) we get 

= H!{l + zr^. (2.91) 

or 2 

f^- 1 = ^(1 + ^)'(^^o-l). (2.92) 

Formula (2.91) gives the evolution of the Hubble parameter in terms of the relative density 
and the redshift and it will be useful in a moment. The last formula can give the relative 
density as a function of present values and the redshift z: combining the evolutions of the 
three states 

Pom I -K 



Pm \Ro 

POA 

Pa 

POR f 



= (1 + z)-' (2.93) 
1 (2.94) 



— =(^ =(l + z)- (2.95) 

PR \Ro, 



with the observation 



Pc 



fn, ) , (2.96) 



we can get the three relative densities as 



n« = fiofl(l + 2)'(^) (2.97) 
= noM{l + zf (^^^ (2.98) 
^^A = f^oA (^) (2.99) 
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and the total Q, as 

= ((1 + zfnoM + (1 + 2)'0oij + Qoa) • (2.100) 

Obtaining the Hubble parameters ratio from this equation and replacing in (2.92) we get 
after some trivial algebra 

n-l = ^° ~ ^ (2 1011 

1 - f^o + (1 + z)^oM + (1 + z)moR + (1 + zy^QoA ■ ^ ' ' 

This in turn allows to express the Hubble parameter solely in terms of present values and 
the redshift. Replacing in (2.91) we have 

H = Ho{l + 2) Vl - 5^0 + (1 + z)noM + (1 + z)^noR + (1 + z)-2(loA. (2.102) 

From the definition of redshift (2.76) we can relate the variation of time with the variation 
of z 

1 dz , , 

* = -Ji—, (2-103) 

SO that finally we can integrate with equation (2.102) to obtain 

1 dx , , 

to = — ^ (2.104) 

with X = {1 + z)~^. Neglecting the radiation contribution and for A; = the integral can 
be performed analytically to give 



t, = 2^-1^ In (2.105) 



3 -s/HoA ^/l — 



OA 



The energy density for a dilute, weakly-interacting gas of particles with g internal 
degrees of freedom is given in terms of its distribution function /(p): 

P=7A-, I E{p)f{p)d'p. (2.106) 



(2^) 

If the specie is in thermal equilibrium then the phase space distribution is given by the 
Fermi-Dirac or Bose-Einstein distribution 

^f"' - e.p((^(p)-,)/T±l) 

with +1 for fermions and —1 for bosons and /i the chemical potential. Moreover, if the 
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particles are relativistic then this expression can be easily integrated to give 



^gT^, bosons 
^gT'^, fermions 



30 

IZL 

8 30 



(2.108) 



Also, for relativistic particles 



P = gP- (2.109) 



Replacing the relativistic energy density in Friedmann equation with = gives the 
evolution of the Hubble parameter as a function of the temperature in the radiation 
dominated era: 



where g^ = (1,7/8) (bosons, fermions). Also, from the fact that p decreases as 
Friedmann equation shows once again that 

E = — (2.111) 
which allows to obtain a relation between time and temperature: 

1 / 90 \ Mpi 



with Mpi the Planck mass in natural units. This will be useful in the discussion of 
Boltzmann equation. 



2.4.2 Boltzmann equation 

Following closely [40] and [41] we now turn to the derivation of Boltzmann equation 
describing the out-of equilibrium processes in the early Universe. 
For an affine parameter r we have 

dT^ = g^^dzf^dx", (2.113) 

such that the four- velocity has norm 1: 

g^^u^u" = 1. (2.114) 

It is convenient to reparametrize r [39] as A = r/m where m is the particle mass in a way 
that 

m^dX^ = g^^dxf'dx" (2.115) 
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(2.116) 



and that the four-momentum p'' is 

= m— — 

(It 

In this case the geodesies equation (2.68) becomes 

^ + r>''p- = 0. (2.117) 
From the FRW metric (2.71) we can obtain the evolution equations for p^: 

^ = -p^i?^ (2.118) 

To get the evolution of some scalar density function / with variables x{X) and p(A), we 
note that 

df _ dx^' df ^ dpf" df 



dX dX dxi' dX dpt" 



The variation in can be dropped on the basis of the FRW metric isotropy. Also, we 
are interested in the variation on mass-shell for a spatially fiat metric. In this way we 
consider the quantity 

f{p,t) = J f{p',p,t)s{po-{p'R' + m'f^yp'. (2.121) 

Then, integrating (2.120), we have 

^^^j{po-{p'R' + m^y^)dp^ 



l^6{po-ip'R' + m')^)dp'^ 



-2^p^^ 

dp' 

+ ^RR^ 

pO QpO 



30 = (p2i?2+m2)l/2 



(2.122) 



Observe that in obtaining this formula we must differentiate the delta function with re- 
spect top' (see below). That is why the last sign is symmetric from (2.120). Differentiating 
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(2.121) with respect to t 

df{p,t) [df{p\v.t) 



dt 



dt 



6{po-{p^R^ + m')^)dp' 



+ J f{p',P,t)p{po-{p'R' + m'f^)dp'. (2.123) 



Noting that 



fip^p,t)p{po-{p'R' + m')^)dp' 
dt 27r 



f{p\p,t)RR- 



p2i?2 + ^2-^2 27r 



/ , ,,. -^J{p\^.t)^ (Po - {p'R' + m')^) dp' 

J {p^R^ + m^)2 dp^ V / 



^RR^ 

pO QpO 



p0 = (p2fl2+m2)l/2 



we see that 



dfip,t) _ fdfip',p,t) 
dt 



dt 

+ ^RR^ 

pO QpO 



S(po- (p"/?" + m')5 dp 



p" = (p2_R2+,n2)l/2 



which allows to write (2.122) as 



d^l 
dXp' 



R^ dpi 



(po - {p'R' + m')^ ] dp" = ::^L^£i^ - 2'-p 



dt 



(2.124) 



(2.125) 



Dropping the tildes, using the isotropy of FRW metric and putting p = |p| we can define 
tlie operator 

(2.126) 



dt ^R^dp 



To see that this is the pertinent operator in studying processes with a variable number 
of particles, we start by defining the local momentum [41] p = Rp. Using the zero 
component of (2.116) it is easy to show that 



(2.127) 

and repeating the computation that led to (2.126) gives the operator L{f) in terms of 
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local momentum: 



From now on we will drop the bars and consider only local momentum. The particle 
current density is 



where / is the momentum distribution function and g the effective degrees of freedom. 
Because of spatial symmetry only the zero component does not vanish, which gives the 
particle number density 

As it is well known the covariant divergence of a four- vector is 

= ^2^,(9'^'^'), (2.131) 

where here and only here g denotes the metric determinant. So, for the flat FRW metric, 
9 d f ^ f ^ d^p 



R^dt\ J ' (27r)3 

R f d^p f df d^p 



R J ' (27r)3 J dt. (27r) 

3^ h r ^fi 



^ R f r d'p ^ R f d d'p f d'p 

= '^R J ^ (2^0^ ^'rJ P'w/Wr ^ ' J W ^ 

where we used (2.128). Integrating by parts cancels the first and the second terms, so we 
end up with 

d^p 9 f r>3 f r d^p 



The operator L is the relativistic Liouville operator and a necessary and sufficient condi- 
tion for the conservation of the number of particles is L{f) = 0. Boltzmann equation can 
be written in a generic way as 

L(/) = C{E)/E, (2.134) 

where C{E) is the collision term and encodes all the specificities of the processes that 
change the number of particles. The evolution of the particle density can then be given 
in the more familiar form 



In particular, C{E) contains the scattering matrix for some specific process. Consider as 
an example the process that destroys the particles i^j and creates the particles k^l. Then 
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the amplitude A for this process must be proportional to the creation operators ol k^l 
and anihilation operators of j 

A oc cilalaiaj. (2.136) 
On the other hand, the normalization of creation and annihilation operators is 

4 \n1n2 . ..rik - ..) = Vl ± rik \n1n2 ...Uk + l...) (2.137) 
tti \n1n2 .. .Hi. . .) = \n1n2 . . .Hi - 1. . .) . (2.138) 

To compute the T matrix, which is the nontrivial scattering matrix after factoring out the 
conservation delta functions, we must take the square of the amplitude and this means 
that for this process 

Toe \{kl\A\ij)\'^ oc {l±nk){l±ni)ninj. (2.139) 

Since we are integrating over the phase-space density what counts is in fact the distribution 
densities fi, fj,--- related with the number density by (2.130). Following [40] we can write 
the Boltzmann equation for the specie i 

+ 3Hni = - 9i j ^^:;;^F^^nj'inA;(in, X {271)^6 {pk +Pi-Pt- Pj) 

± h){i ± ft) - \A\i^^^^^^ /,/,(! ± mi ± m . 

(2.140) 

where dH is the relativistic phase space density 

dn = — (2.141) 
(27r)3 2p0 ^ ' 

The minus sign in (2.140) reflects the fact that if the forward scattering prevails with 
respect to the inverse one, then obviously rii will decrease. Of course, identical equations 
must be introduced to take into account the evolution of the other species. This means 
that often we have a set of coupled Boltzmann equations that can only be solved nu- 
merically. But in some limit cases it is possible to assume that only one specie departs 
significantly from equilibrium and approximate analytical solutions are possible. When 
the i specie is heavy and interacts with some relativistic particles, it is possible to obtain 
an approximate expression for the Boltzmann equation that occurs often in the literature, 
first stated without proof in [44] and derived in [42] . We will follow closely this last refer- 
ence. Assuming a two body inelastic scattering between two heavy fermion L, L and two 
light ones /, I 

I{pr)+J{p2)^ Lip[) + L{p,) (2.142) 
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with Fermi-Dirac distribution functions / and g, Boltzmann equation assumes the form 

X 1^1' (5(pi,t)5(p2,t)(l - - fip'2,t)) 

-fip[,t)fip,,t)il-gip,,t))il-gip2,t))y (2.143) 

with the CP-invariance hypothesis. Here n is the particle number density of L and L. 
The light fermions are kept in thermal equilibrium by other interactions that are very fast 
in comparison with the expansion rate so that its distribution function is 

with p^ = |p|. As for the heavy fermion, the distribution function can be approximated 
[42] by 

exp {a{t) + E{p )/T(t)) + 1 
where a{t) has the meaning of a chemical potential. With (2.144) and (2.145) note that 

(1 - f{p„t)){l - f{p',,t)) = exp(2a(t))exp (^{E[ + E',)/T{t)') f{.Pi,t)f{p^,t). (2.146) 

Energy conservation implies that 

E[ + E'^ = Ei + E2 (2.147) 

and this allows to write 

(1 - /(p'l, t))(l - /(p;, t))g{pr, t)g{p2, t) = 

exp(2a(t))(l - g{p„ t))(l - g{p2, t))f{p[, t)f{p„ t). (2.148) 

In this way Boltzmann equation (2.143) becomes 

^ + 3Hn =2 I ^^El_dUm2dU', x {27r)^5 (pi + P2 - p[ - P2) 

X \Af f{p[, t)f{p'^, t){l - 5(pi, t))(l - 5(P2, t)) (exp(2Q(t)) - 1) (2.149) 



In [42] it is argued that the chemical potential a{t) only departs significantly from zero 
for temperatures well below the mass of the heavy lepton. In this case the Fermi-Dirac 
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n, = 2 I -^.-^(py^. (2.152) 



distribution can be approximated by the Maxwell-Boltzmann distribution 

f{p,t)=e^p(^-a{t)-^y (2.150) 

This means that the actual density n{t) is related with the equilibrium one no by 

n{t) = e-"(*'no (2.151) 

with 

"° = '/(2>r)3- 

Replacing (2.150) in (2.149) and remembering (2.151) we obtain the promised result 

^ + 3Hn = {av) (ng - n{tf) (2.153) 
with the thermally averaged cross-section identified as 

{(jv) =\ I (inidn2(in';^(in2 x {2'kY5 (pi +P2 - p'l -P2 

X 2e-^('''i)/^2e-^(P^)/^(l - g{p^, t))(l - g{p2, t)) \Af . (2.154) 

On an isentropic Universe the entropy density will decrease by a factor . This can 
be used to scale out the Universe expansion efi'ect in Boltzmann equation by defining the 
quantity 



A trivial computation shows that 



y = -. (2.155) 
s 
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Y = - + 3H- (2.156) 
s s 



and Boltzmann equation in terms of Y is 

sY = {av) {Y^ - Y^) (2.157) 

where now (av) is defined with Yq. Usually the evolution of Y is given as a function of the 
temperature instead of time. From (2.110), (2.112) and the reduced variable x = m/T, 
where m is the mass of the particle in study, we get finally 

rIY 

sHx^ = {av) (Fcf - Y^) . (2.158) 

The Boltzmann equation in this form will be used in Chapter 3. With a little rearrange- 
ment it can be used to set a rough estimate for a specie to freeze-out. If the interaction 
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rate is much smaller than the expansion rate determined by the Hubble parameter 

r = no{av) <^ H, (2.159) 

where {av) is again given by (2.154), then the relative density Y may stabilize in a value 
different from the one given by the equilibrium density. In a scattering this just means 
that the expansion of the Universe makes more and more difficult the bath particles to 
meet and interact. On a decay it simply says that the particle's half-life is bigger than 
the age of the Universe. 



Relativistic particles contribute the most to the entropy density, so that in a very good 
approximation 

s = ^9*sT' (2.160) 
with ^ ^ 

i— bosons i^fermions 

When all the Standard Model degrees of freedom are active and the particles are in 
thermal equilibrium, then g^s = 106.75. On the other hand the equilibrium relativistic 
particle density can be obtained by integration of (2.130) with (2.107) and E = |p|. The 
result is 

f bosons 
n = I , 2.162 

\ fermions ^ ^ 

with C, the Riemann's zeta function. Then the equilibrium relative density Y^q for rela- 
tivistic particles will be 



45C(3) geff ^ Q ( 
2vr4 g^s ' \ 9*s 



Y,^ = = 0.278 ( ^ ) (2.163) 



where gefj = g for bosons and g^ff = 3g/4 for fermions. This can be used to obtain a 
cosmological bound on neutrino masses. The present value for the entropy density is [40] 



So = 2.97 X lO-' cm-^ (2.164) 

Neutrinos decouple when they are still relativistic, roughly at temperatures of few MeV 
[40] when g^s = 10.75, so that their relic density is related with Yeq by 

no = SoY^g = 825 ( = 115.12 cm"^ (2.165) 



and then 



^ m^rio 115.12 , 
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which gives 

777 

n,h^ = eV (2.167) 
91.5 ^ ' 

for just one family of a 2-component neutrino. The value quoted at PDG [43] for the 

relative matter density is Qmh^ = 0.133 ±0.006. Assuming 3 neutrino families we get the 

cosmological upper bound 

3 

J] m^i < 12.2 eV. (2.168) 

1-1 



2.4.3 Dark Matter 



If the particle decouples from the thermal bath when it is non-relativistic the precise 
determination of Y^o is more difficult. Following [40] we start by noting that av can be 
expanded in powers of v 

av(xvP, (2.169) 

with p = for s-wave annihilation, p = 2 for p-wave annihilation, etc... Because cc T 
we have that 

{av) = (Tore-", (2.170) 
for n = p/2, and Boltzmann equation (2.153) expressed in terms of Y becomes 

^ = -A.-"- (y^ - y^) (2.171) 



with 



and 



; (av) s(x)rrr' \ / , i /9\ , ■. 

A = I ^li/^ i^,, = 0-264 {g.s/g'J') M„mao (2.172) 

1.679/ m^MpiJ 



Y 

^ eq 



(x)=^=i^n'''— -^/v^ 

^ ' s{x) 27r4 Vs/ g,s 



= 0.145— x^/^e-^ (2.173) 
9*s 

This differential equation can be solved approximately. Start by defining the quantity 

A = y - Feg (2.174) 
and consider the corresponding differential equation 

^ = _ Ax— (2y,, + A) . (2.175) 

ax ax 
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At early times, for 1 < x <^ Xf where Xf is relative to the decoupling point, Y follows 
eg so that A and ^ 



closely Yeq so that A and ^ can be safely ignored. In this way 



rIV 

A = -A-ix"+2^ (2n,)-^ . (2.176) 

Differentiating (2.173) we get 

A ~ (2.177) 

At late times, for x ^ xj, y is very different from Y^q and we can take A ~ y :§> Yeq. 
Then 

^ = -Xx-^-'A' (2.178) 
ax 

which by integration gives 

11 A 1 



A{xf) n + la;"+i' 
Assuming that A(x/) — Aoo ~ ^{xj) we get 



(2.179) 
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A^ ~ -^^T'- (2-180) 

Now we must determine the decoupling point xj that is when Y starts to deviate from Yeq. 
We choose the criteria A{xf) = cYeq{xf) where c is a constant of order 1. Then solving 
the general equation for A (2.175) with the freeze-out criteria we get the approximate 
transcendental equation 

Xc{2 + c)ae-''f = in + ^)\nxf (2.181) 
with a = 0.145 {g/g*s) that can be solved by iteration. The first two iterates give 

Xf = In ((2 + c) Xac) -{n+^) In (In ((2 + c) Xac)) . (2.182) 

The choice c(c + 2) = n + 1 gives the best fit to the numerical integration [40] . With this 
choice we have 

X/ = In (0.038 (n + 1) {9/9^^) Mpima,,) , (2. 183) 

3.79(n+ l)x"+^ 
= \ ^ . (2.184) 

{9*3/9* )MpimoQ 
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Now it is easy to get the number and relative densities for a cold relic: 



(n + l)x]+' 

1 /2 

{g*s/g* )Mpimao 



% = so^'oo = 1.13 X 10^ ^ cm-3 (2.185) 



2 m^n 1.07 X 10^ GeV-^(n + l)xf 



n^h' = ^ = J- (2.186) 

{g*s/g* )Mpiao 

To get an order of magnitude for the mass of a cold relic let's assume that the interaction 
strength is of the order of the electroweak coupling and on dimensional grounds we may 
suppose that 

a, ~ — . (2.187) 



For s-wave and p-wave annihilations and with < 1 we obtain the order of magnitude 

~ 1 TeV. (2.188) 

That is why in the search for Dark Matter candidates, WIMPs (Weakly Interating Massive 
Particles) are well motivated. In this respect mSugra provides an attractive framework 
for these searches, not only because in most regions of parameter space the LSP (Lightest 
Supersymmetric Particle) has a mass on the hundreds of GeVs range but also because 
there are only 5 independent parameters 

tan/3, 1^1, mo, mi/2, Aq (2.189) 

from which the mass spectrum depends and this makes the phenomenological study sim- 
pler. 

Usually in most regions of parameter space the LSP is the neutralino and we will 
focus on this possibility. The dominant channels for the neutralino annihilation or co- 
annihilation depend on its nature. As explained in section 1.1.2 the neutralino mass 
matrix depends on the gaugino masses Mi and M2, that are obtained at the electroweak 
scale by the running of the RGEs with the boundary value mi/2 at the GUT scale and 
also on the /i parameter which is fixed by the electroweak breaking condition (1.40) and 
that depends on rrii and m2. These in turn are the solutions at the electroweak scale 
of the corresponding RGEs, so they are in fact functions of mo, and also of mi/2, since 
these equations are coupled with the gauge part. This means that what determines 
the predominant neutralino component are the parameters mo, mi/2 and tan/3. On large 
regions of mSugra parameter space the fi parameter is high, which means that the lightest 
neutralino will have a big gaugino component. But for tan /3 low (roughly of order 10 or 
below) and low to moderate values of mi/2 the RGEs are relatively insensitive to the 
value of mo, in the sense that for a wide range of initial values for mo, these equations 
give virtually the same negative value for m2, small in absolute value. This is known 
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as the focus point. Since in this region the value of ni^ is small in absolute value, the fi 
parameter obtained by the electroweak breaking condition (1.40) can also be small and the 
lightest neutralino for this region has a big higgsino component. In this case the relevant 
couplings are between neutralinos and the Z boson and neutralino, charginos and the W 
boson. The annihilation proceeds through s and t channels with the exchange of another 
neutralino or a chargino. Since the mass differences between Xi and X2 ^^'^ ^^^o between 
Xi and Xi are small, co-annihilation must also be taken into account. These processes 
can be efficient enough to reduce the neutralino density down to values compatible with 



In the case where Xi is mainly a bino, then the only gauge coupling is by U{1)y and 
the annihilation cross-section will be suppressed by a factor tan^ 9\y- However, for high 
values of mi/2 roughly of order 1 TeV and low values of mo, the sleptons are relatively 
light and co-annihilation between the lightest r and Xi is very active and can also reduce 
the neutralino density to acceptable values. 

For high values of tan /3 there is another process that should be taken into account. 
For these values of tan j3 the CP-odd Higgs boson mass is lowered to a few hundreds of 
GeVs. The amplitude A for the process in figure 2.3 goes as 



where F^ is the decay width of the A scalar, and it is resonant for the condition 



0.105 ±0.008 [43]. 



A 



1 



(2.190) 



4 - (mA/m^)2 -hiFAmA/m2 ' 



rriA 



(2.191) 




which can be fulfilled in some regions of the (mo,mi/2) plane and makes the annihilation 
of neutralinos very effective. 




Figure 2.3: Higgs Funnel 



Some explicit examples of these processes will be given in Chapter 6. 
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2.4.4 The Majoron as a Dark Matter candidate 

In [45] the interplay between spontaneous breaking of lepton number and neutrinos masses 
was fully addressed. We review here the discussion for completeness. The model consists 
of an SU{2) x U{1) gauge theory under which there are n species of neutrinos that are 
doublets for this gauge group and another n species which are singlets. These neutrinos 
are massive. We assume the most general situation where the masses are given by a 
mixed Type-I/Type-II seesaw mechanism and where the heavy neutrino mass matrix is 
obtained by spontaneous breaking of lepton number. In this way the model contains 2 
Higgs fields, a SU{2) triplet T with vev u and a doublet h with vev v. Additionally it 
contains a scalar field (f) with vev a that it is a singlet under this gauge group and that 
carries lepton number / = — 2. The vev of this field is responsible for the heavy neutrino 
mass matrix. The full neutrino mass matrix is 

M= J] (2.192) 

\mjy Mr J 

and the corresponding mass term in the lagrangian is 

Anass = -^P^icr2Mp+ h.c. (2.193) 

Here p refers to a 2n column vector of Weyl spinors, where the first n are SU{2) doublets 
and the last n are singlets. The term Mjj is typically generated by a Type-II seesaw 
mechanism, whereas rrio refers to Dirac terms. The physical fields u are obtained by a 
diagonalization of the mass matrix through a unitary matrix such that 




m 

rriR. 



real positive, diagonal, (2.194) 

(2.195) 

The full diagonalization procedure has been described in the op. cit. In the process the 
expansion parameter 

is of relevance. For Type-I seesaw Mr can be at most at the energy scale of 10^^ GeV 
and because niD is of the order of the electroweak scale, e can be as low as 10"^"^. In 
diagonalizing (2.192) we make the ansatz 
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where S is of order e and Vi and V2 are unitary and of order 1. Replacing (2.197) and 
(2.192) in (2.194) gives 

S = -im*o (M*)"^ (2.198) 



and 



Vi (^—moMff ^m^ + Mil) Vi = m = real positive, diagonal, (2.199) 

V2 (^Ir + -{M^)~^m\jmD + ^2 = "mn = real positive, diagonal 

(2.200) 



The transformation matrix U is then 




^ /(I - '^m},{M*^)-^M-'ml) V, m^,(M^)-V2 \ 

1^ -M-^mlV, {l-lM-V^mUM*i,)-')V2) 

+ 0{e^). (2.201) 

We now turn to the discussion of the Yukawa couplings between the neutrinos and the 
scalar fields. As the neutrino masses are obtained by a Type-I/Type-II seesaw mechanism 
and by the spontaneous breaking of lepton number, the Yukawa couplings are just the 
full neutrino mass matrix divided by the corresponding vevs: 



p+ h.c. (2.202) 



Due to the breaking of a global lepton number symmetry the physical model will contain 
a Goldstone boson, the Majoron. In order to identity the coupling of the Majoron with 
the light neutrinos we must first address the question of the minimization of the Higgs 
potential and identify the massless components of the Higgs fields. To do that we note 
that the scalar potential is invariant under the hyper-charge and lepton number gauge 
transformations. For example, in the first case we have the infinitesimal transformations 

6{Tr) = -2aYm) 
5{T,) = 2aY{Tr) (2.203) 




where r, i refer to the real and imaginary components, ay is the abelian gauge parameter 
and noting that Y(T) = 2. There is a similar expression for h with Y{h) = 1. Recall that 



78 



Baryon Asymmetry, Leptogenesis and Dark Matter 



is a singlet under SU{2). In this way, for hypercharge, we get 

dV dV ^ dV ^ dV , dV , , 

27^^. - 2T^r. + l^h, -^K = (2.204) 



day dTr dTi dhr dh 
and for lepton number we have 

dV dV , dV , dV . dV ^r, , , 

assuming that the triplet has I = —2. Differentiating this second relation with respect to 
and taking the vacuum expectation value we get, with the hypothesis of non spontaneous 
breaking of CP, 

f>-'«(5>-'^"'(^>- 

and using the minimization conditions we obtain 



Proceeding in exactly the same way we can derive four more relations 

+ ul—^) =0 



+ = (2-208) 



I d^v \ / d^v 

d(t)^dhy^''\dT^dh, 

These five equations are enough to express the 3 by 3 symmetric mass matrix of the 
imaginary components of scalar fields in terms of just one quantity, which can be chosen 
to be 

/ d^V \ 

"'{WW/^ (2^209) 

This means that in the basis T-\ h^, <jp^ this matrix is 



/ 1 JL 

' 2a 2m > 

V 

2a iau 
„2 



a. (2.210) 



, V_ _J£_ f / 

\ 2u iau 4u^ 1 

The explicit computation of the determinant shows that there is a null eigenvalue. In fact 
this eigenvalue has double multiplicity which means that there are two Goldstone bosons: 
one should be identified with the Majoron and the other is absorbed by the Z boson. 



2.4 Cosmology and Dark Matter 



79 



Since the Z does not couple to the scalar (f) the second Goldstone boson is proportional 
to 

vh^i+2uTf (2.211) 

while the Majoron is orthogonal to it: 

J = N {-2vu^h°^ + a{v^ + Au^)(j)i + uv^T^) , (2.212) 

with = (4f + (j2(w2 + 4n2)2 + u^v^)~^^^. We see that each of the fields h°, and 
contains a piece of the Majoron: 

(f,* = -iNaiv"^ + 4u^)J + ... (2.213) 
h° = -2iNvu^J + ... (2.214) 
= iNuv^J + ... (2.215) 

In this way, from (2.202), we see that neutrinos will couple to the Majoron. Specifically 
for the light neutrinos we have 

1 /T° 6* \ 

Ant = -;Ti^^i^2 —UlMjjUa + -(/7jm^f/, + V^mlUa) + —UjMnU, v + h.c. 
I \u V a J 

(2.216) 

Using the components (2.213), the expressions (2.201) and the definition of e we find that 
Aj = -]^NJv^v^ia2V^ {Mji - rriDM-^m^) V^u + h.c. + 0{e^) (2.217) 
Note that in obtaining this result it is important to observe that 



Mil ~ O(m^). (2.218) 

Using (2.199) we see that the coupling is proportional to the neutrino masses. 

As the Higgs doublet has some Majoron component, also charged leptons will couple 
to the Majoron. It can be seen that this coupling is 

^ j^^^^. (2.219) 

with m.0 the fermion mass and — (+) corresponds to a positively (negatively) charged 
fermion. There is also a one-loop induced coupling between the Majoron and photons 
that will be explicitly computed on section 4.5.2. This coupling is the most active in the 
Majoron decay since the Majoron may acquire a mass on the keV range and also because 
the coupling with neutrinos is proportional to neutrino masses. The present density nj(to) 
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of Majorons is 

nj{to) = nj{tn)e-'°/\ (2.220) 

where r is the mean half-life ant n{t]j) is the density at the decoupling time to, when 

Majorons start to depart from thermal equilibrium with photons. Assuming constant 

total entropy 5" from tjj to to, the entropy density varies as s oc This means that 
= ni?^ (X n/s and so 

rijM = (2.221) 



The photon density is related with the entropy density as 



(2.222) 



which means that (2.221) can be rewritten as 



njjto) _ g^sjto) nj{tD) 



(2.223) 



Assuming equal temperatures for the relativistic effective degrees of freedom and for 
to > 170 GeV, we have g^to) = 106.75 and g^to) = 3.91. Also n^(to) = 422 cm'^ 
When Majorons are in thermal equilibrium with photons, the ration nj/n^ is equal to 
1/2. Putting all this together gives 

Qjh^ = e-*°/- (2.224) 

1.36 keV ^ ' 

for the relic density. We have considered the simplest possibility, that Majorons where 
in thermal equilibrium when produced. If more complex scenarios are allowed, when can 
encode our ignorance on the production process by a parameter /3 and write 



Clearly the Majoron must be long-lived (r > tg) if it is to explain the Dark Matter density. 



Chapter 3 

Leptogenesis in an Model 



3.1 The Model 



In [46] it was proposed a model for neutrino masses with mixed Type-I and Type-II 
seesaw, based on the A4 symmetry group. The A4 or ahernating group is the rotational 
symmetry group of the tetrahedron or equivalently the group of even permutations of four 
objects. Even permutations are obtained from pairs of transpositions, so the generators 
of A4 can be taken as = (21) (43) and T = (312) which form a presentation of the group. 
Note that (231) is conjugate of (312) by the single transposition (32) which means that 
these two elements belong to different classes in A4, which are the inverse of each other. 
In summary, A4 has four classes 

Ci : / = (1234) (3.1) 
C2 : T = (2314), ST = (4132), = (3241), STS = (1423) 
C3 : = (3124), ST^ = (4213), T^^ = (2431), T^T = (1342) 
C4 : 5 = (4321), T^^T = (3412), T^T^ = (2143) 



Then it has four irreducible representations and a well-known relation for discrete groups 

N = ^dl (3.2) 

i 

where is the number of group elements and di is the dimension of the i-th irreducible 
representation, shows that it has three one-dimensional and 1 three-dimensional represen- 
tations. The characters are shown in table 3.1. The procedure to obtain the characters 
can be seen in [47]. Note that as expected C2 and C3, being inverse classes of each other, 
have characters that are complex conjugate. Note also the important relation 



1 + w + = 0. 



(3.3) 
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Table 3.1: Characters of A4 



Class 


x' 


x'' 




X^ 


Ci 


1 


1 


1 


3 


C2 


1 









Cs 


1 




iO 





C4 


1 


1 


1 


-1 





Li 


L2 


L3 




I'm 




A 


SU{2) 


2 


2 


2 


1 


1 


2 


3 


U{1) 


-1 


-1 


-1 


-2 





1 


2 


A, 


1 


V 


1" 


3 


3 


3 


1' or 1" 



Table 3.2: Lepton multiplet structure of the model 



The model of [46] is an extension of the Standard Model with right-handed neutrinos 
in the singlet representation of SU (2) and with an extra Higgs boson in the triplet rep- 
resentation. In table 3.2 the quantum numbers for the leptons and the Higgs bosons in 
the model are indicated. The Higgs potential for the Higgs doublets has been calculated 
in [48] and we present it here for completeness. Since the Higgs doublets belong to the 3 
dimensional representation of A4, we can produce invariants from 3 (8) 3: 

3 (8 3 = 1 © 1' © 1" © 3, © 3a (3.4) 

This means that in condensed form the Higgs potential can be written as 

y ($) = • ^ + Ai • J + A2 • ^, • ^>) ^„ (3.5) 

By a redefinition of the coupling constants we obtain the expression in 



V = m 



+ A2(<I>i^>l + UJ^^l^2 + t^^3^'3)('I'l^l + UJ^2^2 + ^^^^l^s) 
+ A3[($|<^3)(^J^>2) + (^'^$l)(<^t$3) + ($I$2)(^|^l)] 

(3.6) 



+ {^Hi'^l^sf + i'^l'^if + {^{'^2)'] + /i.e.} 
It has been shown in [48] that the vacuum alignment 



(<!>?> = ($°> = ($°> = j= (3.7) 
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is a solution for the minimization of (3.6). The Lagrangian pertinent to neutrino masses 
is 



+A (^L^Cia2a ■ ALs^ + A (^LlCia2a ■ ALi j 



+A LlCia2a • AL3 + h.c. 



(3.8) 



where the ^4 representations for each term are shown explicitly. Assuming that the 
triplet acquires a vacuum expectation value we obtain a mixed Type-I/Type-II seesaw 
mechanism for neutrino masses. The charged lepton and Dirac neutrino masses are 

Ml = v diag(/ii, /i2, hs)U 
mo = v diag(/ii£,, h2D, hsD)U , 



with 



/ 1 1 1 \ 



= e 3 



As a result, the Type I neutrino mass matrix after electroweak symmetry breaking is 



(3.9) 



M 



I II'ID 



\ 

h2DhiD 

\ h2Dh-iD j 



(3.10) 



On the other hand, a small vev u for the triplet induces a Type II contribution to the 
light neutrino mass matrix 



M 



II 



/ An \ 

An 
\ X'u J 



(3.11) 



where A, A' are two Yukawa couplings and for definiteness we are assuming that the triplet 
belongs to the l" representation of A4. The total neutrino mass matrix is given by the 
sum of equation (3.10) and (3.11) and has the form 



/ a X 0\ 

X b 
\0 b y J 



(3.12) 



where a, b and x, y refer to the type-I and type-II contributions, respectively. This provides 
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a simple derivation of the two-zero texture classified as Bi in Ref . [49] . As it is referred 
in [46], the main feature of two-zero texture models, such as the ones derived here, is 
their power in predicting the as yet undetermined neutrino parameters. Current neutrino 
oscillation experiments determine two mass splittings Aml^.^ and Am^Qj and the corre- 
sponding mixing angles 612 and 623, with some sensitivity on ^13 which is bounded [28]. 
The Dirac CP phase will be probed in future oscillation experiments. Similarly, the abso- 
lute neutrino mass scale will be probed by future cosmological observations [50], tritium 
beta decays [51] and neutrinoless double beta decay experiments [52] with improved sen- 
sitivity. The latter will also shed light on the two Major ana CP phases which are hard to 
test otherwise, as they do not affect lepton number conserving processes. The general 3x3 
light neutrino mass matrix A^j, in the flavour basis contains a priori nine independent real 
parameters, once the three unphysical phases associated with the charged lepton fields 
are removed. In contrast, in the proposed model all the above nine parameters are given 
in terms of only five unknowns. Hence the number of physical parameters characterizing 
the charged current weak interaction is reduced with respect to what is expected in the 
general case [53]. 

3.2 Triplet Decays and Leptogenesis 

Now we turn to the study of leptogenesis in this model. As our model has not only 
right-handed neutrinos but also Higgs boson triplets, the analysis of leptogenesis is more 
complicated than in the usual type-I seesaw case. This case has been studied, in a generic 
way, in Ref. [54] and we will use their results. To make contact with their notation we need 
to fully specify our model, giving the couplings between the SU (2) triplet and doublets. 
These are 



where i = 1,2,3 is an A4 indices and fi = /(I, u;^). These two Lagrangians, Eq. (3.8) 
and Eq. (3.13), are gauge invariant for the quantum numbers given in Table 3.2. With 
the identifications Ai + 1/^2 = A°, Ai - 1/^2 = A++ and A3 = A+/v^, and the A in the 
1" representation, the Lagrangian Cl is, explicitly, 




(3.13) 



= hiD VilVircIP* - hiD llLl^lR<Pi + ho T^lLt^2R(pf " ho Iil1^2r4'2 

+ hiD T^1L1^3R<PT - hiD hL'^3R(p3 + ^2£i T^2LI^Ir4>^i " h2D kLl^lR<Pl 

+ U?h2D 'i^2L'^2R(pf - W^/i2£l 'l2Ll^2R<p2 + ^^^D 1^2L'^3R(I>T ~ ^2Li^3fl^3 

+ h3D T^3L1^1R<PT - ho hLl^lRfPl + ^ho T^3Ll^2R<p2* " hL'^2R(p2 

+ LU^ h3D I^SL^^SiJ^r - ^'^ ho 13L^3R(I)3+^^^^2lC^1L " ^A+I^J^C/iL 
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V2 



On the other hand, Ca/s can be written as 



(3.14) 



-Ca/^ = fA-4>tct>t + ujfA--(f>tct>t + uj'fA-(t>t(t>t 



/A°>V? - ^/A°>°</)° - u;'fA'*cl>l4 + h.c. 



(3.15) 



The expressions for the asymmetries are normahy written in the basis where the right- 
handed neutrinos are diagonal. In our case tliis already happens because the right handed 
neutrinos are fully degenerate as a consequence of the A4 symmetry. Translating our 
notation into the one of [54] we get 



/ A \ 
1a = A 

V A' y 



(3.16) 
(3.17) 
(3.18) 



The diagrams contributing to the CP-asymmetry are shown in Fig. 3.1. The asym- 
metries for these diagrams have been computed generically in [54]. The first two do not 
depend on the triplet and give 



Stt ^ 



Im 


(i'ivi"/If)fcj4j 




{YN)ki 


2 



(3.19) 



where f{M) is some kinematical factor. The third diagram depends on the triplet and is 
given by 



1 T.ii^'m[{yN)uyN)ki(yi)ii{f^k 



27r 



N ki\ 



X 1 



Ml 
M2 



log(l + MVMA2) , (3.20) 



while for the fourth diagram we have 



E,:,l(^A),,PAfi + EJ^,:|^ 



xlog(l + Mi/M2), 



(3.21) 
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Figure 3.1: One Loop diagrams contributing to asymmetry from tlie i/ji decay. 



With the couphngs in Eq. (3.16) we can see immediately that the diagrams with 
only right-handed neutrinos do not contribute to the asymmetry due to the A4 group 
properties. The same is not true for the other asymmetries. As we see that without the 
triplet there is no asymmetry, we will concentrate in the case where all the asymmetry 
comes from the fourth diagram in Fig. 3.1, given in Eq. (3.21). As it was discussed in 
Ref. [54] this corresponds to the case where A^a ^ 

We have performed a scan of the parameter space subject to the condition that all 
the points satisfy the neutrino data [14]. The results for the asymmetry ea are shown in 
Fig. 3.2 as a function of the triplet parameters, the mass A/a and an effective coupling 
defined as 

\l = V2^Ti (y^yI^ (3.22) 

We see that the asymmetry can easily be high. However, to be able to discuss if the 
washout effects are significant, we have to solve the Boltzmann equations for our model. 

For their study we will also need the decay rates of the triplet into leptons and Higgs 
bosons. The tree-level decay rates are. 



r(A ^ LL) 
r(A ^ HH) 



(3.23) 
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where we have defined the dimensionless effective coupHngs 



X,^V2jTr{Y^Yl), ^ V2 ^^SIMUK (3.24) 
In terms of these coupUngs the branching ratios into leptons and Higgs bosons, B^, Bh, 

We solve tlie Boltzmann equation following the setup of Ref.[55]. Using their notation 
we recall that the Boltzmann equations describe the evolution as function of z = M/^/T 
of the total triplet density, = (n^ + njT)/s, and of the asymmetries Ap = {up — np)/s 
for the species p = L, H, where Up is the number density of the type p particles, and s 
is the total entropy density ^. With this notation the Boltzmann equations are 

sHz'-^ = X- 27.5.(^-0) (3.26) 



^In order to follow as close as possible the notation of Ref.[55], in this section we denote the triplet 
fields by T instead of A. 
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where H is the Hubble constant at temperature T,Yp = Up/s, a suffix '^'^ denotes equilib- 
rium, 7p is the spacetime density of type p computed in thermal equilibrium. The other 
quantities are 

X = ^n^L - 1) - 2 + ^) (7r ' + 7t*) (3.27) 

where and 7Tt, are the spacetime densities for the s-channel (properly subtracted to 
avoid double counting of the processes already counted on the decay, see below) and t- 
channel AL = 2 scatterings, respectively. We recall that the reaction densities are defined 
through 

7 = ^ dsVTsK, (4) -is) (3.28) 



647r4 7, . " Mr 
where Ki is a Bessel function and a the reduced cross-section defined as 

a{s) = Y,I rfi^ (3.29) 

^-^ J OTTS 

where the sum runs over initial and final spins and gauge indices of the amplitude A. 

The expressions for 7^ and 7^) are the same as in Ref.[55] (with the obvious exception 
that we have three Higgs bosons) and we do not reproduce them here. There are compu- 
tational differences for AL = 2 scatterings but the order of magnitude for these processes 
remains similar. 

The baryon asymmetry is then obtained as 

riB 



Vb= — 

where 



today S 



s 

today ^7 



S 



TT^ 43 
today " 4503) n 



Ai (3.30) 

today 



(3.31) 



and Osph = 36/105 is the sphaleron conversion factor for three doublets of Higgs bosons. 
We therefore have to solve the Boltzmann equations, Eqs. (3.26) to find the lepton asym- 
metry Ai and then we obtain rjB using Eq. (3.30). Due the complexity of Eqs. (3.26), 
they have to be solved numerically. Although our model differs in the details from that 
of Ref. [55], it shares many of the same features that are due to the presence of the triplet 
in the Boltzmann equations. Therefore we recover their result that the efficiency can be 
high in this type of models. 

We have performed a scan of the parameter space with the free parameters in the 
following ranges 

Ma e [10^ 10^^] GeV, |/| e [10^ 10^^] GeV, M G [10, 100] Ma (3.32) 
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The other parameters, including the coupling matrix Ya, were obtained from the require- 
ment that the neutrino masses and mixings were within the present allowed 3(t range [14] . 
The results are summarized in Fig. 3.3, where we plot contours of equal baryon asymme- 




10"^ 10""^ 10"^ 10"^ 10"'' 10° 

Figure 3.3: Contour plots for Xogi^rjE as a function of -A/a and Al, after the effect of washout 
has been taken in account by solving the Boltzmann equations. For comparison we note that 
(at 95% CL) we have: -9.35 < \ogn^r]B < -9.18. 

try, 77b, as a function of the triplet mass Ma and the effective coupling, A^, defined in 
Eq. (3.24). In Fig. 3.4 the contours are shown as functions of the effective couplings 
and Xh ■ We conclude that in a large region of parameter space the present value obtained 
from Big Bang Nucleosynthesis (@ 95% CL) [56], 

4.7 X 10-1° <r]B< 6.5 x IQ-^" (3.33) 

can be easily achieved, accounting at the same time for the present experimental limits 
on neutrino masses and mixings[14]. 

3.3 Higgs Potential Minimization 

Until now we have just assumed that the triplet acquires a small vev, breaking in this 
way A^. However, a careful analysis of the Higgs potential in the presence of the triplet 
shows that this is a non trivial matter^. Following a suggestion by A. S. Joshipura^, let 

^We thank F. Joaquim for pointing this out to us. 
^A. S. Joshipura, private communication. See also [22]. 
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Log-|o|riBl as a function of and 




10"^ 10"^ 10"'' 10' 



Figure 3.4: Contour plots for log^QriB as a function of and Xjj, after the effect of wasiiout 
iias been taken in account by solving the Boltzmann equations. For comparison we note that 
(at 95% CL) we have: -9.35 < logior/B < -9.18. 



us write the full Higgs potential 

V = V{^) + T/i($, T) + F2(^>, T) + V{T) (3.34) 
with V{^) given by (3.6) and Vi(<I>,T) the trilinear coupling 

yi($, T) = Art ^(c,$f $,) + h.c. (3.35) 

i 

in a simplified notation omitting SU{2) indices, with Ci = {l,w,w'^). The terms V2($,T) 
and V{T) refer to quartic couplings and the triplet mass term and its explicit form will 
not play an essential role in this discussion. The minimum condition with respect to T is 



5V 
6T 



The first term is essential to give a vev to the triplet, but with the Higgs doublets vacuum 
alignment (3.7) it vanishes because of the A4 relation 



1 + w + = 0. 



(3.37) 
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There is also a second problem. The minimization with respect to the gives 

— = AXu^a + ' + 3.38 

V3 

Noting that the last two terms are independent of i at the minimum and summing the 
three equations for ii = 1, 2, 3 we get at the minimum, using (3.37), 

5V 6V ^ _ Q 

(5<l>2 

O 3a = (3.39) 

with 



Replacing in (3.38) we conclude that 

Xuvc, = (3.41) 

which means that the two vevs cannot be simultaneously non zero. 

To solve the first problem A. S. Joshipura proposed to introduce an explicit breaking 
term 

A'Tt J]($f$,)i+ h.c. (3.42) 

i 

with a' small, where we choose the product of the two Higgs SU (2) doublets A/^ triplets 
to be in the trivial representation of A^. The triplet vev would break ^4 anyway, although 
spontaneously, so this just amounts to break ^4 already in the tree level Lagrangian. This 
avoids the relation (3.37) and it allows the triplet to acquire a vev. For the second prob- 
lem we propose to introduce a second triplet in the l' representation and with opposite 
hypercharge of T. Let us call this new triplet T2 and the initial one Ti. Then the A4 
invariant cubic couplings between the Higgs triplets and doublets are 

AiT/ 5^(c,:$f $,) + \2T2 J^(Q$f $,) + h.c. (3.43) 

i i 

which permits to replace (3.41) by 

XiUiVCi + \2U2VCi = 0. (3.44) 

This no longer requires the vanishing of one of the vevs. 

This second solution turns the model less predictive because there are more param- 
eters. But in a supersymmetric framework the second triplet would not couple to the 
leptons because of the superpotential holomorphy, which means that the Type-II seesaw 
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would still be determined only by one triplet-leptons coupling. On the other hand it is 
necessary to introduce a second family of SU{2) Higgs doublets for the cancelation of 
anomalies, which turns the model more elaborated. Alternatively, since the two triplets 
have different quantum numbers, it is licit to make the assumption that one is much 
heavier than the other and Type-II seesaw would still be determined mainly by the light- 
est triplet, assuming small to moderate Yukawa couplings. Of course, leptogenesis in 
this framework would be much more involved as there are now more processes engaging 
the two Higgs triplets. This will be the subject of future work. Another possibility of 
reconciling the A4 symmetry with the Bi zero texture is explored in the next chapter [22] . 



Chapter 4 



744-based neutrino masses with Majoron decaying 
dark matter 

4.1 Introduction 

As we saw in Chapter 1, the discovery of neutrino oscillations [10, 57-60], now confirmed 
at reactors and accelerators [61-63], has brought neutrino physics to the center of particle 
physics research. Global analysis of current oscillation data indicate that the pattern of 
lepton mixing differs sharply from that characterizing quarks [14]. Understanding the 
origin of neutrino mass and the pattern of neutrino mixing angles from basic princi- 
ples constitutes a major challenge [64, 65]. A paradigm framework to generate neutrino 
masses is provided by the seesaw mechanism, for which several realizations have been pro- 
posed [66]. The observed pattern of neutrino mixing may arise from suitable non-abelian 
flavour symmetries, as those based on the A4 group [48, 67-69]. 

Elucidating the nature of dark matter constitutes another intriguing problem of mod- 
ern physics which has so far defied all efforts. It is therefore crucial to build a fundamental 
particle physics theory of dark matter and, since the Standard Model of elementary parti- 
cles (SM) fails to provide a dark matter candidate, such theory necessarily requires physics 
beyond the SM. 

In the work [22], that it is described in this Chapter, we suggest a version of the 
seesaw mechanism containing both type-I [53, 70-76] and type-II contributions [53, 75, 77- 
80] in which we implement an A4 flavor symmetry with spontaneous violation of lepton 
number [75, 76]. We study the resulting pattern of vacuum expectation values (vevs) and 
show that the model reproduces the phenomenologically consistent and predictive two- 
zero texture proposed in Ref. [46], avoiding the problem of the vanishing of the triplet 
vev discussed in section 3.3. 

In the presence of explicit global symmetry breaking effects, as might follow from grav- 
itational interactions, the resulting pseudo-Goldstone boson - Majoron - may constitute 
a viable candidate for decaying dark matter if it acquires mass in the keV-MeV range. 
Indeed, this is not in conflict with the lifetime constraints which follow from current cos- 
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mic microwave background (CMB) observations provided by the Wilkinson Microwave 
Anisotropy Probe (WMAP) [81]. We also show how the corresponding mono-energetic 
emission line arising from the sub-leading one-loop induced electromagnetic decay of the 
Majoron may be observed in future X-ray missions [15]. 

4.2 The Model 

Our model is described by the multiplet content specified in Table 4.1 where the transfor- 
mation properties under the SM and groups are shown (as well as the corresponding 
lepton number L) . The Li and Ijn fields are the usual SM lepton doublets and singlets and 
Uj^ the right-handed neutrinos. The scalar sector contains an SU(2) triplet A, three Higgs 
doublets (which transform as a triplet of A4) and a scalar singlet a. Three additional 
fermion singlets Si are also included. 



Table 4.1: Lepton multiplet structure {Q = T3 + Y/2) 





Li 


L2 


^3 








A 


a 




SU{2) 


2 


2 


2 


1 


1 


2 


3 


1 


1 


U{1)y 


-1 


-1 


-1 


-2 





-1 


2 








A, 


V 


1 


1" 


3 


3 


3 


1" 


1" 


3 


L 


1 


1 


1 


1 


1 





-2 


-2 


1 



Taking into account the information displayed in Table 4.1, and imposing lepton num- 
ber conservation, the Lagrangian responsible for neutrino masses reads 

+ ALfCALa + ALf CALi + A'L^CALg 

+ Mr (SIur) ^ + h {SICSl)[ a + h.c. , (4.1) 
where h and A are adimensional couplings, Mr is a mass scale and 

A=( -A+/v^ 
\ -A+/V2 A++ 

Note that the term (//^Cz/^J'^a is allowed by the imposed symmetry. This term however 
does not contribute to the light neutrino masses to the leading order in the seesaw ex- 
pansion and we omit it. Alternatively, such term may be forbidden by holomorphy in a 
supersymmetric framework with the following superpotential terms 

w = ■ ■ • + Xcabh'^L^o'Ht + Mru's + hssa 

where the hats denote superfields and the last term replaces the corresponding bilinear 
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employed in Ref. [82, 83]. Assuming that the Higgs bosons $i, A" and a acquire the 
following vevs (see section 4.3 below) 



(4.3) 



we obtain an extended seesaw neutrino mass matrix M. [82-84] in the (z/^, i^^, S) basis 

/ 1 u;2 ^ \ 



/ 



M 



\ 



M 
\ n ) 



niD = V diag(/ii, h) U, U = —j= 

v3 



1 1 1 



V 1 



J 



(4.4) 

with UJ = e^'^*/'^, M = M^j diag(l, 1, 1) and /i = u^jh disig{l,w^, w). This leads to an 
effective light neutrino mass matrix given by 



( hi \ 

h2h 

\ h2hs / 



(4.5) 



On the other hand the vev of the triplet, ua, will induce an effective mass matrix for the 
light neutrinos from type-II seesaw mechanism 



/ A \ 



M]} = 2 



A 

\o A' y 

and the total effective light neutrino mass matrix will then be 



(4.6) 



(4.7) 



In Ref. [46] it was shown that the neutrino mass matrix given by Eq. (4.7) could explain 
the currently available neutrino data. In section 4.4 we will present an update of that 
analysis taking into account the latest neutrino oscillation data. 



4.3 A4 Invariant Higgs Potential 

As in Chapter 3, 3.3, we now address the question of the minimization of the neutral 
Higgs scalar potential in this new framework, which, as we saw, is a necessary condition 
to reproduce the structure of the neutrino mass matrix presented in the previous section. 
With the assignments of Table 4.1, the Higgs potential consistent with gauge and A4 
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invariance and lepton number conservation reads, 

l/ = y($) + y($,A,a), (4.8) 

where V{^) is given as"^: 

= ml ^ + Ai ($t$) ^ ($t$) ^ + A2 ($^«>) ^, ($^$) ^„ 

+ A3 (^^^)3. • (^^^)3. + A4 ($^$)3, • ($^$)3, + A5 ($t$)3^ . ($t.5)^^ , (4.9) 

and V{^, A, cr) contains pure A , a terms, together with others involving mixed invariant 
combinations of the scalar fields. Assuming the so-called seesaw hierarchy <C v <C 
u„ [75] ^, the relevant terms in A, cr) are ^ 

y(«>, A, a) = {Ml + p|(j|2) Tr(A^A) + A^|a|^ +[ml + ^ ($^$) J \af - (5$^A$a* + h.c), 

(4.10) 

Taking the vacuum alignment for the Higgs doublets <J>a given in equation (4.3) the min- 
imization of the Higgs potential with respect to A gives 

SV 

— = 0^(^Ml + pul)uA-5v\, = 0. (4.11) 

We stress that the A4 symmetry, together with the doublet vev alignment assumed in 
Eq. (4.3), requires that the product <l> (g) $ ~ 1 under A4. If ^> (g) <I> ~ 1', 1", then the 
second term in the above equation would reduce to 25{l+uj + uj^)ucr = implying ~ 0. 
Moreover, as a direct consequence of the requirement <!> (g ^> ~ 1 under A and a must 
have the same (singlet) transformation properties under that group. 

The above equation leads to the following solution for the triplet vev 

5v'^Ucr Sv^ . . 

where the last approximation holds for Ma <C u„. This result shows that the "vev- seesaw" 
relation u^Ug. ~ v"^ is fulfilled. The minimization with respect to the gives 

^ = ^ + 2^vul - AdvuAU^ = 0. (4.13) 



^The decomposition of the tensorial product of two triplets in A4 is shown in (3.4) 
^In contrast to the inverse seesaw models used in Refs. [83, 84] here we consider large values of u^r, 
u„ > 10^ GeV or so. 

^Notice that the scalar potential contains other invariant terms such as $t(j)Xr(A'''A), Tr(AtA)|(Tp, 
[Tr(AtA)]^, etc. Assuming the vev hierarchy ma << v << and that the adimensional coefficients of 
these terms are of the same order of the ones in V{^, A, cr), then V{^, A, cr) is enough for our purposes. 
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Finally, 

AT/ 

= ^ 2X^ul + (ml + + pul) - 25vW = 0. (4.14) 
which, in the limit U/^,v « u^, has the approximate solution 




(4.15) 



as it is typical from spontaneous symmetry breaking scenarios. In summary, we have 
shown that in our framework it is possible to achieve a consistent minimization of the 
scalar potential with non-zero vevs satisfying the "vev-seesaw" relation u^Ucr ~ v'^. 



4.4 Neutrino parameter analysis 

Given the two contributions to the light neutrino mass matrix discussed in Eqs. (4.5) and 
(4.6) one finds that the total neutrino mass matrix has the following structure: 
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a 


b 





\ 






b 





c 






I 





c 


d 


J 



This matrix with two-zero texture has been classified as Bl in [49]. One can show that con- 
sidering the [Li, L2, L^) transformation properties under A4 as being (1', 1", 1) or (1", 1', 1) 
an effective neutrino mass matrix with A4^{1,2) = A^y(3,3) = is obtained (type B2 
in [49]). Moreover, by choosing A, a ~ 1' and appropriate transformation properties of 
the Li doublets, we could obtain the textures Bl and B2 as well. Still, the configuration 
A, a ~ 1 would lead to textures which are incompatible with neutrino data since, in this 
case, both type I and type II contributions to the effective neutrino mass matrix would 
have the same form. Since the textures of the type Bl and B2 are very similar in what 
concerns to neutrino parameter predictions, we will restrict our analysis to Bl, shown in 
(4.16). 

As was already mentioned in section 1.2.3, the neutrino mass matrix is described by 
nine parameters: three masses, three mixing angles and three phases (one Dirac -|- two 
Majorana). From neutrino oscillation experiments we have good determinations for two 
of the mass parameters (mass squared differences) and for two of the mixing angles {612 
and ^23) as well as an upper-bound on the third mixing angle ^13. Using the 3a allowed 
ranges for these five parameters and the structure of the mass matrix in Eq. (4.16) we 
can determine the remaining four parameters. The phenomenological implications of this 
kind of mass matrix have been analyzed in Refs. [46] and [85]. Here we will update the 
results in light of the recently determined neutrino oscillation parameters [14]. 
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Figure 4.1: Correlation between the neutrinoless double beta decay amplitude parameter \mee\ 
and the atmospheric mixing parameter. Experimental sensitivities are also given for comparison. 



The main results are shown in Figs. 4.1 and 4.2. In figure 4.1 we plot the correlation 
of the mass parameter characterizing the neutrinoless double beta decay amplitude: 

|mee| = \cl3cl2m1 + c?3Si2m2e^*" + sj^mse^'^l , (4.17) 

with the atmospheric mixing angle 6*23. Here and Sij stand for cos 6ij and sin6'y 
respectively. At the zeroth order approximation rrii/ms = tan^ 6*23, and therefore ^23 < 45° 
for normal hierarchy (NH), while ^23 > 45° for inverted hierarchy (IH). The main result 
from this plot is a lower bound on the effective neutrino mass:|meel > 0.03 eV. For 
comparison the range of sensitivities of planned experiments as well as current bounds is 
also given. Note that the lower bound we obtain lies within reach of the future generation 
of neutrinoless double beta decay experiments. 




sin sin 9^^ 

Figure 4.2: CP violating phase 6 and CP-invariant J in terms of the reactor mixing parameter. 
The 3 (T-excluded range for sin"^ 9ij is given for comparison. 

The panels in Fig. 4.2 show the CP-violating phase 5 and the corresponding CP- 
violating invariant in neutrino oscillations, already discussed in section 1.2.1 - see formula 
(1.146): 

J = Sl2S23.Sl3Cl2C23C?3 siu 6 , (4.18) 
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versus sin^ ^13 . Note that these hold both for normal and inverted hierarchy spectra. 



Moreover, for large Ois values, where CP violation is likely to be probed in neutrino os- 
cillations, one can see that our model predicts maximal violation of CP. Quantitatively, 
from the right panel one sees that the 3a bound on Oi^: sin^ ^^13 < 0.053 implies an upper 
bound: \J\ < 0.06 on the CP-invariant. 

In addition, the two-zero texture structure of our neutrino mass matrix may have other 
implications, for example for the expected pattern of lepton flavor violating decays. In 
fact, thanks to the strong renormalization effects due to the presence of the triplet states, 
the latter are quite sizeable in sypersymmetric models [5, 86, 87]. 



4.5 Majoron Dark Matter 

In models where neutrinos acquire mass through spontaneous breaking of an ungauged 
lepton number [75, 76] one expects that, due to non-perturbative effects, the Nambu- 
Goldstone boson (Majoron) may pick up a mass that we assume to lie in the kilovolt 
range [88]. This implies that the Majorons will decay, mainly in neutrinos. As the coupling 
gji^iy is proportional to ^ [75], the corresponding mean lifetime can be extremely long, 
even longer than the age of the Universe. As a result the Majoron can, in principle, 
account for the observed cosmological dark matter (DM). 

This possibility was explored in Refs. [89, 90] in a general context. Here, we just 
summarize the results. It was found that the relic Majorons can account for the observed 
cosmological dark matter abundance provided 



where Pj^i^ is the decay width oi J ^ uv and mj is the Majoron mass. The parameter /3 
encodes our ignorance about the number density of Majorons, being normalized to /? = 1 
if the Majoron was in thermal equilibrium in the early Universe decoupling sufficiently 
early, when all other degrees of freedom of the standard model were excited [90]. In the 
following we will follow their choice and will take 




Pj^^ < 1.3 X 10"^^ s"^ , 0.12 keV < Pnij < 0.17 keV, 



(4.19) 



10"^ < /3 < 1 



(4.20) 



and calculate both the width into neutrinos as well as the subleading one-loop induced 
decay into photons. 
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4.5.1 Decay into neutrinos 

We now proceed with the computation of the Majoron decay width into neutrinos, which 
win be useful to obtain the allowed parameter space for which the Majoron can be a viable 
DM candidate. In order to calculate the decay amplitude we recall that the coupling gj^^,^j 
is defined through 

= --g.ji>,u,Ji'iVj + ^-c. (4.21) 
2 

For the evaluation of Qjv^Uji we follow the steps developed in Ref. [75]. First we notice that 
with scalar potential defined in section 4.3, the Majoron, in the basis [lm(0°), Im(A''), Im(cr°)] 
is given by 



T 



J = Nj 

and 



2^^a;^> ^"^71' ^"^71' ^'"'^^"^ + 



Nj = [4«2ni + v\l + ul{Aul + v'f] c -±- , (4.23) 

where the last equality follows from the assumed hierarchy <C w <C tio- implied by the 
vev-seesaw relation. Using this, one can obtain 

5..., = -^, (4.24) 



(4.22) 



leading to the decay width 



(4.25) 



327r 2ul 

It is worth mentioning that the sum X^,;(?Ti^)^ is in our framework constrained by the 
special form of the effective neutrino mass matrix shown in Eq. (4.16). In particular, 
there is a lower bound on the mass of the lightest neutrino: m > 0.03 eV, as we saw in 
section 4.4. 



4.5.2 Decay into photons 

The Majoron also couples with photons (at the quantum level) and therefore the radiative 
decay J — >■ 77 is expected to occur with a photon energy ~ 'n%j/2. Consequently, this 
decay exhibits a mono-energetic emission line which could be detected in a variety of 
X-ray observatories, see for example the discussion given in Refs. [15, 90]. 
The effective Majoron-photon interaction can be written as 

C = gj^^e^'^^F^^F^p , (4.26) 

resulting from the one- loop diagrams shown in Fig. 4.3 (top diagrams). The effective 
coupling gj^^ (bottom graph in Fig. 4.3) is 
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f{k-p) 




fik+p) 



•C 8z gj^^ e^yc^p grg2e"(9i)e (92) 



7 



Figure 4.3: Top: One loop diagrams for the decay J — > 77. Bottom: Effective J77 vertex. 



5^777 



(4.27) 



with X/ = —2rn^j:Co{0,Q,rn'^j/nij^rn'j,rn^j:) ~ 1 + mj/(12m^) where Co is the invariant 
Passarino-Veltman loop function [91]. The last approximation is valid for rrij <C ruf. 
Tg^, Qf and Nj denote the weak isospin, the electric charge and the colour factor of 
the corresponding charged fermion /, respectively. The coupling of the Majoron to the 
charged fermions gjfj is given by [90] 



9Jff 



2ul 



-mf{-2Ti) 



(4.28) 



We then get for the decay width, 



m , 



f 



77 



647r3 
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NfQjffQjXf 



9 3 
a rri j 

647r3 



/V^S^(-2T/) 



m , 



I2m) 



(4.29) 



where the cancelation of the anomalous contribution has been taken into account. The 
above calculation was performed using the Mathematica package FeynCalc [92] . Here we 
give a more pedestrian derivation of (4.29). Take the diagrams in figure 4.4. The integral 
associated with the left diagram is 



d^q Tr (^75(^ + ^ + mf)Yif^ + ^ + mf)Y{^ + rUf) 
(2^)4 ((p + qf _ mj) ((p2 + qf - mj) {q^ - m}) 



(4.30) 
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and the right diagram gives an integral T^j that is obtained from this one with (/i, u) — > 
Pi, (z^, fJ')- 

The trace is easily computed: 

Tr {j5{f + i + "^/)7''(|>2 + i + rnf)Y{i + m^)) = (4.31) 

m/Tr (75(^ + ^)7^7'^) + 
m/Tr(75(^ + ^)7''(^2 + ^)^'') = 

- Umfe^'^P" Hp + q)pq, - (p2 + q)pqa - {p + q)p{p2 + q)<y) = 
which gives for Tj" 

Tr = -4^m,e'^^^>l,p2. / ,% , . ^ \ ^y— ^. (4.32) 

J (^TTj'' ((p + g) _ ((p2 + g) -mj){q^-mj) 

We see that the integral is in fact convergent, although naively we would expect it to be lin- 
early divergent. Also, this expression is symmetric under the interchange {pi,P2), (m, z^) — > 
(P2,Pi), (z^, m) and the change of variables q — > —q — p in the integration, which gives 
■ So, we have 

T^'' + (4.33) 

=2Tr 

d^q 1 



(27r)4 ((p + qf - m)) {{p^ + qf - m}) (q^ - ru' 



Now we combine denominators with Feynman parameters to perform the integration: let 
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us call / the scalar integral 



We have 

"•l — X 

(2vr; JO ^0 



1=2 / / / dy (4.35) 



(x ((p + qf - mj) + y [{p2 + qf - mj) + {1 - x - y){q^ - mj))^ 
Using = m?j and P2 = the denominator can be cast in the form 

denom. = q^ — + x(l — x — y)m^j = (f' — C, (4.36) 
with C = rnj — x{l — x — y)m^j- Then the scalar integral is 



Jo 



(27r)4 {q2-Cf 

Now we change to Euclidean momentum by a Wick rotation q^ — > iqoE to get 



1 r'^-x r ,4 



/ =2i I dx I dy I ^ (4.37) 



JO 



10 Jo 

nl nl—x r roo 

-2i dx dy / dfl^ / du 
Jo Jo J Jo 

f'l f'l—x n 

-2i dx dy / dfl^ 



1 u 



Jo 



2(27r)4 (n + C)3 

1 



Using the result 



we get 



4(2^)4^" 

J dn^ = 2n^ (4.38) 
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Remembering that C = rnj — x{l — x — y)m^j we then have 



C-^ = nif il-x{l-x-y)^\ (4.40) 



m 



mjUl + x{l-x-y)-^ + ... I (4.41) 



my 



and puttmg this in (4.39) we obtain 



1 ri-^ / \ 
dx dy\l + x{l-x-y)^ + ...\ (4.42) 



m 



l + + (4.43) 



whicli gives for T^'^ 



T^"' = -8imfe^''P''pipP2aI (4.44) 

flipper 



4 



Noting that the Jff couphng is gjff, that the electromagnetic one is proportional to eQj 
and summing over all fermions we are able to write the M matrix. The Feynman rules 
applied to the diagrams give 

M = -Y.9JffNfQyT'^''e,{p,)e,{p2), (4.45) 
/ 

where e^(p) is the photon polarization vector and we can add the amplitudes for the two 
diagrams because the analytic expressions are the same. The coupling gjfj is 

9Jff = -^"^/(-2^3/) (4.46) 
which allows us to eliminate the anomalous like contribution, giving 

M = '^f"^M>AP2)P.,P2. E ^fQ]{-1T,,) [^^^ (4.47) 
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This gives for the square of M 



\Mf = ^^2ip^.p,r (4.48) 

2 4 



2A2 ^2 



with 

A = \ (4 49) 

In obtaining (4.48) we used tlie replacement rule 

Y,^^*{p)e^{p)^-g^'' (4.50) 



and the identity 
The decay width is 
and so finally 



^^''''''^.uap = -2 - 6$6l) . (4.51) 



11 ,1 

r = — — -|A/|'- (4.52) 



647r3 A2 ■ 
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4.5.3 Numerical results 

In this section we discuss some numerical results regarding the implementation of the 
decaying Majoron dark matter hypothesis in our scenario. In Ref. [90] it was shown that 
the experimental limit in the Majoron decay rate into photons is of the order of 10"'^° s~^. 
It was also shown that, in a generic seesaw model, a sizeable triplet vev plays a crucial 
role in bringing the decay rate close to this experimental bound. Here we have computed 
the width of the Majoron into neutrinos and photons in our extended seesaw model which 
incorporates the A4 flavor symmetry, generalizing the models of Ref. [46] . The results are 
shown in Fig. 4.5. These take into account the current neutrino oscillation data, discussed 
in section 4.4. We chose flve values for the triplet vev, =1 eV (turquoise), 100 eV 
(dark green), 10 keV (magenta), 1 MeV (grey) and 10 MeV (dark blue) and 100 MeV 
(black). For the right panel we consider only points that satisfy the WMAP constraint 
(4.19) indicated by the red horizontal band on the top of the left plot. 

In order to be able to probe our decaying Majoron dark matter scenario through the 
mono-energetic emission line one must be close to the present experimental limits on the 
photon decay channel, discussed in Ref. [90] and references therein. As mentioned, this 
requires the triplet vev to be sizeable, as shown on the right panel of Fig. 4.5 for the 
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mj (keV) mj (keV) 

Figure 4.5: Left panel: Fjjyi, as function of tlie Majoron mass respecting Eq. (4.19) for u/^ =1 
eV (turquoise), 100 eV (dark green), lOkeV (magenta), IMeV (grey), lOMeV (dark blue) and 
100 MeV (black). Right panel: T j^^ as function of the Majoron mass for the same values of 
the triplet vev as in the left panel. The upper orange shaded region is the excluded region from 
X-ray observations taken from Ref.[15]. 

same choices of -ua- In principle there is an additional lower bound on the Majoron mass 
coming from the Tremaine-Gunn argument [93] , which, for fermionic dark matter would 
be around 500 eV. Under certain assumptions this bound could be extended to bosons, 
and is expected to be somewhat weaker [94]. The upper orange shaded region is the 
excluded region from X-ray observations given in Ref. [15]. One should point out that, in 
this model, because of the vev seesaw relation UaU^ ~ one cannot arbitrarily take large 
values for ua to enhance Tj^^ because then the singlet vev gets correspondingly smaller 
values, hence reducing the lifetime of the Majoron to values in conflict with the WMAP 
constraint. This interplay between the CMB bounds and the detectability of the gamma 
line is illustrated in Fig. 4.5, where the dark-blue points corresponding to Ua = 10 MeV 
illustrate the experimental sensitivity to our signal. 

4.6 Conclusions 

We have studied the possibility that the seesaw model with spontaneously broken un- 
gauged lepton number may simultaneously account for the observed neutrino masses and 
mixing as well as the dark matter of the Universe. We have presented a two-texture 
structure for the neutrino mass which arises in a specific seesaw scheme implementing an 
^4 flavor symmetry. A predictive pattern of neutrino masses emerges from the interplay 
of type-I and type-II seesaw contributions, with a lower bound on the neutrinoless double 
beta decay rate, which correlates with the deviation from maximality of the atmospheric 
mixing angle ^23, as well as nearly maximal CP violation, correlated with the reactor 
angle 6*13. 
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On the other hand, assuming that associated Majoron picks up a mass due to exphcit 
lepton number violating effects that may arise, say, from quantum gravity, we showed 
how it can constitute a viable candidate for decaying dark matter, consistent with cosmic 
microwave background lifetime constraints that follow from current WMAP observations. 
We have also shown how the Higgs boson triplet, the existence of which is required by 
the consistency of the model, plays a key role in providing a test of the decaying Majoron 
dark matter hypothesis, implying the existence of a mono-energetic emission line which 
arises from the sub-leading one-loop-induced decay of the Majoron into photons. We 
also discussed the possibility of probing its existence in future X-ray observations such as 
expected in NASA's Xenia mission [95]. The presence of the type-II seesaw Higgs triplet 
would also have other particle physics implications, such as lepton flavor violating decay 
rate enhancements due to the strong renormalization effects of the triplet, quite sizeable 
in a supersymmetric model. 



Chapter 5 

Seesaw Models and Lepton Flavour Violation 



At "low" energies one cannot decide whether tree-level or loop physics generates the op- 
erator equation (1.176) responsible for neutrino masses, nor can any measurements of 
neutrino angles, phases or masses distinguish between the different tree-level seesaw re- 
alizations. Under the assumption of a pure type-I or pure type-II minimal supergravity 
seesaw mechanisms, we reconsider here the prospects for reconstructing the underlying 
high energy parameters from a combination of different measurements. Clearly, as men- 
tioned at the end of section 1.2.3, observables outside the neutrino sector are needed in 
order to ultimately learn about the high energy parameters characterizing the seesaw. If 
the CERN LHC, due to take first data, finds signs of electroweak scale supersymmetry, 
indirect insight into the high-energy world might become possible through the search for 
flavour violation effects [96, 97]. 

Starting from flavour diagonal soft supersymmetry breaking terms at some high energy 
"uniflcation" scale, flavour violation appears at lower energies due to the renormalization 
group evolution of the soft breaking parameters [98, 99]. If the seesaw mechanism is 
responsible for the observed neutrino masses, the neutrino Yukawa couplings leave their 
imprint in the slepton mass matrices as flrst shown in [100]. Potentially large LFV is 
then induced by the flavour off-diagonal structure in the Yukawa couplings required by 
the large mixing angles observed in oscillation experiments [28]. Expectations for LFV 
decays such as li Ij + 'j and k — >■ 3lj in the supersymmetric seesaw have been studied 
in [84, 101-105]. For the related process of fi — e conversion in nuclei see, for example 
[106, 107]. The potential of LHC experiments in probing the allowed seesaw parameters 
through measurements of masses and branching ratios of supersymmetric particles has 
also been discussed in Refs. [108-112]. 

In two previous studies [113, 114] it was pointed out that ratios of branching ratios 
are especially useful for learning about the unknown seesaw parameters. In [114] the case 
of type-I seesaw was discussed, whereas [113] addresses the case of seesaw type-II. For the 
type-I seesaw, there are in general too many unknown parameters that preclude making 
any deflnite predictions for LFV decays. In contrast, in the simplest type-II seesaw model 
(with only one triplet coupling to standard model leptons) neutrino mixing angles can be 
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related to ratios such as Br('r2 — > ex'j*)/Br(r2 — > fJ-Xi)- 

It has been shown that, to a good approximation, such ratios do not depend on the 
mSUGRA parameter values. However, from an experimental point of view, calculations 
of absolute event rates are needed, before ratios of different final state channels can be 
studied. In [113, 114] were taken as reference just a few benchmark mSUGRA points, for 
which we have made detailed studies. In the work [5], that is described here, we have 
calculated branching ratios and event rates over a large region of mSUGRA parameter 
space, in order to identify the maximal number of events one can expect in experiments 
at the LHC, while still respecting all low-energy constraints. 

5.1 Theory setup 

In order to fix the notation, we will briefly review the main features of the seesaw mech- 
anism and mSUGRA. As already discussed in section 1.2.3, the type-I supersymmetric 
seesaw consists in extending the particle content of the MSSM by three gauge singlet 
"right-handed" neutrino superfields. The leptonic part of the superpotential is then 

w = Yi%Hd^ + yj%h^n:; + m.n^ni j = i, . . . , 3, (5.1) 

where Y^, and Y^, denote the charged lepton and neutrino Yukawa couplings, while Af are 
the "right-handed" neutrino superfields with M; Majorana mass terms. One can always 
choose a basis on which the Majorana mass matrix of the "right-handed" neutrinos is 
brought to diagonal form Af^ = diag(Mi, M2, M3). Without loss of generality we will 
also assume that equation (5.1) is written on the basis where the charged lepton Yukawa 
matrix is already diagonal. In this simple setup, the type-I seesaw model, as defined 
by equation (5.1), is characterized by a total of 21 parameters, from which only 12 are 
measurable in the low-energy theory, as we discuss below. 

The effective mass matrix of the "left-handed" neutrinos at low energies is then given 
as (see section 1.2.3) 

m. = -'^Yj.M-^.Y,, (5.2) 

so that, for each "right-handed" neutrino, there is one non-zero eigenvalue in rriy. In 
equation (5.2) we use the notation {Hu^d) = for the vacuum expectation values of the 
neutral components of the Higgs boson doublets. 

The parameters of equation(5.1) are defined at the Grand Unified Theory (GUT) scale, 
whereas the entries of equation (5.2) are measured at low energies. In order to connect 
these two scales we numerically solve the full set of renormalization group equations 
(RGE) [104, 115]. 

As already discussed on Chapter 1, the light Majorana neutrino mass matrix in equa- 
tion (5.2), being complex symmetric, is diagonalized by a unitary 3x3 matrix U [53] as 
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shown in (1.184). Inverting the seesaw equation, equation (5.2), allows to express as 
[21] 

Y^ = V2—\^R-R-J^^-U\ (5.3) 

where rhi, is the diagonal matrix with mi eigenvalues and in general i? is a complex 
orthogonal matrix. Note that, in the special case R = 1^ Yy contains only "diagonal" 
products y'Mi'mi. In this simplified case the 18 parameters in Y^ are reduced to 12. Note 
also that in general type-I seesaw schemes, the unitary matrix diagonalizing the effective 
neutrino mass matrix differs from the lepton mixing matrix by terms of order D/Mji, 
where the D = Y^v^. For the high-scale schemes considered here one can safely neglect 
these deviations ^ . In this case we can set the diagonalization matrix as the lepton mixing 
matrix (partially) determined in neutrino oscillation measurements. 

Implementing the type-II seesaw mechanism within supersymmetry requires at least 
two SU (2) triplet states Ti^2 for the cancelation of anomalies and because of the superpo- 
tential holomorphy. A scalar triplet with mass below the GUT scale changes the running 
of gi and g2 in an unwanted way and gauge coupling unification shown in section 1.1.4 is 
lost. If instead one adds only complete SU{5) multiplets (or GUT multiplets which can 
be decomposed into complete SU{5) multiplets) to the standard model particle content, 
the scale where couplings unify remains the same (at one loop level), only the value of 
the GUT coupling itself changes [116]. 

Our numerical calculation uses an SU{5) inspired model [86, 87], which adds a pair of 
15 and 15 to the Minimal Supersymmetric Standard Model (MSSM) particle spectrum. 
This variant of the type-II seesaw mechanism, as discussed above, allows us to maintain 
gauge coupling unification even for Mt <C Mq, Mq being the unification scale. Under 
SU{?>) X 5'f/L(2) X Uy{1) the 15 decomposes as 

15 = S + T + Z (5.4) 
S ~ (6,1,-^), r~ (1,3,1), Z~ (3,2,1). 

T has the correct quantum numbers to generate the dimension-5 operator of equation 
(1.176). The SU{5) invariant superpotential reads 



1 = 1 . = = 1 - - ^ ....... . 

• Oh 



W = -=Ylf 5i ■ 15 • 5j + -^\i5h • 15 • 5h + -^Xi^n • 15 • 5^ + Yf. 10, • 5, • 5 

V2 V2 v2 

+ Yj!^ 10, • lOj • 5h + M15I5 ■ T5 + MsSh • 5h . (5.5) 
Here, 5 = (d^ L), 10 = (n^ e^ Q), 5h = {t, H2) and 5h = {t, H^). Below the GUT scale. 



^However for other type-I schemes, hke the inverse seesaw [84, 107] this approximation fails and leads 
to large LFV from right-handed neutrino exchange, even in the absence of supersymmetric contributions. 
For a systematic perturbative seesaw diagonalization method that covers all cases see Ref. [53]. 
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in the SU{5)-hroken phase, the superpotential contains the terms 

-^{Y^' L,T,L, + Y^' dfiSd^ + Y'^' d'rZL, + Y^' d^Q.H, + Y^^ u'rQ.H^ + Y:^ e^L^Ha 
+^{X^HdT^Hd + XiHuT^H^) + MtT^T^ + MzZ^Z^ + MsS,S2 + iiH^H^ . (5.6) 

As long as Mz ~ Ms ~ My, gauge couphng unification wiU be preserved. Note that 
exact equahty is not required for a successful unification. In our numerical studies we 
have taken into account the different running of these mass parameters. 

Integrating out the heavy triplets at their mass scale, the dimension-5 operator of equa- 
tion (1.176) is generated and after electroweak symmetry breaking the resulting neutrino 
mass matrix can be written as 

As in the case of the type-I seesaw, equation (5.7) depends on the energy scale. In order 
to compute the neutrino mass rrii, measured at low energies, one needs to know A2, Ir 
and Mt as input parameters at the high energy scale. As will be discussed in section 5.2, 
one can use an iterative procedure in order to find the high scale parameters from the low 
energy measured quantities. 

Note that, without loss of generality, we have the freedom to write eqs. (5.1) and (5.5) 
on the basis where the charged lepton mass matrix is diagonal, fitting the corresponding 
Yukawa couplings so as to reproduce the three measured charged lepton masses. However 
there are important differences between the type-I and type-II seesaw schemes. For ex- 
ample, in contrast to type-I, in a pure type-II seesaw scheme the unitary matrix U that 
diagonalizes equation (5.7) coincides with the lepton mixing matrix studied in neutrino 
oscillations. Moreover, in sequential type-I seesaw for each "right-handed" neutrino added 
there is one non-zero light neutrino mass eigenstate ^. In contrast, in type-II seesaw one 
can produce three neutrino masses with just one pair of triplet superfields, with only one 
triplet directly coupling to leptons. This implies that in the minimal type-II seesaw one 
has less parameters than in the sequential type-I seesaw. Indeed, as already mentioned 
in section 1.2.3, from the 12 parameters in the complex symmetric Yt matrix, one can 
remove 3 phases by redefining the charged leptons [53] . In addition, from the 3 complex 
parameters A1.2 and Mt, one does not enter, as only one of the triplets couples to leptons, 
and finally, two of the three phases can also be removed by field redefinitions. The net 
result is that there are only 11 physical parameters governing neutrino physics [113]. This 
number is substantially smaller than the 18 free parameters describing the simplest type-I 



^We do not consider here the possibility of having just two right-handed neutrino states in the type-I 
seesaw, called (3,2) in Ref. [53]. This could well account for the current neutrino data with just 12 
parameters, instead of the 18 characterizing the sequential (3,3) seesaw considered here. 
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seesaw scheme containing three "right-handed" neutrinos [117] ^. 

At low energies a maximum of 9 neutrino parameters can be fixed by measuring lepton 
properties: 3 neutrino masses, 3 mixing angles and 3 CP phases. Thus from neutrino data 
only, neither type-I nor type-II seesaw schemes can be completely reconstructed, even in 
their simplest realizations. However, especially important in the following is the fact that 
low-energy neutrino angles are directly related to the high-energy Yukawa matrix in the 
type-II seesaw, whereas no such simple connection exists in the seesaw type-I (see also 
the discussion in [118]). 

As already commented above and in Chapter 1, to a good approximation the lepton 
mixing matrix may be taken in unitary form, with three mixing angles Oij, and three 
physical CP phases (f)ij [53]. Of these only the leptonic analogue of the Kobayashi-Maskawa 
phase 6, taken as the invariant combination 6 = (pu — (pis + </)23 would enter the class of 
LFV processes discussed in this paper, so that we get the standard form (1.186). Since 
no current experiment is sensitive enough to probe leptonic CP violation we take, for 
simplicity, 6 = 0. Neutrino oscillation experiments can be fitted with either a normal 
hierarchical spectrum (NH), or with an inverted hierarchy (IH) one. If one does not 
insist in ordering the neutrino mass eigenstates m^^, i = 1, 2, 3 with respect to increasing 
mass, the matrix U can describe both possibilities without re-ordering of angles. In 
this convention, which we will use in the following, m^-^ ~ {m^.^ ~ 0) corresponds to 
normal (inverse) hierarchy and Si2, Si3 and S23 are the angles in both types of spectra. 
Basically S12 is measured in solar -|- reactor experiments, S23 in atmospheric -|- accelerator 
experiments and S13 is constrained by reactor neutrino oscillation data. 

In the general MSSM, LFV off-diagonal entries in the slepton mass matrices involve 
additional free parameters which arise from the mechanism of supersymmetry breaking. 
In order to relate LFV in the slepton sector with the LFV encoded in or Yt one must 
assume some particular scheme for supersymmetry breaking. For simplicity and definite- 
ness we will adopt mSUGRA boundary conditions, characterized by four continuous real 
and one discrete free parameter, usually denoted as 



Here, mo is the common scalar mass, M1/2 the gaugino mass and Aq the common trilinear 
parameter, all defined at the grand unification scale, Mq ~ 2 • 10^^ GeV. The remaining 
two parameters are tan (3 = v^/vd and the sign of the Higgs mixing parameter /i. 

In order to have a qualitative understanding of the magnitudes of the LFV rates we 
first present approximate leading-log analytical solutions for the renormalization group 
equations ^ . For the case of type-I seesaw, the LFV elements induced in the charged left- 

^We are treating the three charged lepton masses as experimentally determined parameters. 
^Note that in the numerical code that leads to the results presented in our plots we have numerically 
solved the full set of RGEs. 



mo, Ml/2, ^0, tan/3, Sgn(;u) . 
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slepton mass matrix by renormalization group evolution can be approximated as [104] 

(AM|),, = -^{5ml + Al){Y}LY,),, , (5.9) 

where is given in terms of the neutrino parameters by equation (5.3) and the factor L 
is defined as 

Lki= log (^^)6ki. (5.10) 

Similarly, one can get an analogous approximate expression for the off-diagonal el- 
ements of the charged left-slepton mass matrix characterizing LFV in type-II seesaw 
schemes [86]. 

5.2 Numerical results 

Due to the non-trivial structure of the neutrino Yukawa matrix Y^ in equation (5.3) 
and of Yt in equation (5.7) for type-I and type-II seesaw, respectively, the slepton mass 
matrices contain calculable LFV entries [99, 100]. In order to determine their magnitude 
we solve the complete set of renormalization group equations, given in [86, 104, 115]. All 
results presented below have been obtained with the lepton flavour violating version of 
the program package SPheno [119], where the RGBs for the MSSM part are implemented 
at the 2-loop level. For definiteness we set neutrino mass squared differences to their 
current best fit values [28] and fix the angles to the Tri-Bi-Maximal (TBM) values [29]. 

Fixing the values of other mSUGRA parameters, we used SPheno to perform a nu- 
merical scan over the mQ-Mi/2 plane. For each point in this plane, we adjust the value 
of Mji (Mt) in order to keep the low energy LFV observable BR(/i — )> 67) within its 
present experimental upper bound or within the expected sensitivity of the upcoming 
experiments [120]. 

For type-I seesaw our numerical procedure to fit these masses is as follows. As we have 
already commented, the large number of free parameters characterizing even the simplest 
type-I seesaw schemes forces us to make simplifying assumptions in inverting the seesaw 
equation, equation (5.3). As a first step we assume degenerate "right-handed" neutrinos 
and the simplest possible, flavourless, structure for the matrix /?, i.e. 

R=l, Mr,, = Mr 5ij . (5.11) 

Moreover, we fix the values of the light neutrino masses and Yukawa couplings to repro- 
duce the TBM angle values. In order to determine the resulting LFV observables we 
numerically integrate the renormalization group equations taking into account the fiavour 
structure of the Vj, matrix. We integrate out every "right-handed" neutrino and its super- 
partner at the scale associated to its mass, and calculate the corresponding contribution 
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to the dimension-five operator which is evolved to the electroweak scale. This way we 
obtain the exact neutrino masses and mixing angles for this first guess. The difference 
between the results numerically obtained and the input numbers is then minimized in an 
iterative procedure until convergence is achieved. 

For the type-II seesaw the calculations are performed for the 15-plet case, under the 
assumption Yz = Yt = Ys at Mq, as discussed above, and including the one- loop RGEs 
for the new parameters in SPheno. For consistency, we have also included 1-loop threshold 
corrections for gauge couplings and gaugino mass parameters at the scale corresponding 
to the mass of the triplet, My. The MSSM part is implemented at the 2-loop level and, 
thus, in principle one should also consistently include the effect of the 15-plets for all 
parameters at this level. However, as discussed in [113], the correct fit of the neutrino 
data requires that either the triplet (15-plet) Yukawa couplings are small and/or that Mt 
is close to Mg, implying that the ratio Mt/Mq is significantly smaller than Ma/rriz and 
thus one expects only small effects. Inverting the seesaw equation for any fixed value of 
A2 in equation (5.7), one can get a first guess of the Yukawa couplings for any fixed values 
of the light neutrino masses as a function of the corresponding triplet mass. This first 
guess will not give the correct Yukawa couplings, since the neutrino masses and mixing 
angles are measured at low energy, whereas for the calculation of we need to insert the 
parameters at the high-energy scale. However, we can use this first guess to numerically 
run the RGEs to obtain the exact neutrino parameters (at low energies) for these input 
values. The difference between the results obtained numerically and the input can then 
be minimized in a simple iterative procedure until convergence is achieved. As long as 
neutrino Yukawas are not large, convergence is reached in a few steps. However, in type-II 
seesaw schemes, the Yukawa couplings run stronger than in the type-I seesaw. Thus our 
initial guess can sizeable deviate from the exact Yukawa coupling values. Since neutrino 
oscillation data requires at least one neutrino mass to be larger than about 0.05 eV, we 
do not find any solutions for Mt > 10^^ GeV. 

Finally, the calculation of cross sections for the production of supersymmetric particles 
was done using Prospino [121-125]. The input data was taken from SPheno using the 
SUSY Les Houches Accord standard format [126]. 

5.2.1 LFV stau decays 

The eagerly awaited production of supersymmetric particles at the LHC would open new 
opportunities for the study of flavour violation in the supersymmetric sector [96] . Here 
we study how the LFV decays of staus may provide valuable cross-checks of neutrino 
properties determined at low energies as well as complementary information on the origin 
of neutrino mass. 

The expected LFV branching ratios for f2 — >■ ^ + Xi ^iid ^^2 ^ e + Xi depend on 
the choice of the mSUGRA parameters. After a full scan over the mSUGRA parameter 
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Figure 5.1: Br(f2 — > ^ + Xi) (left panel) and Br(f2 — > e + Xi) (right panel), in the mo,Mi/2 
plane for our standard choice of parameters: fi > 0, tan jS — 10 and Aq = GeV, for type-I 
seesaw, imposing Br(/z — > e + 7) < 1.2 • 10~^^. 
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Figure 5.2: Br(f2 — > ^ + x?) (left panel) and Br(f2 — >• e + Xi) (right panel), in the mo,Mi/2 
plane for om' standard choice of parameters: fi > 0, tan jS — 10 and Aq = GeV, for type-I 
seesaw, imposing Br(/i —7-6-1-7) ^ 10~^^. 

space we found that the dependence on Aq and on the sign of /i is weaker, but that the 
rates decreased with increasing values of tan /3. Therefore, we chose our standard point 
with a relatively low value of tan (3 = 10, and for definiteness took /i > 0, and Aq = 0. 
In Fig. 5.1 we show the contour plots for the LFV decays f2^fi + Xi (l^ft panel) and 
f2 — >■ e -h Xi (right panel) in the mo,Mi/2 plane for our standard choice of mSUGRA 
parameters for the simplest pure type-I seesaw scheme. One sees that there are regions 
in parameter space where the LFV decays of the f2 can be as large as of order 10~^. In 
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these plots the values of ilf^ were chosen as to obtain the maximum LFV compatible 
with the present experimental limit of Br(/i — )> e + 7) < 1.2 • 10~^^ [43]. Also shown in 
these plots are the exclusion regions coming from the LEP constraints on SUSY masses 
and also the exclusion obtained when the neutralino is not the LSP ^. In Fig. 5.2 we 
show the same contour plots for Br(/_f — >■ e + 7) < 10~^^, which will be achievable in the 
coming experiments [120]. Also in this case one observes in Fig. 5.2 that the LFV stau 
decay rates may exceed the 10% level. Notice also that the nontrivial features present 
in both Figs. 5.1 and Fig. 5.2 reflect the well-known cancellations between chargino and 
neutralino contributions to — )> e + 7 already discussed above. 
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Figure 5.3: Br(f2 — > /j, + Xi) (left panel) and Br(f2 — > e + Xi) (right panel), in the mo, Afi/2 
plane for Ai — 0.02 and A2 = 0.5 and our standard choice of parameters: /i > 0, tan /S — 10 and 
Aq — Q GeV, for type-II seesaw, imposing Br(/i — )• e + 7) < 1.2 ■ 10~^^. 



In Fig. 5.3 and Fig. 5.4 the same type of plots are shown for type-II seesaw. A 
comparison of these figures shows that, qualitatively, the behavior is very similar for 
the two types of seesaw. In both cases, the larger rates for f2 — >■ e -|- Xi more 
constrained in parameter space than those for f2 — ^J^+Xv Notice however that there is an 
important difference between type-I and type-II seesaw, coming from the presence of the 
Higgs triplets that contribute sizeably to the running of the type-II beta functions. This 
gets reflected in the supersymmetric particle spectra and hence in the shapes of the red 
(shaded) regions in Fig. 5.3 and Fig. 5.4. One can observe, indeed, that the regions where 
the stau is the lightest supersymmetric particle, as well as the regions already excluded 
by LEP2 are substantially different for type-II seesaw, as compared to the corresponding 
ones for type-I. This follows from the modification in the beta functions introduced by the 
addition of the Higgs triplets, making Mi and M2 smaller in type-II than in type-I seesaw 

^Note that we did not display the constraints coming from Dark Matter (DM) rehc abundance. 
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Figure 5.4: Br(f2 — > ^ + Xi) (left panel) and Br(f2 — > e + Xi) (right panel), in the mo,Mi/2 
plane, for Ai = 0.02 and A2 = 0.5 and our standard choice of parameters: /J. > 0, tan /3 = 10 and 
Aq — Q GeV, for type-II seesaw, imposing Br(/i — )• e + 7) < 10~^^. 



for the same value of Mi/2. The variation with the mSUGRA parameters is illustrated 
in Fig. 5.5 (type-I) and Fig. 5.6 (type-II) for the parameter Aq and in Fig. 5.7 (type-I) 
and Fig. 5.8 (type-II) for tan/?. We can see that there is not much variation with Aq, 
while the rates decrease rapidly with increasing values of tan/3. The reason for this is 
that BR(/i — > e -|- 7) increases along with tan^ /3, thus constraining more strongly the 
maximum attainable stau LFV rates. This effect is stronger for type-I as can be seen by 
noting the different values for the contour levels in Fig. 5.7 and Fig. 5.8. The variation 
with the sign of /i is weak and we do not show it here. So, in summary, large LFV rates 
prefer moderate values of tan jS and this explains a posteriori the choice of our standard 
parameters. 

5.2.2 Total production cross section of X2 

As important as having a large branching ratio into a LFV final state, it is to be able 
to produce a large enough event sample. In order to estimate the number of LFV events 
expected at the LHC, one notes that, from Figs. 5.1 - 5.8, in the regions where the LFV 
is sizeable, the direct production of staus at the LHC is negligible compared to that 
which arises from cascade decays of heavier neutralinos, mainly X2- We focus on the 
X2, because decays such as X2 ~^ f^T~Xi ^^'^ sensitive to flavour violation, whereas in the 
corresponding chargino decays the flavour information is lost. Hence we flrst compute 
the total X2 production cross section. In the left panel of Fig. 5.9 we show the results 
for the cross- section for X2 production as a function of Afi/2, for different choices of niQ 
and for our standard choice of mSUGRA parameters: /i > 0, tan /3 = 10 and ^40 = 
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Figure 5.5: Br(f2 — > + Xi) (left panel) and Br(f2 — > e + x^) (right panel), in the mo,Mi/2 
plane for standard choice of parameters: /i > 0, tan 13 — 10 but different Aq = —300 GeV, for 
type-I seesaw, imposing Br(/z — > e + 7) < 1.2 ■ 10"-^^. 
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Figure 5.6: Br(f2 — > + x?) (left panel) and Br(f2 — )• e + Xi) (right panel), in the mo,Mi/2 
plane for Ai = 0.02 and A2 — 0.5 and standard choice of parameters: /.t > 0, tan /3 = 10 but 
different Aq = —300 GeV, for type-II seesaw, imposing Bv{^ — )• e + 7) < 1.2 • 10~^^. 



GeV, for the pure type-I mSUGRA seesaw scheme. This choice of mSUGRA parameters 
corresponds, as will be discussed below, to the case where the branching ratios of the LFV 
stau decays are the largest. This result was obtained using the Prospino code [121-125] 
at Leading Order (LO) approximation. We have checked that the Next to Leading Order 
(NLO) calculation only changes the results slightly, due to an appropriate choice of the 
renormalization scale [121-125]. So, in all cross sections presented here, we only used the 
LO approximation. The corresponding results for type-II seesaw are shown in the right 
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Figure 5.7: Br(f2 — > + Xi) (^^^^ panel) and Br(f2 — > e + Xi) (right panel), in the mo,Mi/2 
plane for standard choice of parameters: /i > 0, Aq = but different tan/3 — 30, for type-I 
seesaw, imposing Br(/z — > e + 7) < 1.2 • 10~^^. 
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Figure 5.8: Br(f2 + Xi) (left panel) and Br(f2 — )■ e + Xi) (right panel), for Ai = 0.02 and 
A2 = 0.5, in the mo, M1/2 plane for standard choice of parameters: > 0, Aq = 0, but different 
tan/3 — 30, for type-II seesaw, imposing Br(// — s- e + 7) < 1.2 • 10~^^. 



panel of Fig. 5.9, for the same choice of mSUGRA parameters and for Ai = 0.02 and 
A2 = 0.5. 



5.2.3 Total production of X2 times BR to /-l-t lepton pair 

In order to get an estimate of the expected number of LFV events at the LHC we now 
use a combination of the Prospino and SPheno codes to evaluate the product of the X2 
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Figure 5.9: Production cross section (at leading order) of X2 versus for varying mg, and 

for our standard choice of parameters: > 0, tan /3 = 10 and ^0 = GeV, in type-I seesaw (left 
panel) and type-II seesaw (right panel) for Ai — 0.02 and A2 ~ 0.5. 



production cross section times the branching ratios into LFV processes. Once we know 
the luminosity at LHC we can multiply it by the above product to get the number of 
events. 

In Fig. 5.10, we have plotted, for type-I seesaw (left panel) and type-II (right panel), 
the production cross section at leading order of the second lightest neutralino cr(x2) times 
the BR of X2 goii^g to the opposite-sign dilepton signal x'j'yur as a function of M1/2, for 
different values of niQ. We have fixed the rest of the mSUGRA parameters to our standard 
mSUGRA point and imposed an upper limit on Br(/i — > e -|- 7) < 1.2 • 10~^^. In type-I 
seesaw, the number of events of the opposite-sign dilepton signal X2 ~^ Xi ^^'^ can be of 
the order of 10"^ for mo ~ 100 GeV and M1/2 ~ [450, 600] GeV, assuming a luminosity 
C = 100 fb~^ and = 14 TeV. In type-II seesaw, there can be a maximum number of 
events of the order of 10^ for mo ~ 100 GeV and M1/2 ~ [600, 800] GeV. 




Figure 5.10: Production cross section (at leading order) of X2 times BR of X2 goi^iS to /i-T 
lepton pair versus A/1/2 for mo = 100 GeV (red), 200 GeV (green), 300 GeV (blue) and 500 GeV 
(magenta), and for our standard choice of parameters: fJ, > 0, tan /? = 10 and ^Iq = GeV, for 
type-I (left panel) and for type-II seesaw (right panel) with Ai — 0.02 and A2 = 0.5, imposing 
Br(^^ e-h7) < 1.2 ■ 10"^^ 
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For type-II seesaw where we have less parameters, we can look at variations of the 
result with the values of the triplet Higgs boson coupling A2, a parameter that can not 
be determined from neutrino data alone as it appears only in the ratio A2/MT, see equa- 
tion (5.7). In Fig. 5.11 we show the dependence of the product of cross section times LFV 
branching ratios as function of A2 for our standard point. We should mention that the 
other Higgs boson triplet coupling Ai, does not contribute to LFV decays, and hence is 
left undetermined by this analysis. 
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Figure 5.11: Production cross section (at leading order) of X2 times BR of X2 goii^S to fi-r 
lepton pair versus M1/2, for our standard choice of parameters: > 0, tan /3 = 10 and Aq — 



GeV, for type-II seesaw, imposing Br(yU —7-6-1-7) < 1.2 -10 , for a fixed value of mo = 100 
GeV and different values of A2 = 0.1 (green), 0.5 (red), 0.9 (blue). 



As has been discussed in [127], the dominant standard model backgrounds for the 
process considered are expected to be WW and tt production. The cuts necessary to 
reduce this background will depend on the details of the SUSY spectrum and a detailed 
investigation is beyond the scope of this paper. The results of [127] suggest that the signal 
should be visible for a(x2)xBR of order 0(10) fb. 



5.3 Conclusions 



Low energy neutrino experiments, including oscillation studies and neutrinoless double- 
beta decay searches, may optimistically determine at most 9 neutrino parameters: the 3 
neutrino masses, the 3 mixing angles and potentially the 3 CP violating phases. This is 
insufficient to fully reconstruct the underlying mechanism of neutrino mass generation. 
Under the assumption that neutrino masses arise a la seesaw, we have considered the 
simplest pure type-I or pure type-II seesaw schemes in mSUGRA. 

We have performed a full scan over the mSUGRA parameter space in order to identify 
regions where LFV decays of X2 ^^^"^ be maximal, while still respecting low-energy con- 
straints that follow from the upper bounds on Br(/i — )> 67). We have also estimated the 
expected number of events for X2 ~^ Xi + + foi" ^ sample luminosity of C = 100 fb~^. 
The expected number of events for the other channel X2 ~^ Xi + + ^ is always smaller. 
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as can be seen from the LVF branching ratios presented in section 5.2.1. We have found 
that the pure seesaw-II scheme is substantiahy simpler and comes closer to beeing fully 
reconstructable, provided additional LFV decays are detected and some supersymmetric 
particles are discovered at the Large Hadron Collider, providing the necessary parameter 
reconstruction information of the supersymmetric lagrangian. 

Note that in what concerns the expected number of events both type-I and type- 
II schemes give similar results. However, as we have seen, given their smaller number 
of paramaters, type-II seesaw schemes are more likely to be reconstructable through a 
combination of low energy neutrino measurements, with the possible detection of super- 
symmetric states and lepton flavour violation at the LHC. This should encourage one to 
perform full-fledged dedicated simulations, in order to ascertain their feasibility within 
realistic experimental conditions [96]. 

Finally we note that we have not analyzed in detail the fact that LFV might induce 
new "edge variables", giving additional information [128]. We have focused here on LHC, 
but we should mention that a future ILC would be much more suited for measuring LFV 
SUSY processes [112, 129-134]. 



Chapter 6 

Dark Matter in Seesaw Type II Model 



Standard cosmology requires the existence of a non-baryonic dark matter (DM) contribu- 
tion to the total energy budget of the universe [135, 136]. In the past few years estimates 
of the DM abundance have become increasingly precise. Indeed, the Particle Data Group 
now quotes at 1 o" c.l. [43] 

= 0.105 ±0.008. (6.1) 

Since the data from the WMAP satellite [81, 137] and large scale structure formation [138] 
is best fitted if the DM is cold, weakly interacting mass particles (WIMP) are currently the 
preferred explanation, as discussed in section 2.4.3 - see formulas (2.186), (2.188) and the 
comment that follows. While there is certainly no shortage of WIMP candidates (lists can 
be found in many reviews, see for example [135, 136, 139, 140]), the literature is completely 
dominated by studies of the lightest neutralino. Neutrino oscillation experiments have 
shown that neutrinos have non-zero mass and mixing angles [10, 58, 59, 141, 142] (see 
the discussion in Chapter 1) and the most recent global fits to all data [28] confirm 
again that the mixing angles are surprisingly close to the so-called tri-bimaximal mixing 
(TBM) values [29], as in equation (1.228). In the minimal supersymmetric extension of 
the standard model (MSSM) with conserved R-parity neutrino masses are zero for the 
same reasons as in the SM. However, it was shown long ago that if neutrinos are Majorana 
particles, their mass is described by a unique dimension-5 operator [19] given by equation 
(1.176). All (Majorana) neutrino mass models reduce to this operator at low energies. 
If / is a coefficient 0(1), current neutrino data indicates A < 0(10^^) GeV. This is the 
essence of the "seesaw" mechanism. There are three different tree-level realizations of the 
seesaw, classified as type-I, type-II and type-Ill in [143]. As already discussed on Chapter 
1, Type-I is the well-known case of the exchange of a heavy fermionic singlet [70, 144-146] 
and Type-II corresponds to the exchange of a scalar triplet [53, 77]. One could also add 
one (or more) fermionic triplets to the field content of the SM [147] . This is called seesaw 
type-Ill in [143]. 

Neutrino experiments at low energies measure only fa/s/A, thus observables outside 
the neutrino sector will ultimately be needed to learn about the origin of equation (1.176). 



126 



Dark Matter in Seesaw Type II Model 



Augmenting the SM with a high-scale seesaw mechanism does not lead to any conceiv- 
able phenomenology apart from neutrino masses, but if weak scale supersymmetry exists 
indirect probes into the high energy world might be possible. Two kind of measurements 
containing such indirect information exist in principle, lepton flavour violating (LFV) 
observables and sparticle masses. 

Assuming complete flavour blindness in the soft supersymmetry breaking parameters 
at some large scale, the neutrino Yukawa matrices will, in general, lead to non-zero flavour 
violating entries in the slepton mass matrices, if the seesaw scale is lower than the scale at 
which SUSY is broken. This was first pointed out in [100]. The resulting LFV processes 
have been studied in many publications, for low-energy observables such as /i — >■ 67 
and /i — e conversion in seesaw type-I see for example [101-106, 118], for seesaw type-II 
[86, 113]. LFV collider observables have also been studied in a number of papers, see for 
example [5, 96, 112-114, 127-130, 132-134, 148]. 

Mass measurements in the sparticle sector will not only be necessary to learn about the 
mechanism of SUSY breaking in general, but might also reveal indications about the scale 
of the seesaw mechanism. However, very precise knowledge of masses will be necessary 
before one can learn about the high scale parameters [108, 111]. Especially interesting in 
this context is the observation that from the different soft scalar and gaugino masses one 
can define certain combinations ( "invariants" ) which are nearly constant over large parts of 
mSugra space. Adding a seesaw mechanism of type-II or type-Ill these invariants change 
in a characteristic way as a function of the seesaw scale and are thus especially suited 
to extract information about the high energy parameters [109]. Note, however, that the 
"invariants" are constants in mSugra space only in leading order and that quantitatively 
important 2-loop corrections exist [113]. 

In the work [6], that is reproduced here, we study neutralino dark matter [149-151] 
within a supersymmetric type-II seesaw model with mSugra boundary conditions. For 
definiteness, the model we consider consists of the MSSM particle spectrum to which 
we add a single pair of 15- and 15-plets. This is the simplest supersymmetric type-II 
setup, which allows one to maintain gauge coupling unification [86] and explains measured 
neutrino oscillation data. 

In mSugra - assuming a standard thermal history of the early universe ^ - only four very 
specific regions in parameter space can correctly explain the most recent WMAP data [81]. 
As described in section 2.4.3, these are (i) the bulk region; (ii) the co-annihilation line; 
(iii) the "focus point" line and (iv) the "Higgs funnel" region. In the bulk region there are 
no specific relations among the sparticle masses. However, all sparticles are rather light in 
this region, so it is already very constrained from the view point of low-energy data [153]. 
In the co-annihilation line the lightest scalar tau is nearly degenerate with the lightest 

^In models with non-standard thermal history the relation between sparticle masses and relic density 
can be lost completely [152]. 
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neutralino, thus reducing the neutrahno rehc density with respect to naive expectations 
[151, 154]. In the "focus point" hne [154, 155] ^xl^^ small enough to explain QjjAih'^ 
due to a rather small value of fi leading to an enhanced higgsino component in the lightest 
neutralino and thus an enhanced coupling to the boson. Lastly, at large tan j3 an s- 
channel resonance pair annihilation of neutralinos through the CP-odd Higgs boson can 
become important. This is called the "Higgs funnel" region [149]. 

The addition of the 15 and 15 pair at the high scale does not, in general, lead to 
the appearance of new allowed regions. However, the deformed sparticle spectrum with 
respect to mSugra expectations leads to characteristic changes in the allowed regions as 
a function of the unknown seesaw scale. We discuss these changes in detail and compare 
the results to other indirect constraints, namely, the observed neutrino masses and upper 
limits on LFV processes. We concentrate on the seesaw type-II scheme, since for mSugra 
-|- seesaw type-I the changes in the DM allowed regions with respect to pure mSugra are, 
in general, expected to be tiny. ^ 



6.1 Theory setup: mSugra and SU (5) motivated type-II seesaw 

In this section we summarize the main features of the model we will use in the numerical 
calculation. We will always refer to minimal Supergravity (mSugra) as the "standard" 
against which we compare all our results. The model consists in extending the MSSM 
particle spectrum by a pair of 15 and 15. It is the minimal supersymmetric seesaw type-II 
model which maintains gauge coupling unification [86] . 

mSugra is specified by 4 continuous and one discrete parameter [157]. These are 
usually chosen to be mo, the common scalar mass, M1/2, the gaugino mass parameter, 
Aq, the common trilinear parameter, tan/3 = ^ and the sign of fi. mo, M1/2 and Aq are 



defined at the GUT scale, the RGBs are known at the 2-loop level [158] 
Under SU{3) x SUl{2) x U{1)y the 15 decomposes as 

15 = S + T + Z (6.2) 

S ~ (6,1,-^), r~ (1,3,1), Z~(3,2,^). 

3 o 

The SU (5) invariant superpotential reads as 

W = ^Yi55-15-5 + ^Ai5jf-15-5H + ^A25^f-T5-5H + Y5l0-5-5^(6.3) 
v2 v2 v2 

+ YiolO- 10-5h + Mi515- I5 + M55 H ■ 5h 
Here, 5 = (d^ L), 10 = (n^ Q), 5h = (t, H2) and 5h = (t, Hi). Below the GUT scale 



^We have confirmed this general expectation with some sample calculations. However, an exceptional 
case has been presented recently in [156], see the more detailed discussion in section (6.2). 
^For reviews on mSugra and MSSM, see for example [157, 159, 160]. 
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in the SU (5)-broken phase the potential contains the terms 

-^{YtLTiL + Ysd'Sd') + Yz(fZL + YdcfQH^ + Y^u'QH2 + Y.e'LHi (6.4) 
v2 

+ -^{XiHiTiHi + X2H2T2H2) + MtTiTs + MzZiZs + MsSiS2 + IXH1H2 
v2 

Yd, Yu and Y^ generate quark and charged lepton masses in the usual manner. In addition 
there are the matrices Fy, Ys and Yz- For the case of a complete 15, apart from calculable 
threshold corrections, Yt = Ys = Yz and Mt, Ms and Mz are determined from M15 by the 
RGEs. As long as Mz ~ Ms ~ Mt ~ M15 gauge coupling unification will be maintained. 
The equality need not be exact for successful unification. 

The triplet Ti has the correct quantum numbers to generate neutrino masses via the 
first term in equation (6.4). Integrating out the heavy triplets at their mass scale a 
dimension-5 operator of the form equation (1.176) is generated. This can be seen as 
follows: first we compute the F term associated with T2 which in abbreviated notation is 

Ft, = -^\2H2H2 + MtTi. (6.5) 
v2 

Next we compute the corresponding scalar potential term 



F = ... + |Fr/ (6.6) 

+ h^2\^\H2f + M^\T,f + ^Re{\2MTTiH2H2) 
2 v2 



where in this expression Ti , Hi and H2 refer to the scalar components of the respective 
superfields. Then we assume that at this low energy scales the triplet is almost at rest 
such that we can ignore the kinetic term in Lagrange equations for Ti, leaving us with 

which, solving it in order to Ti gives 

Finally, the leptonic Yukawa couplings are computed from the superpotential term Wl = 
j^YrLTiL: 

2 d'^Wr 
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and after replacing (6.8) we end up with the dimension 5 operator 



LniLL = ^LH2H2L. (6.10) 



After electro- weak symmetry breaking the resulting neutrino mass matrix can be writ- 
ten as 

Here V2 is the vacuum expectation value of Higgs doublet H2 and we use the convention 
(Hi) = can be diagonalized in the standard way with a unitary matrix U, 

containing in general 3 angles and 3 phases. Note that Yt = ■ Yt ■ U is diagonalized 
by the same matrix as m^. This means that if all neutrino eigenvalues, angles and phases 
were known, Yt would be completely fixed up to an overall constant, which can be written 
as ^ ~ lO^^GeV ^ ^ "^ j . Thus, current neutrino data requires Mt to be lower than 
the GUT scale by (at least) an order or magnitude. 



The full set of RGEs for the 15 + 15 can be found in [86] and in the numerical 
calculation, presented in the next section, we solve the exact RGEs. However, for a 
qualitative understanding of the results, the following approximative solutions are quite 
helpful. 



For the gaugino masses one finds in leading order 



Mi{msusY) = TTT^^y^- 

a[Mo) 



(6.12) 



Eq. (6.12) implies that the ratio M2/M1, which is measured at low-energies, has the usual 
mSugra value, but the relationship to M1/2 is changed. Neglecting the Yukawa couplings 
Yi5 (see below), for the soft mass parameters of the first two generations one gets 



fl 



i=i 



ai{Mj 



i 1 



1 

bi + A6, 



1 + 



a{MG 



An 



b, log 



Af2 



^(6, + A6,)log^ 
47r ^ * 



(6.13) 



(6.14) 
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The various coefficients c{ can be found in [113]. Tlie gauge couplings are given as 



ai{mz) 
aiimsusv) 



3 cos^ 6w ' 

ai{mz) 



a2{mz) 



sin^ 9w 



(6.15) 



4iT i ° m-^ 



usy) 



^SUSY 



4^ -bt log ; 



47r 



(6, + A6,)log§ 



with bf^ and bf^^^^ being the usual standard model and MSSM coefficients. A6,; = 7 
for all i in case of a complete 15-plet. 

We can estimate the soft mass parameters given the above formulas for a given choice 
of mo, Ml/2 and Mi^ = Mt- We show some arbitrarily chosen examples in fig. (6.1). 
Note that the result shown is approximate, since we are (a) using the leading log ap- 
proximation and (b) two loop effects are numerically important, especially for mg, but 
not included. The figure serves to show that for any Afis < Mqut the resulting mass 
parameters are always smaller than the mSugra expectations for the same choice of initial 
parameters (mo,Afi/2)- While the exact values depend on (mo,Mi/2) and on the other 
mSugra parameters, this feature is quite generally true in all of the (mo,Mi/2) plane. 
Note, that the running is different for the different scalar mass parameters, but the ratio 
of the gaugino mass parameters M1/M2 always stays close to the mSugra expectation, 
Ml ~ |tan2 6'vi/M2. 
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Figure 6.1: Analytically calculated running of scalar (to the left) and gaugino mass parameters 
(to the right), leading order only. The mass parameters are calculated as a function of M15 for 



the mSugra parameters ttiq — 70 GeV and Mj 



1/2 



250 GeV. For M15 2 x 10^*' GeV the 



mSugra values are recovered. Smaller AI15 lead to smaller soft masses in all cases. Note that 
the running is different for the different mass parameters with gaugino masses running faster 
than slepton mass parameters. 
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6.2 Numerical results 

In this section we discuss our numerical results. All the plots shown below are based on 
the program packages SPheno [119] and micrOMEGAs [161, 162]. We use SPheno V3 [7], 
including the RGEs for the 15 + 15 case [86, 113] at the 2-loop level for gauge couplings 
and gaugino masses and at one-loop level for the remaining MSSM parameters and the 
15-plet parameters, for a discussion see [113]. For any given set of mSugra and 15-plet 
parameters SPheno calculates the supersymmetric particle spectrum at the electro-weak 
scale, which is then interfaced with micrOMEGAs2.2 [163] to calculate the relic density 
of the lightest neutralino, Q^oh"^. 

For the standard model parameters we use the PDG 2008 values [43], unless specified 
otherwise. As discussed below, especially important are the values (and errors) of the 
bottom and top quark masses, mi, = 4.2 -|- 0.17 — 0.07 GeV and rrit = 171.2 ± 2.1 GeV. 
Note, the rrit is understood to be the pole-mass and mb{mi,) is the MS mass. As the 
allowed range for Vtohih'^ we always use the 3 a c.l. boundaries as given in [43], i.e. 
VLoAih^ = [0.081,0.129]. Note, however, that the use of 1 (J contours results in very 
similar plots, due to the small error bars. 

In the "seesaw sector" we have the parameters connected with the 15-plets, i.e. M15, 
1x5, Ai and A2. For the calculation of the dark matter abundance the most important pa- 
rameter is M15. It has turned out that the effects of Fis, Ai and A2 on the relic abundance 
of neutralinos are very minor. Note, however, that as discussed in the previous section, 
atmospheric neutrino oscillation data can not be explained in our setup, if the triplet 
mass is larger than approximately Afis = Mt = 10^^ GeV. Also, the non-observation of 
lepton flavour violating (LFV) decays puts an upper bound on M^^. The latter, however, 
is strongly dependent on tan/3 and depends also on rriQ and Mi/2- We will first show 
results using different values of Mt as free parameter, without paying attention to neu- 
trino masses and LFV. We will discuss how our results change for correctly fitted neutrino 
masses and angles towards the end of this section, where we also discuss and compare 
LFV excluded regions with DM allowed ones. 

We define our "standard choice" of mSugra parameters as tan/3 = 10, Aq = Q and 
fi > and use these values in all plots, unless specified otherwise. We then show our results 
in the plane of the remaining two free parameters, (mo, M1/2). Fig. (6.2) shows in the top 
panel contours of equal dark matter density, fl^oh^. The lines are constant ^x'l^^ with 
Vt^ah? = 0.1, 0.2, 0.5, 1,2. In the bottom panel we show the range of parameters allowed 
by the DM constraint at 3 o" c.l. In both cases, to the left a pure mSugra calculation, 
whereas the plot to the right shows mSugra -|- 15-plet with Mt = 10"^^ GeV. In each plot 
the yellow regions are eluded either by the lighter scalar tau being the LSP (to the bottom 
right) or by the LEP limit on the mass of the lighter chargino (to the left), > 105 
GeV. In addition, we show two lines of constant lightest Higgs boson mass, niho = 110 
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Figure 6.2: Top: Contours of equal dark matter density (i7^o/i^) in tlie {ttiq, M1/2) plane for 
the "standard choice" tan (3 — 10, Aq — and /x > 0, for mSugra (left panel) and type-II seesaw 
with Mt = 10^'' GeV (right panel). The lines are constant il^oh"^ with Qnh'^ ^ 0.1,0.2,0.5, 1, 2. 
Bottom: Range of parameters allowed by the DM constraint at 3 cr c.l. To the left: mSugra; to 
the right: Mt — 10^^ GeV. For a discussion see text. 

GeV (dotted) and rrifio = 114.4 GeV (dashed), as calculated by SPheno, see the discussion 
below. 

The plots show three of the different allowed regions discussed in the introduction. To 
the right the co-annihilation region, here the lightest neutralino and the lighter scalar tau 
are nearly degenerate in mass. The line going nearly vertically upwards at constant M1/2 
is the "focus point" line. The small region connecting the two lines are the remains of the 
bulk region, which has shrunk considerably due to the reduced error bars on QoAih^ after 
the most recent WMAP data [81]. The focus point line is excluded by the LEP constraint 
on the lighter chargino mass at low and moderate values of ruQ. It becomes allowed only 
at values of ttiq larger than (very roughly) 1-1.5 TeV. However, note that the exact value 
of mo at which the focus point line becomes allowed is extremely sensitive to errors in 
m + , both from the experimental bound and the error in the theoretical calculation. 
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Comparing the results for the pure mSugra case to the mSugra+15-plet calculation, 
two differences are immediately visible in fig. (6.2). First, the focus point line is shifted 
towards larger values of M1/2. This is due to the fact that for the 15-plet at M15 = 10^^ 
GeV the neutralino is lighter than in the mSugra case at the same value of M1/2, compare 
to fig. (6.1). Maintaining the same relation between Mi and fi as in the mSugra case 
requires a then a larger value of Mi/2- Note that for the same reason the excluded region 
from the LEP bound on the chargino mass is larger than in the mSugra case. Second one 
finds that the co-annihilation line is shifted towards smaller values of rriQ. The latter can 
be understood from fig. (6.3). 
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Figure 6.3: Allowed region for dark matter density (0.081 < Vi^oh^ < 0.129) in the (mo,Mi/2) 
plane for the "standard clioice" tan/3 = 10, Aq = and fJ, > 0, for five values from Mt, 
Mt = 10^^ GeV (red), to Mt = 10^^ GeV (cyan), to the left. To the right: Variation of the 

(top lines) and of J^/i^ (bottom lines), as a function of Mt for four 



mass difference m^^ 



m 



xo 



different values of mo: (cyan), 50 (magenta), 100 (blue) and 150 GeV (green) for one fixed 
value of Ml/2 = 800 GeV. The yellow region corresponds to the experimentally allowed DM 



region. 



Fig. (6.3) shows the allowed region for the dark matter density in the (mo, M1/2) plane 
for our "standard choice" of other mSugra parameters for a number of different Mt (to 
the left). The plot shows how the co-annihilation line moves towards smaller values of 
niQ for smaller values of Mt- The plot on the right in fig. (6.3) explains this behaviour. 
It shows the variation of the mass difference m^j — m^^ (top lines) and of Qh^ (bottom 
lines), as a function of Mt for four different values of mo: (cyan), 50 (magenta), 100 
(blue) and 150 GeV (green) for one fixed value of M1/2 = 800 GeV. The yellow region 
corresponds to the experimentally allowed DM region. Co-annihilation requires a small 
value of rrif^ — m^^, typically smaller than a few GeV. With decreasing values of Mt the 
gaugino masses run down to smaller values faster than the slepton masses, thus effectively 
increasing rrif^ — m^g in these examples with respect to mSugra. To compensate for this 
effect at constant M1/2 smaller values of mo are required to get the ruf^ — m^^ in the 
required range. 



134 



Dark Matter in Seesaw Type II Model 




200 400 600 800 1000 200 400 600 800 1000 

My2 (GeV) M„2 (GeV) 



Figure 6.4: Limits for mSugra with tan/3 = 10, and > for Aq — —300 GeV (left panel) 
and Aq — —500 GeV (rigiit panel). The blue regions are allowed by the DM constraint, for the 
explanation of the bounds see fig. (6.2) and text. 

At this point a short discussion of the Higgs boson mass bound might be in order. 
LEP excluded a light Higgs boson with SM couplings with masses below rrih < 114.4 
GeV [43]. For reduced coupling of the Higgs boson to bb the bound is less severe, so this 
bound is not strictly valid in all of MSSM space. More important for us, however, is 
the theoretical uncertainty in the calculation of the lightest Higgs boson mass. SPheno 
calculates m/jo at two-loop level using DR renormalization. Expected errors for this kind 
of calculation, including a comparison of different public codes, have been discussed in 
[164]. As discussed in [164, 165] even at the 2-loop level uncertainties in the calculation 
of niho can be of the order of 3 — 5 GeV. In this context it is interesting to note that 
FeynHiggs [166], which calculates the Higgs masses in a diagrammatic approach within 
the OS renormalization scheme tends to predict Higgs masses which are systematically 
larger by 3 — 4 GeV, when compared with the DR calculation. We therefore showed in 
fig. (6.2) two lines of constant Higgs boson masses. The value of ruho = 114.4 GeV is 
taking the LEP bound at face value, while the lower value of m/jo = 110 GeV estimates 
the parameter region which is excluded conservatively, including the theoretical error. 
Since the lightest Higgs boson mass varies slowly with rriQ and M1/2, even a relatively tiny 
change in rriho of, say 1 GeV, shifts the extreme values of the excluded region by ~ 50 
GeV in M1/2 (at small ttiq) and by ~ 150 GeV in niQ (at small Afi/2). 

Moreover, it is well known that the calculated Higgs boson masses are strongly de- 
pendent on the mixing in the stop sector and thus, indirectly, on the value of Aq. This is 
shown for the case of a pure mSugra calculation in fig. (6.4). Here we show two examples 
for the DM allowed region and the regions disfavoured by the Higgs boson mass bound 
at mfto = 114.4 GeV and niho = 110 GeV. Larger negative Aq leads to a less stringent 
constraint (for /i > 0). Note, that all of the bulk region becomes allowed at Aq = —500 
GeV, once the theoretical uncertainty in the Higgs boson mass calculation is taken into 
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account. We have checked for a few values of Mt that for the case of mSugra+15 the 
resulting Higgs boson bounds are very similar. We thus do not repeat the corresponding 
plots here. Comparing the calculations shown in fig. (6.4) and the mSugra calculation 
in fig. (6.2) with each other, one finds that the DM allowed regions are actually affected 
very little by the choice of Aq. We have checked that this is also the case for mSugra + 
seesaw type-II. 
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Figure 6.5: Logarithmically scaled zoom into the focus point region. In red the allowed region 
for 0.081 < Qh^ < 0.129 and in cyan the allowed region due the variation of mtop — 171.2 ± 2.1 
GeV. The left panel is for mSugra case and the right panel for Mt = 10^^ GeV. The other 
parameters are taken at our "standard" values. 

As mentioned above the uncertainty in the top mass is important for the calculation 
of the relic density. At low and moderate values of tan j3 the exact value of mt affects 
mainly the focus point region. As fig. (6.2) demonstrates near the focus point line 
the relic density changes very abruptly even for tiny changes of M1/2. This is because 
a comparatively small value of fi is required to get a sufficiently enhanced coupling of 
the neutralino to the boson. In mSugra the value of /i is determined from all other 
parameters by the condition of having correct electro-weak symmetry breaking (EWSB) 
and usually leads to Mi, M2 <C fJ,- In the focus point region fi varies abruptly, points to 
the "left" of the focus point region are usually ruled out by the fact that EWSB can not 
be achieved. Since rrit is the largest fermion mass, its exact value influences the value 
of /i required to achieve EWSB most. The change of fi with respect to a change of rrit 
then can lead to a significant shift in the DM allowed region of parameter space. This 
is demonstrated in fig. (6.5), which shows a zoom into the focus point region for pure 
mSugra (to the left) and mSugra + 15 (to the right). The variation of the top mass 
shown corresponds to the current 1 a allowed range [43]. The pure mSugra is especially 
sensitive to a change of nit . At large values of rriQ the uncertainty in "fixing" M1/2 from 
the DM constraint can be larger than 100 GeV in the case of mSugra. Given this large 
uncertainty it would be impossible at present to distinguish the pure mSugra case from 
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mSugra + seesaw, if the focus point region is the correct explanation of the observed 
DM. Note, however, that in the future the top mass will be measured more precisely. At 
the LHC one expects an uncertainty of 1-2 GeV [167]. At a linear collider rrit could be 
determined down to an uncertainty of 100 MeV [168]. 
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Figure 6.6: Allowed region for dark matter density in the {ttiq, M1/2) plane for Aq = 0, ^ > 
and tan/3 = 45, for (from top to bottom) Mt ~ 5 x 10^^ GeV (red), Mt = lO^^(green) and 
Mr = 10^5 GeV (blue). 

We now turn to a discussion of large tan /3. At large values of tan /3 the width of the 
CP-odd Higgs boson A becomes large, Ta ~ Ma tsm^ j3 [ml + 'rn^), and a wide s-channel 
resonance occurs in the region m^o ~ Ma/2. The enhanced annihilation cross section 
reduces ^x°^^ ^° acceptable levels, the resulting region is known as the "Higgs funnel" 
region. In fig. (6.6) we show the allowed range of parameters in the (mo, Afi/2) plane for 
one specific value of tan/? = 45 and three different values of Mq-. As demonstrated, the 
Higgs funnel region is very sensitive to the choice of Mt . It is fairly obvious that varying 
Mt one can cover nearly all of the plane, even for fixed values of all other parameters. 
We have calculated the DM allowed region for various values of tan /3 and found that the 
funnel appears for all tan /3 > 40, approximately. 

The strong dependence of the Higgs funnel region on Mt unfortunately does not 
imply automatically that if large tan j3 is realized in nature one could get a very sensitive 
indirect "measurement" of the seesaw scale by determining (mo, M1/2). The reason is that 
the Higgs funnel is also very sensitive to the exact value of tan/3 and to the values (and 
errors) of the top and bottom quark mass. The latter is demonstrated in fig. (6.7), where 
we show the DM allowed range of parameters for a fixed choice of tan /3 and Mt varying 
to the left (to the right) rrit [nify) within their current 1 a c.l. error band. The position of 
the funnel is especially sensitive to the exact value of m{,. Comparing fig. (6.7) with fig. 
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Figure 6.7: Allowed region for the dark matter density in the (mo, plane for Aq = 0, /.t > 

and tan /3 = 45, for Mr = 5 x 10^^ GeV and (to the left) for three values of rutop = 169.1GeV 
(blue), nitop = 171.2 GeV (red) and rritop ~ 173.3 GeV (green). To the right: The same, but 
varying mi,, nihot = 4.13 GeV (blue), rribot ~ 4.2 GeV (red) and rrihot ~ 4.37 GeV (green). 

(6.6) one can see that the uncertainty in and rrit currently severely limit any sensitivity 
one could get on Mt- However, future determinations of mi, and mt could improve the 
situation considerably. For future uncertainties in rrit see the discussion above for the 
focus point region. For mi, reference [169] estimates that mi, could be fixed to 4.17 ± 0.05 
GeV, which might even be improved to an accuracy of Ami, ~ 16 MeV according to [170]. 
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Figure 6.8: Allowed region for dark matter density in the (rriQ, M1/2) plane for the "standard 
choice" of mSugra parameters for Mt = 10^^ GeV. To the left: For one fixed value of A2 = 0.5 
the allowed range for negligibly small neutrino Yukawa couplings (red) and Yt fitted to correctly 
explain solar and atmospheric neutrino data (blue lines). To the right: the DM allowed range 
of parameters for 3 different values of A2, A2 = 0.5 (red), A2 = 0.75 (green) and A2 = 1 (blue). 
Note the logarithmic scale. 

All of the above figures have been calculated using fixed values for Ai and A2 and neg- 
ligibly small Yukawa couplings Yt- This choice in general does not affect the calculation 
of the DM allowed regions much. However, a fully consistent calculation can not vary 
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Mt, Yt and A2 independently, since this will lead to neutrino masses and angles outside 
the experimentally allowed ranges. Since Yt is diagonalized by the same matrix as the 
effective neutrino mass matrix, (see the previous section) the measured neutrino an- 
gles provide constraints on the relative size of the entries in Yt- The absolute size of Yt 
is then fixed for any fixed choice of A2 and M^, once the neutrino spectrum is chosen to 
be hierarchical or quasi-degenerate. In the numerical calculation shown in fig. (6.8) we 
have chosen neutrino masses to be of the normal hierarchical type and fitted the neutrino 
angles to exact tri-bimaximal (TBM) values [29], i.e. tan^ 6'Atm = 1, tan^ ^0 = 1/2 and 
sin^ 6'r = 0. This has to be done in a simple iterative procedure, since the triplet param- 
eters are defined at the high scale, whereas neutrino masses and angles are measured at 
low scale. For more details on the fit procedure see [113]. 

In fig. (6.8) to the left we show two calculations of the DM allowed regions. The 
allowed range for negligibly small neutrino Yukawa couplings is shown by the filled (red) 
region, while the calculation with Yt fitted to correctly explain solar and atmospheric neu- 
trino data is the one inside the (blue) lines. Note the logarithmic scale. As demonstrated, 
the exact values of Yt are of minor importance for the determination of the parameter 
region allowed by the DM constraint. Slightly larger differences between the fitted and 
unfitted calculations are found pushing Mt to larger values (see, however, below). For 
smaller values of Mt, the entries in Yt needed to correctly explain neutrino data are 
smaller and, thus, Yt affects the DM allowed region even less for Mt < 10^^ GeV. 

In fig. (6.8) to the right we compare three different calculations for A2, A2 = 0.5 
(red), A2 = 0.75 (green) and A2 = 1 (blue), for fixed choice of other parameters. This 
plot serves to show that also the exact choice of A2 is of rather minor importance for the 
determination of the DM allowed region. Very similar results have been found for Ai, we 
therefore do not repeat plots varying Ai here. 

In fig. (6.9) we show the DM allowed parameter regions for tan (3 = 10 and two values 
of Mt, Mt = GeV (to the left) and Mt = 10" GeV (to the right), for a fixed choice 

of all other parameters. Superimposed on this plot are lines of constant branching ratio for 
Br{ii eq). The latter have been calculated requiring neutrino masses being hierarchical 
and fitted to solar and atmospheric neutrino mass squared differences and neutrino angles 
fitted to TBM values. Within the (mo, M1/2) region shown, i?r(/i — > 67) can vary by two 
orders of magnitude, depending on the exact combination of (mo, Afi/2), even for all other 
parameters fixed. The most important parameter determining i?r(/i — > 67), once neutrino 
data is fixed, however, is My, as can be seen comparing the figure to the left with the plot 
on the right. While for Mt = GeV about "half" of the plane is ruled out by the non- 
observation of /i — >■ 67, for Mt = 5- 10^"^ GeV with the current upper limit nearly all of the 
plane becomes allowed. The strong dependence of — > 67 on Mt can be understood from 
the analytical formulas presented in [113]. In this paper it was shown that i?r(/i 67) 
scales very roughly as Br{fi — > 67) oc M|.log(Mr), if neutrino masses are to be explained 
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correctly. For tan /3 = 10 one thus concludes that with present data values of Mt larger 
than (few) 10^^ GeV - (few) 10^^ GeV are excluded by Br{fi — > 67), to be compared 
with Mt/X2 ^ 10^^ GeV from the measured neutrino masses. Note, however, that (i) the 
constraint from neutrino masses is relatively independent of tan/3, rriQ and M1/2, while 
/i — )> e7 shows strong dependence on these parameters; and (ii) allowing the value of the 
reactor angle sin^ to vary up to its experimental upper limit, sin^ 6^ = 0.056 [28], leads 
to larger values of Br{fi — >■ 67) and thus to a tighter upper limit on Mt- 



Mj= 5 X 1 0^^ (GeV) tanp=1 0, A(,=0 (GeV) Mj= 1 0^'' (GeV) tan(3=1 0, A(,=0 (GeV) 




100 200 300 400 500 600 700 800 900 1000 100 200 300 400 500 600 700 800 900 1000 

M,/2 (GeV) Mi/2 (GeV) 

Figure 6.9: Allowed region for dark matter density in the {mo, M1/2) plane for our "standard 

choice" of mSugra parameters and for two values of Mt: Mt = 5 x 10^"^ (left panel) and for 

Mt = 10"'^'* (right panel). Superimposed are the contour lines for the Br{iJ. — > 67). 

6.3 Conclusions 

In conclusion, we have calculated the neutralino relic density in a supersymmetric model 
with mSugra boundary conditions including a type-II seesaw mechanism to explain current 
neutrino data. We have discussed how the allowed ranges in mSugra parameter space 
change as a function of the seesaw scale. The stau co-annihilation region is shifted towards 
smaller mo for smaller values of the triplet mass Mt , while the bulk region and the focus 
point line are shifted towards larger values of M1/2 for Mt sufficiently below the GUT scale. 
The Higgs funnel, which appears at large values of tan j3 has turned out to be especially 
sensitive to the value of Mt- Determining M1/2 from the mass of any gaugino and rriQ 
from a sparticle which is not important for the DM calculation, one could, therefore, get 
a constraint on Mt from the requirement that the observed flDAih^ is correctly explained 
by the calculated Q^o/i^. 

On the positive side, we can remark that current data on neutrino masses put an upper 
bound on Mt of the order of 0(10"^^) GeV. Since this is at least one order of magnitude 
smaller than the GUT scale, the characteristic shifts in the DM regions are necessarily 
non-zero if our setup is the correct explanation of the observed neutrino oscillation data. 
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Even more stringent upper limits on Mt follow, in principle, from the non-observation of 
LFV decays. A smaller Mt implies larger shifts of the DM region. However, the "exact" 
upper limit on Mt from LFV decays depends strongly on tan/3, mo and M1/2, and thus 
can be quantified only once at least some information on these parameters is available. 

On the down side, we need to add a word of caution. We have found that the DM 
calculation suffers from a number of uncertainties, even if we assume the soft masses to be 
perfectly known. The most important SM parameters turn out to be the bottom and the 
top quark mass. The focus point line depends extremely sensitively on the exact value of 
the top mass, the Higgs funnel shows a strong sensitivity on both, mt, and mt. 

Finally, it is clear that quite accurate sparticle mass measurements will be necessary, 
before any quantitative conclusions can be taken from the effects we have discussed. 
Unfortunately, such accurate mass measurements might be very difficult to come by for 
different reasons. In the focus point region all scalars will be heavy, leading to small 
production cross section at the LHC. In the co- annihilation line with a nearly degenerate 
stau and a neutralino, the stau decays produce very soft taus, which are hard for the 
LHC to measure. And the Higgs funnel extends, depending on tan/3 and My, to very 
large values of {mQ, M1/2), at least partially outside the LHC reach. Nevertheless, DM 
provides in principle an interesting constraint on the (supersymmetric) seesaw explanation 
of neutrino masses, if seesaw type-II is realized in nature, a fact which to our knowledge 
has not been discussed before in the literature. 



Chapter 7 



Lepton Flavour Violation and Dark Matter in 
Seesaw Type III Model 

7.1 Introduction 

This chapter is based on tlie work [171]. We study a supersymmetric version of the 
seesaw mechanism type-Ill. The model consists of the MSSM particle content plus three 
copies of 24 superfields. The fermionic part of the SU{2) triplet contained in the 24 
is responsible for the type-Ill seesaw, which is used to explain the observed neutrino 
masses and mixings. Complete copies of 24 are introduced to maintain gauge coupling 
unification. These additional states change the beta functions of the gauge couplings 
above the seesaw scale. Using mSUGRA boundary conditions we calculate the resulting 
supersymmetric mass spectra at the electro-weak scale using full 2-loop renormalization 
group equations. We show that the resulting spectrum can be quite different compared to 
the usual mSUGRA spectrum. We discuss how this might be used to obtain information 
on the seesaw scale from mass measurements. Constraints on the model space due to 
limits on lepton flavour violating decays are discussed. The main constraints come from 
the bounds on fi ^ but there are also regions where the decay r — > /i7 gives stronger 
constraints. We also calculate the regions allowed by the dark matter constraint. For the 
sake of completeness, we compare our results with those for the supersymmetric seesaw 
type-II and, to some extent, with type-I. 

7.2 Supersymmetric seesaw type-Ill 

In the case of a seesaw model type-Ill one needs new fermions S at the high scale belonging 
to the adjoint representation of SU{2). This has to be embedded in a 24-plet to obtain 
a complete SU{5) representation. The superpotential of the unbroken SU{5) relevant for 
our discussion is 

W = V2 ^mY'Wm^h - jWmY'^IOm^h + 5h2AmY^"5m + ]^2AmM2a2Am ■ (7.1) 
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We have not specified the Higgs sector responsible for the SU{5) breaking. The new 
parts, which will give the seesaw mechanism, comes from the 24jv/. It decomposes under 
SU{3) X SU{2) X U{1) as 

24m = (1, 1,0) + (8, 1,0) + (1,3,0) + (3, 2, -5/6) + (3*, 2, 5/6) , (7.2) 
= Bm + Gm + Wm + Xm + Xm ■ 

The fermionic components of (1, 1, 0) and (1, 3, 0) have exactly the same quantum numbers 
as A^'^ and S. Thus, the 24jv/ always produces a combination of the type-I and type-Ill 
seesaw. 

In the SU{5) broken phase the superpotential becomes 

Win = Wmssm + Hu{WmYn- '^BmYb)L + hSmYxD'' 

+^BmMbB,j + ^GmMgGm + ^WmMwWm + XmMxXm (7.3) 

As before we use at the GUT scale the boundary condition Yat = Yb = Yx and Mb = 
Mg = Mw = Mx- Integrating out the heavy fields yields the following formula for the 
neutrino masses at the low scale: 

= "f {^^bMb^Yb + \y^M^'Yw^ . (7.4) 

As mentioned above there are two contributions stemming from the gauge singlet as well 
as from the SU{2) triplet. In this case the calculation of the Yukawa couplings in terms 
of a given high scale spectrum is more complicated than in the other two types of seesaw 
models. However, as we start from universal couplings and masses at Mqut we find that 
at the seesaw scale one still has Mb — My/ and Yb — Yw so that one can write in a good 
approximation 

= -vl^Y^M^'Yw (7.5) 

and one can use the corresponding decomposition for Yy/ as discussed in section 5.1 up 
to the overall factor 4/5. 

7.3 Effects of the heavy particles on the MSSM spectrum 

The appearance of charged particles at scales between the electro-weak scale and the 
GUT scale leads to changes in the beta functions of the gauge couplings [86, 109]. In 
the MSSM the corresponding values at 1-loop level are (61,^2,^3) = (33/5,1,-3). In 
case of one 15-plet the additional contribution is /Shi = 7/2 whereas in case of 24- 
plet it is Afo,; = 5. This results in case of type-II in a total shift of Aft^ = 7 for the 
minimal model and in case of type-Ill in A6j = 15 assuming 3 generations of 24-plets. 
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Figure 7.1: Mass parameters at Q = 1 TeV versus the seesaw scale for fixed high scale 
parameters mo = A/1/2 = 1 TeV, Aq — 0, tan/3 — 10 and /i > 0. The full lines correspond to 
seesaw type-I, the dashed ones to type-II and the dash-dotted ones to type-Ill. In all cases a 
degenerate spectrum of the seesaw particles has been assumed. 



This does not only change the evolution of the gauge couplings but also the evolution 
of the gaugino and scalar mass parameters with profound implications on the spectrum 
[109, 113]. Additional effects on the spectrum of the scalars can be present if some of the 
Yuka-wa couplings get large [113, 172, 173]. In Fig. 7.1 we exemplify this by sho-wing the 
values of selected mass parameters at Q = 1 TeV versus the seesaw scale for fixed high 
scale parameters rriQ = M1/2 = 1 TeV and we have set the additional Yukawa couplings 
to zero. As expected, the effects in case of models of type-II and III are larger the smaller 
the corresponding seesaw-scale is. The scalar mass parameters shown are of the first 
generation and, thus, the results are nearly independent of tan/3 and Aq. For illustration 
we show in Fig. 7.2 the corresponding spectrum where we have fixed tan/3 = 10 and 
Ao = 0. 

We note that in all three model types the ratio of the gaugino mass parameters is 
nearly the same as in the usual mSUGRA scenarios but the ratios of the sfermion mass 
parameters change [109, 113]. One can form four 'invariants' for which at least at the 
1-loop level the dependence on M1/2 and mo is rather weak, e.g. (m| — m\)/Ml, (rriQ — 
m%)/Mf, {m?£, — m\)/Ml and (mg — m\j)/Ml. Here one could replace Mi by any of 
the other two gaugino masses which simply would amount in an overall rescaling. In 
Fig. 7.3 we show these 'invariants' in the leading-log approximation at 1-loop order to 
demonstrate the principal behaviour for seesaw type-II with a pair of 15-plets and seesaw 
type-Ill with three 24-plets. From this one concludes that in principle one has a handle to 
obtain information on the seesaw scale for given assumptions on the underlying neutrino 
mass model, if universal boundary conditions are assumed. For the type-I, i.e. singlets 
only, of course A&j = and no change with respect to mSUGRA are expected. If, for 
example, the seesaw III model would be realized in nature with three 24-plets having 
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Figure 7.2: Example of spectra at Q = 1 TeV versus the seesaw scale for fixed high scale 
parameters mo — M1/2 ~ 1 TeV, tan/3 = 10 and u > 0. On left panel Mh,m^o ,m-+ while on 
the right panel we have M^,m^o,m^+. The line codes are as in Fig. 7.1. 



similar masses around 10^^ GeV one could e.g. sho-w that the corresponding ratios cannot 
be obtained with one pair of 15-plets in the seesaw II model, thus excluding this possibility. 
However, taking the seesaw II with two pairs of 15-plets one would obtain similar ratios 
as in this case the corresponding additional beta-functions at 1-loop would be Abi = 14, 
e.g. nearly equal to our seesaw III model. 

The leading-log approximation gives only the general trend, but there is an important 
dependence on the SUSY point chosen. In Fig. 7.4 we show as illustration {m\ — rn\)/M'l 
and {m?Q — m\)/Ml for different mSUGRA points and at different loop orders: the dashed 
lines are at 1-loop level whereas the solid ones are at 2-loop level. The points considered 
are SPSS [174] with mo = 90 GeV, M1/2 = 400 GeV, Aq = 0, tan/3 = 10, /x > and 
for the same values of Aq and tan/3 two points with M1/2 = 1 TeV: mo = 500 GeV and 
mo = 1 TeV. The black line shows for comparison the leading-log approximation. We 
observe that usually the approximation gets worse for lower values of M24 and this is even 
stronger at the 2-loop level which is a consequence of the large coefficient in the beta 
functions at the 2-loop level. Nevertheless, one sees that in general it gives the correct 
trend, but it might even fail completely, e.g. in the case of Afi/2 = mo = 1 TeV. The 
reason for the drop around — 3.5 x 10^"^ is that the difference between the parameters 
goes to zero as can also be seen from the right of Fig. 7.1, see also discussion below. 

Last but not least we note that the use of the 2-loop RGBs leads to a shift of Mqut 
from about 2x10^^ GeV for 24-plet mass of 10^^ GeV to about 4xW^ GeV for 24-plet 
mass of 10^"^ GeV, which is part of the differences between the 1-loop and 2-loop results 
in Fig. 7.4. Here Mqut is defined as the scale where the electro-weak couplings meet, 
e.g. 5(7(1) = 55(7(2) • This implies also that there is some difference for the strong coupling 
which is, however, in the order of 5-10% which can easily be accounted for by threshold 
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JQl2 JQl3 JQl4 JQl5 JQl6 

Mi5 = M24 [GeV] 

Figure 7.3: Four different "invariant" combinations of soft masses versus tiie mass of the 
15-plet or 24-plet, Afis — Af24- The plot assumes that the Yukawa couplings are negligibly 
small. The calculation is at 1-loop order in the leading-log approximation. The lines running 
faster up towards smaller M are for type-Ill seesaw, the values for type-II seesaw are shown for 
comparison. 



effects of the new GUT particles, e.g. the missing members of the gauge fields and the 
Higgs fields responsible for the breaking of the GUT group [175]. A second reason why the 
deviations between the leading log calculation, the case of 1-loop and 2-loop RGEs gets 
larger for smaller seesaw scale is that the increase of the beta coefficients implies larger 
values of the gauge couplings at the GUT scale. This implies that one reaches a Landau 
pole for sufficiently low values of the seesaw scale. As an example we show in Fig. 7.5 
the value of the gauge coupling at Mqut = 2 x 10"^^ GeV as a function of the seesaw 
scale for type-II with a pair of 15-plets (black lines) and type-Ill with three degenerate 
24-plets (green lines). In both cases the 2-loop RGEs imply a larger gauge coupling for 
a fixed seesaw scale. One sees that in case of type-II (type-Ill) in principle one could 
reach a seesaw scale of about 10^ GeV (10^"^ GeV). However, we believe that we can no 
longer trust even the 2-loop calculation for such large values of the gi, as the neglected 
higher order terms become more and more important. Especially, we should not trust the 
"turn-over" of the invariants in Fig. 7.4 for very low values of the seesaw scale, since the 
numerical calculation at these points is already very close to breaking down. 

We would also like to mention that, in the numerical calculation we find very often 
that one of the scalar masses squared, in particular staus and/or sbottoms, gets large 
negative values already for values of the seesaw scale larger than the Landau pole and 
thus we can not go to values of the seesaw scale as low as the examples shown in Fig. 7.4 
in many SUSY points. 
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A/24 [GeV] A/24 [GeV] 

Figure 7.4: The limits of tlie invariants in seesaw type-Ill models. Left: (m| — m^)/Mf, 
right {rriQ — m|.)/Mf . The blue lines are for SPSS, the light blue one for mo = 500 GeV and 
A/1/2 — 1 TeV, and the red one for mg = = 1 TeV; full (dashed) lines are 2-loop (1-loop) 

results. The black line is the analytical approximation, for comparison. 



7.4 Lepton flavour violation in the slepton sector 



From a one-step integration of the RGEs one gets assuming mSUGRA boundary condi- 
tions a first rough estimate for the lepton flavour violating entries in the slepton mass 
parameters: 

- -^^i^^ + ^Di^N^LYj^) , (7.6) 

for i 7^ j in the basis where is diagonal, = Iu^Mgut / Mi)Sij and is the additional 
Yukawa coupling of the type-A; seesaw at Mqut {k = /, //, ///). We obtain 

9 

a/ = 1 , a// = 6 and aju = - . (7.8) 



Note, that in case of the type-11 the matrix L is degenerate and thus can be factored 

out. All models have in common that they predict negligible flavour violation for the 
right-sleptons 

m%,^ - 0. (7.9) 

We know that these approximations work well only in case of the type-I models. Never- 
theless they give a rough idea on the relative size one has to expect for the rare lepton 
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Figure 7.5: Values of the gauge coupling at Mqut = 2 x 10^^ GeV as a function of the 
seesaw scale, black lines seesaw type-II and green lines seesaw type-Ill with three 24-plets with 
degenerate mass spectrum; full (dashed) lines are 2-loop (1-loop) results. For the calculation of 
the electroweak threshold the spectrum corresponds to rriQ = A/1/2 = 1 TeV, Aq = 0, tan/3 — 10 
and /i > 0. 

decays — >■ Ij-^ which very roughly scale like 

I fii^ j ^ 

Br{li Ij-i) oc a^ml^-^-tan^ p. (7.10) 

where fh is the average of the SUSY masses involved in the loops. Note, that for a given 
set of high scale parameters both, the different size of the flavour mixing entries and the 
changed mass spectrum, play a role. 

7.5 Numerical results 

In this section we present our numerical calculations. All results presented below have 
been obtained with the lepton flavour violating version of the program package SPheno 
[7, 119]. The RGEs of the seesaw II and seesaw III models have been calculated with 
SARAH [176-178]. All seesaw parameters are defined at Mqut and as mentioned in the 
previous section we require for models of type-II the boundary condition Yz = Ys = Yt 
and Mz = Ms = Mt and in case of type-Ill models Y^ = Yb = Y\y and Mb = Mq = 
My/ = Mx- We evolve the RGEs to the scale(s) corresponding to the GUT scale values 
of the masses of the heavy particles. The RGE evolution implies also a splitting of the 
heavy masses. We therefore add at the corresponding scale the threshold effects due 
to the heavy particles to account for the different masses. In case of type-Ill models 
off-diagonal elements are induced in the mass matrices. This implies that one has to 
go the corresponding mass eigenbasis before calculating the threshold effects. We use 
2-loop RGEs everywhere except stated otherwise. In the appendix we give the necessary 
ingredients on how to obtain them in the seesaw type-II and III models. The analogous 
anomalous dimensions for the type-I model can be found in [115]. 

Unless mentioned otherwise, we fit neutrino mass squared differences to their best fit 
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values [14] and the angles to tri-bi-maximal (TBM) values [29] . Our numerical procedure 
is as follows. Inverting the seesaw equation, see eqs. (6.11) and (7.4), one can get a 
first guess of the Yukawa couplings for any fixed values of the light neutrino masses (and 
angles) as a function of the corresponding triplet mass for any fixed value of the couplings. 
This first guess will not give the correct Yukawa couplings, since the neutrino masses and 
mixing angles are measured at low energy, whereas for the calculation of we need to 
insert the parameters at the high energy scale. However, we can use this first guess to run 
numerically the RGEs to obtain the exact neutrino masses and angles (at low energies) for 
these input parameters. The difference between the results obtained numerically and the 
input numbers can then be minimized in a simple iterative procedure until convergence is 
achieved. As long as neutrino Yukawas are VF^j < 1 we reach convergence in a few steps. 
However, in seesaw type-II and type-Ill the Yukawas run stronger than in seesaw type-I, 
so our initial guess can deviate sizable from the correct Yukawas, implying in general 
also more iterations until full convergence is reached. Since neutrino data requires at 
least one neutrino mass to be larger than about 0.05 eV, we do not find any solutions for 
^ A2 X 10^^ GeV and M24 > 8 x lO"^"* GeV, respectively. In the latter case we have 
assumed that all 24-plets have similar masses. For sake of completeness we note that one 
can also satisfy all neutrino data by giving one of the 24-plets a large mass in the order 
of Mgut or larger having a model with effectively only two 24-plets. 

7.5.1 Lepton flavour violation 

We have seen in eq. (7.10) that rates for the lepton fiavour violating decays of /i and r 
scale like the LFV entries in the slepton mass squared matrix squared and inverse to the 
overall SUSY mass to the power eight. From this one immediately concludes the rates for 
the rare lepton decays are in general larger in seesaw models of type-II and HI than in 
type-I models for fixed SUSY masses and seesaw scales except if one arranges for special 
cancellations. 

Comparing the type-II with the type-Ill model one finds that LFV decays are larger 
for type-Ill, as shown for the case of — )> 67 in Fig. 7.6. From eqs. (7.7) and (7.8), 
however, one would expect that type-II should have larger LFV. Numerically we find 
the opposite for two reasons, (i) Br{li — >■ Ij'y) strongly depends on the SUSY masses, 
see eq. (7.10) and type-Ill has a lighter spectrum than type-II (for the same mSUGRA 
input parameters). And (ii) 2-loop effects are very important in type-Ill, due to the large 
coefficients, in general leading to large fiavor violating soft SUSY breaking parameters. 

In Fig. 7.6 we compare Br{li — > Ijj) for the three seesaw models taking degenerate 
seesaw spectra in case of type-I and type-Ill. Note that in case of seesaw type-Ill we 
can only show a relatively short interval for the seesaw scale which is mainly due to two 
reasons: (i) for scales below approximately 10^'^ GeV the gauge couplings get large at 
Mqut as a consequence of the large beta functions and, thus, perturbation theory breaks 
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mo=M„2=300 (GeV), tanp=10, Ao=0 (GeV) m(,=M„2=1000 (GeV), tanp=10, Ao=0 (GeV) 




Figure 7.6: Br{fi — > 67) as a function of the seesaw scale for seesaw type-I (red line), seesaw 
type-II (blue line) and seesaw type-III (magenta line). In case of type-I and type-Ill a degenerate 
spectrum has been assumed. On the left panel mo — = 300 (GeV), on the right panel 

mo — rni/2 = 1000 (GeV). In both cases we take tan/3 — 10, Ao=0 and /_i > 0. 



mo 




m-+ 




m-g 










500 


178 


333 


617 


1029 


535 


543 


600 


772 


1000 


180 


338 


642 


1057 


1008 


1020 


1043 


925 



Table 7.1: Examples masses in GeV for M1/2 = 1000 GeV, tan fi = IQ, Aq = Q GeV and /i > 0, 
for seesaw type-III for a degenerate seesaw spectrum with M24 = 10^'^ GeV. 

down, (ii) One encounters negative mass squares for the scalars, in particular for the 
lighter stau and/or lighter sbottom. The latter point is also the reason why the possible 
range is larger in case of the larger soft SUSY breaking parameters. 

The values for i?r(/i — > 67) in Fig. 7.6 are larger than the current experimental bound 
[56], so one might worry if in case of type-III models only SUSY spectra beyond the 
reach of the LHC are allowed. (Note, that even for the examples shown the masses of 
the sfermions are already in the range of several hundred GeVs as can be seen from 
table 7.1.) Indeed we find that by putting generic Yukawa couplings which are able to 
explain neutrino data one needs a heavy spectrum to be consistent with bounds on the 
rare lepton decays. However, this is strictly true only for the TBM angles and R =1. 
Accidental cancelations due to different contributions to the flavor violating soft masses 
and thus to the rare lepton decays are possible in type-III (and in type-I). As an example 
we show in Fig. 7.7 Br{^ — >■ 67) as a function of the reactor angle 5^3 for different values 
of the Dirac phase 5. For comparison we also show the calculation for a type-I model. 
For 6 = Ti there is a range of where this branching ratio is below the experimental 
constraint. 

At flrst glance this seem to require some flne-tuning of the underlying parameters. 
However, one can look at this from a different perspective; Assume that the MEG col- 
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m(,=M„2=1000 (GeV), tanp=10, A(,=0 (GeV) 




„2 
= 13 

Figure 7.7: Br{^j, -> ej) versus sfg for mo = M1/2 = 1000 GeV, tan/3 = 10, Aq = GeV 
and /i > 0, for seesaw type-I (solid lines) and seesaw type-Ill (dashed lines), for il/geesaw = 
GeV. The curves shown are for 2 values of the Dirac phase: 6 — (red) and 6 = n (blue), both 
for normal hierarchy. 

laboration has found a non- vanishing value for i?r(/i e-y) and from LHC data one has 
found that the spectrum is consistent with the type-Ill seesaw model. For a fixed R- 
matrix, e.g. R=l one would obtain in this case a relation between s^g and M24. This can 
be exploited to put a bound on M24 or even to determine it depending on the outcome 
of measurements of reactor angle and, thus, the model assumptions can be tested. In 
Fig. 7.8 we show the corresponding rare tau decays. Note that also for r — > 67 such a 
cancelation exists in principle but the corresponding range is excluded hy fi ^ 67. In 
contrast r — >■ //7 is insensitive to the reactor angle and should be measurable in the near 
future. 

Up to now we have assumed that the seesaw spectrum is nearly degenerate which 
is of course a strong assumption. We show in Fig. 7.9 two examples where we keep in 
each case two masses fixed and vary the third one. Note, that in contrast to SUSY 
particles the indices of the heavy particles are generation indices and do not correspond 
to a particular mass ordering, e.g. Mji^ corresponds to the 'solar neutrino scale' and Mr^ 
to the 'atmospheric neutrino scale'. In case that the mass of the first generation state is 
varied, e.g. the left plot of this figure, one finds a decrease of the branching ratios with 
increasing seesaw mass M/j^ . This is mainly caused by an increase of the SUSY spectrum 
while at the same time neutrino physics is only affected mildly requiring only a light 
increase of the corresponding Yukawa couplings to obtain the correct neutrino masses. If, 
on the other hand, the mass M^j^ of the third generation seesaw particles is increased on 
needs also a sizable increase of the Yukawa couplings to obtain the correct neutrino mass 
difference squared for the atmospheric sector. This leads to the observed behaviour that 
the branching ratios for r — )> /i7 and r — )> 3yU increases while the other ones decrease. 
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mo=Mi,2=1 000 (GeV), tan|3=1 0, Ao=0 (GeV) m(,=M„2=1 000 (GeV), tanp=1 0, Ao=0 (GeV) 




Figure 7.8: i3r(r — > ej) versus 5^3 (left) and i?r(r — > /17) versus 5^3 (right) for mo = M1/2 ~ 
1000 GeV, tan (3 = 10, Aq — GeV and > 0, for seesaw type-I (solid lines) and seesaw type-III 
(dashed lines), for Mseesaw = 10^"* GeV. The curves shown are for 6 — (red) and 5 = tt (blue) 
for normal hierarchy. 

7.5.2 Dark Matter 

The changes in the spectrum induced by the new heavy states also impact on the pre- 
dictions with respect to the relic density which we have calculated using the program 
micrOMEGAs [161]. As is well-known, within mSUGRA there are 4 regions in parameter 
space, in which the constraint from dark matter can be satisfied. These are (i) the bulk 
region; (ii) the stau co-annihilation region; (iii) the focus point line and (iv) the Higgs 
funnel. Below we will show usually the range of Qh^ allowed at 3 a according to [56] 

0.081 < fl/i^ < 0.129. (7.11) 

In particular, the co-annihilation region is very sensitive to the difference between 
the masses of the lightest stau and the lightest neutralino. In Fig. 7.10 we observe that 
this difference depends strongly on the seesaw scale in both models. For a fixed M1/2 
and mo lowering the seesaw scale increases this mass difference, which then leads to a 
larger calculated flh"^. To compensate for this effect one needs to lower mo, with the 
value depending on the seesaw scale chosen. For certain seesaw scales then mo needs to 
be lowered below mo = and the co- annihilation region disappears. In this region of 
parameter space both models behave in a qualitatively similar way. However, recall that 
spectra run faster towards smaller masses in seesaw type-III. 

Also the focus point region is very sensitive to the precise values of the input pa- 
rameters. The focus point region appears in mSUGRA for large values of mo and 
small/moderate values of M1/2 of the order of O(IOO) GeV, the exact value depending 
on mo. This can be seen in figs. 7.11 and 7.12 where we show m^o, the higgsino content 
A'^isP + |A^i4p and the corresponding Qh^ as a function of mo for a fixed seesaw scale 
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mo=M„2=1000 (GeV), tanp=10, Ao=0 (GeV) mo=Mi,2=1000 (GeV), tan(i=10, Ao=0 (GeV) 




Mr, (GeV) Mr^ (GeV) 



Figure 7.9: Branching ratios for U — )• Ijj (solid lines) and U, — )• 3lj (dashed lines) versus the 
seesaw scale for tan jS = 10, fi > 0, Aq = GeV, M^j^ = 'tiq = 1000 GeV. On the left panel we 
scan on with M^^ = M^.^ = 2 x 10^^ GeV while on the right panel we scan on M^^ with 
Mfi^ — = 2 y. 10^"^ GeV. The color code is red for ^ — > 67 or /i — > 3e, blue for r — > /i7 or 
T — 7- 3/i and green for r — )■ 67 or r — )■ 3e. 



Mt,w = 10^^ GeV, Ao = 0, tan/3 = 10, /i > and various values of M1/2. Note, that 
we take different values of M1/2 for the two models in such a way that we obtain similar 
values for m^o . We find that both models behave differently in this region of parameter 
space, e.g. the higgsino content |A'^i3p + \ Nii\^ decreases (increases) with increasing values 
mo for seesaw type-II (type-Ill). However, also for type-II the higgsino content increases 
for increasing niQ once we reach the multi-TeV range but we did not get correct elec- 
troweak symmetry breaking in case of multi-TeV values for mo in case of type-Ill models. 
The increased higgsino content of the lightest neutralino leads to on increase (decrease) 
of its couplings to the Z-boson and the light Higgs boson (to sfermions) resulting in the 
observed dependence of O/i^ for mo close to the 1-TeV region. 

With these observations it is clear that the DM allowed regions will be shifted in the 
mo-Mi/2 plane compared to the usual mSUGRA expectations. We fix in the following 
fTitop = 171.2 GeV, tan/3 = 10, Aq = {) and /i > as well as the seesaw scale to 10^^ 
GeV. For comparison we show in Fig. 7.13 the usual mSUGRA case without any heavy 
intermediate particles (left plot) as well as the case of a seesaw type-I scenario (right plot). 
The blue bands show the 3a range according to [56] and we see the three usual regions: 
the stau co-annihilation with a lighter stau mass close to the LSP mass for Afi/2 < 300 
GeV, the bulk region for moderate values of Afi/2 and mo resulting in small sfermion 
masses as well as the focus point region for M1/2 — 170 GeV and large values of mQ. In 
addition, we show the lines corresponding to Mh = 110 GeV and 114 GeV. Note, that the 
theoretical uncertainty on Mh is still of the order of 3-5 GeV [164, 165]. Moreover, the 
value of the Higgs boson mass also depends strongly on Aq and in particular for negative 
values of Aq one can easily increase the value of Mh while the DM allowed regions hardly 
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Figure 7.10: Difference between tiie masses and ttie liglitest stau and tlie liglitest neutralino 
(upper row) as well as the corresponding flh^ (lower row) as a function of the seesaw scale. The 
left (right) plots are for seesaw type-II (III). A degenerate seesaw spectrum has been assumed in 
case of seesaw type-Ill. = 800 GeV, Aq = 0, tan/3 = 10 and n > 0. The lines correspond 

to full blue line mo = 0, red dashed line toq = 50 GeV, green dashed dotted line mo = 100 GeV, 
black dashed line mo = 150 GeV and orange full line mo — 200 GeV. The gray band shows the 
preferred range according to eq. (7.11). 



change. 

The part of parameter space most affected is the one at large niQ. Since in mSUGRA fi 
is calculated from the requirement of correct electroweak symmetry breaking, /x changes 
rapidly in this region. With the Higgsino content in the lightest neutralino changing 
rapidly as a function of /i, this region is then very sensitive to any changes of parameters. 
Since the Yi, also impacts on the running of the Higgs mass parameters and thus slightly 
affects the value predicted for ^, some small changes are found relative to mSUGRA here. 
Note, however, that this region is highly constrained by the lower bound on the lightest 
chargino mass of the order of 103 GeV [179]. 

In case of the other two seesaw models the shift of the allowed regions is much more 
pronounced, as discussed above. In Figs. 7.14 and 7.15 we show to regions for type-II 
(left plot) and type-Ill (right plot) and two different values for Aq. As claimed above, 
tlie Higgs mass bounds gets shifted significantly while the DM allowed regions are hardly 
affected. As expected the effects are much more pronounced in case of type-Ill as the 
effects of the heavy particles on the spectrum is much stronger. Note, that in particular 
the bending of the allowed region for large rriQ is due to the changed higgsino content as 
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Figure 7.11: Mass of the lightest neutralino (left plot), its higgsino content (middle plot) 
and the corresponding Qh"^ (right plot) as a function of mo for a seesaw type-II model with 
Mt = 10^^ GeV, mtop = 171.2 GeV, Aq = 0, tan/3 = 10 and fi > 0. The lines correspond 
to full blue line A/1/2 = 195 GeV, red dashed line M1/2 ~ 200 GeV, green dashed dotted line 
My2 = 205 GeV, black dashed line M1/2 = 210 GeV and orange full hue M1/2 = 215 GeV. The 
gray band shows the range eq. (7.11). 




Figure 7.12: Mass of the lightest neutralino (left plot), its higgsino content (middle plot) 
and the corresponding O/i^ (right plot) as a function of mo for a seesaw type-Ill model with a 
degenerate seesaw scale Mw = 10^* GeV, mtop = 171.2 GeV, Aq — 0, tan/3 = 10 and 1^1 > 0. 
The lines correspond to full blue line M1/2 ~ 400 GeV, red dashed line M1/2 ~ 405 GeV, green 
dashed dotted line M1/2 = 410 GeV, black dashed line M1/2 ~ 415 GeV and orange full line 
Ml/2 = 420 GeV. The gray band shows the range eq. (7.11). 



discussed in case of figs. 7.11 and 7.12. Moreover, the case of stau co-annihilation is not 
viable anymore in case of the type-Ill model already for this value of the seesaw scale. For 
completeness we mention that for the type-II the stau co-annihilation region disappears 
(below Ml/2 = 1500 GeV) for Mt < 10^"^ GeV. For completeness we note that the results 
here differ slightly from the ones of our previous work [6] because (i) of the corrections 
of the 1-loop RGBs of ref. [86] by [180] and (ii) the complete set of 2-loop RGEs are now 
used. 

In the case of large tan j3 an additional region, usually called the Higgs funnel, opens 
up. This region is characterized by Ma — 2m^o. Also here the regions gets shifted 
compared to usual mSUGRA scenario. However, this region is very sensitive to higher 
order corrections and therefore it is quite important to use full 2-loop RGEs as can be 
seen in Fig. 7.16. We have again fixed Aq = 0, fi > 0, rritop = 171.2 GeV and the seesaw 
scale to lO^'' GeV, with a degenerate spectrum in case of the type-Ill model. The main 
reason for the observed and rather surprisingly large differences between the different 
calculations is that the 2-loop contributions decrease the neutralino mass compared to 



7.6 Conclusions 



155 



tanp=10, Ao=0 (GeV) Mr= 10^" (GeV) tanp=10, Ao=0 (GeV) 




My2 (GeV) My2 (GeV) 

Figure 7.13: Dark matter allowed region (in blue) for mSUGRA (left panel) and for type-I 
seesaw (right panel). The parameters are tan /3 = 10, Aq = 0, fj, > and Mt = 10^"^ GeV for 
iTT-top ~ 171.2 GeV. Also shown (in yellow) are the regions excluded by LEP (small values of 
M1/2), and by LSP constraint (small values of mo). Also shown are the Higgs boson mass curves 
for Mh = 110 GeV (in red) and for M/, = 114.4 GeV (in magenta). 



the 1-loop case while at the same time increasing Ma- For example, in case of seesaw II 
and for fixed values of rriQ = M1/2 = 1500 GeV we get in case of 1-loop RGEs m^o = 560 
GeV, Ma = 1090 GeV and in case of 2-loop RGEs m^o = 498 GeV, Ma = 1100 GeV. For 
completeness we note that this region is also very sensitive to input values for mt and mh 
[6]. 



7.6 Conclusions 

To summarize, we have investigated in detail a supersymmetric version of a seesaw model 
of type-Ill and compared it to seesaw models of type-I and type-II. In case of type-II 
and type-Ill models we have embedded the SU{2) triplets in the corresponding SU{5) 
representations to maintain gauge coupling unification, e.g. 15-plets in case of type-II and 
24-plets in case of type-Ill models. For definiteness we have assumed mSUGRA boundary 
conditions for the soft SUSY breaking parameters. 

The additional heavy charged states lead to changes in the beta-functions and, thus, 
also in the running of the SUSY mass parameters. We have calculated the soft masses 
as a function of the seesaw parameters. As discussed in some detail, there are certain 
combinations of soft masses, which are approximately constants over large regions of 
mSUGRA space. These "invariants" contain indirect information about the seesaw scale 
assuming the type of seesaw model. In certain parts of the parameter space, e.g. for 
low seesaw scales, one might even be able to exclude certain seesaw models by combining 
mass measurements at the LHC with the mSUGRA paradigm. We note, that using 2-loop 
RGEs will be crucial to obtain reliable results. 
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Mj= 10^" (GeV) tanp=10, Ao=0 (GeV) M„^= 10^'* (GeV) tan|3=10, Ao=0 (GeV) 
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M„2 (GeV) My2 (GeV) 

Figure 7.14: Like in Fig. 7.13 but for seesaw type-II (left panel) and type-Ill (right panel). 




200 400 600 800 1000 1200 250 500 750 1000 1250 1500 

My2 (GeV) Mi,2 (GeV) 

Figure 7.15: Like in Fig. 7.14 but for — —300. Seesaw type-11 (left panel) and type-Ill 
(right panel). 

The changes in the spectrum leads obviously to changes in the phenomenology. We 
have calculated lepton flavour violating observables, such as -Br(/j — )> /j +7). We find that 
for fixed (degenerate) seesaw scale these branching ratios are in general largest for type-Ill 
models followed by type-II and type-I. This is a consequence of the fact that for a given 
set of mSUGRA parameters the spectrum in type-Ill is lighter than for type-II models 
which is again lighter than in type-I models. However, the difference in the predictions of 
type-II and type-Ill is somewhat smaller than expected from these considerations because 
in type-II models the flavour violating entries are larger compared to the case of type-Ill 
models. 

We also investigated the predictions for the relic density Qh^ in the type-Ill model and 
compared them with the other models. We flnd the usual four regions in the mSUGRA 
parameter space but of course they are shifted due to the changes in the spectrum. It has 
been found that in particular in case of the Higgs-funnel the use of 2-loop RGEs is crucial 
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Mj= 10^" (GeV) tanp=52, Ao=0 (GeV) Mw^= 10^" (GeV) tanp=49, Ao=0 (GeV) 




500 1000 1500 2000 2500 500 1000 1500 2000 2500 3000 3500 

M„2 (GeV) My2 (GeV) 

Figure 7.16: Comparison between using 1-loop (red) or 2-loop (blue) RGEs on the dark matter 
allowed region for type-II (left panel) and type-III (right panel). The parameters are: Aq — 0, 
/i > and Mseesaw = 10^"* GeV, mtop = 171.2 GeV and tan/3 = 52 for type-II and tan/3 = 49 
for type-III. 

to identify the correct allowed region. Last but not least we note, that for low seesaw 
scales the co-annihilation region vanishes for both, the type-II and the type-III models, 
as the required mass difference between the lightest neutralino and the stau cannot be 
obtained anymore. 



Chapter 8 

Conclusions 



It is now an established fact that neutrinos have mass. This can have important conse- 
quences in the development of Physics in the near future. The quest for an explanation 
why these masses are so tiny puts some stress on the usual Higgs mechanism and asks 
for new Physical frameworks that, in turn, can have important phenomenological conse- 
quences beyond the neutrino sector. We have explored this in the present thesis, con- 
straining the various predictions on Lepton Flavour Violation, Dark Matter and Baryon 
Asymmetry from the neutrino phenomenology. 

In Chapter 3 we started by adopting the possibility of an A4 symmetry in the leptonic 
and Higgs sector, as proposed in [46]. The Higgs sector is extended with an SU{2) 
triplet and therefore a mixed Type-I/Type-II seesaw mechanism was proposed that could 
reproduce a zero texture classified as Bi for the neutrino mass matrix and that is well 
motivated in explaining neutrino data. Then we explored the triplet decay into leptons 
and the consequent production of a lepton asymmetry to obtain the phenomenological 
baryonic asymmetry through the leptogenesis mechanism. To fully address the question 
of the possible washout of the leptonic asymmetry we have solved numerically the full 
set of Boltzmann equations and showed that it is possible to obtain the correct baryonic 
density with the assumptions of this model. Finally we discussed the problem of the 
Higgs potential minimization and substantiated with a full computation demonstrating 
that there is a limitation in obtaining the correct vacuum expectation values for the scalar 
fields within the A4 symmetry and with the representation assignments made therein. We 
ended by proposing two solutions: breaking explicitly the symmetry group and adding a 
second SU{2) triplet. This could find a natural framework in a supersymmetrization of 
the model, although it would make it more elaborate. 

In Chapter 4, based on [22], we discussed an alternative model to the one proposed in 
the previous chapter. We have studied the possibility that the seesaw model with spon- 
taneously broken ungauged lepton number may simultaneously account for the observed 
neutrino masses and mixing as well as the dark matter of the Universe. We have presented 
a two-texture structure for the neutrino mass which arises in a inverse seesaw mechanism 
implementing an A4 flavor symmetry. A predictive pattern of neutrino masses emerges 
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from the interplay of type-I and type-II seesaw contributions, with a lower bound on the 
neutrinoless double beta decay rate, which correlates with the deviation from maximality 
of the atmospheric mixing angle 623, as well as nearly maximal CP violation, correlated 
with the reactor angle Oi^. 

On the other hand, assuming that associated Majoron picks up a mass due to explicit 
lepton number violating effects that may arise, say, from quantum gravity, we showed 
how it can constitute a viable candidate for decaying dark matter, consistent with cosmic 
microwave background lifetime constraints that follow from current WMAP observations. 
We have also shown how the Higgs boson triplet, whose existence is required by the 
consistency of the model, plays a key role in providing a test of the decaying Majoron 
dark matter hypothesis, implying the existence of a mono-energetic emission line which 
arises from the sub-leading one-loop-induced decay of the Majoron into photons. We 
also discussed the possibility of probing its existence in future X-ray observations such as 
expected in NASA's Xenia mission [95]. The presence of the type-II seesaw Higgs triplet 
would also have other particle physics implications, such as lepton flavor violating decay 
rate enhancements due to the strong renormalization effects of the triplet, quite sizeable 
in a supersymmetric model. 

In Chapter 5, based on [5], we started by stressing that low energy neutrino exper- 
iments, including oscillation studies and neutrinoless double-beta decay searches may, 
optimistically, determine at most 9 neutrino parameters: the 3 neutrino masses, the 3 
mixing angles and potentially the 3 CP violating phases. This is insufficient to fully 
reconstruct the underlying mechanism of neutrino mass generation. Then, under the as- 
sumption that neutrino masses arise a la seesaw, we have considered the simplest pure 
type-I or pure type-II seesaw schemes in mSugra. We performed a full scan over the 
mSugra parameter space in order to identify regions where LFV decays of X2 can be max- 
imal, while still respecting low-energy constraints that follow from the upper bounds on 
Br(/x — >• 67). We have also estimated the expected number of events for X2 ~^ Xi + + A*) 
for a sample luminosity of C = 100 fb~^. The expected number of events for the other 
channel X2 ~^ Xi + + ^ is always smaller, as can be seen from the LFV branching ratios 
presented in section 5.2.1. We have found that the pure seesaw-II scheme is substantially 
simpler and comes closer to being fully reconstructible, provided additional LFV decays 
are detected and some supersymmetric particles are discovered at the Large Hadron Col- 
lider, providing the necessary parameter reconstruction information of the supersymmetric 
lagrangian. 

Note that in what concerns the expected number of events both type-I and type- 
II schemes give similar results. However, as we have seen, given their smaller number 
of parameters, type-II seesaw schemes are more likely to be reconstructible through a 
combination of low energy neutrino measurements, with the possible detection of super- 
symmetric states and lepton flavour violation at the LHC. This should encourage one to 
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perform full-fledged dedicated simulations, in order to ascertain their feasibility within 
realistic experimental conditions [96]. 

We note that we have not done a thorough analysis of the fact that the LFV might 
induce new "edge variables", giving additional information [128]. We have focused here on 
LHC, but we should mention that a future ILC would be much more suited for measuring 
LFV SUSY processes [112, 129-134]. 

In Chapter 6, based on [6], we have calculated the neutralino relic density in a super- 
symmetric model with mSugra boundary conditions including a type-II seesaw mechanism 
to explain current neutrino data. We have discussed how the allowed ranges in mSugra 
parameter space change as a function of the seesaw scale. The stau co-annihilation region 
is shifted towards smaller rriQ for smaller values of the triplet mass My, while the bulk 
region and the focus point line are shifted towards larger values of M1/2 for My sufficiently 
below the GUT scale. The Higgs funnel, which appears at large values of tan /? has turned 
out to be especially sensitive to the value of A/y. Determining M1/2 from the mass of 
any gaugino and nio from a sparticle which is not important for the DM calculation, one 
could, therefore, get a constraint on Mt from the requirement that the observed ^l^Aih^ 
is correctly explained by the calculated fl^oh^. 

On the positive side, we can remark that current data on neutrino masses put an upper 
bound on Mt of the order of 0(10"^^) GeV. Since this is at least one order of magnitude 
smaller than the GUT scale, the characteristic shifts in the DM regions are necessarily 
non-zero if our setup is the correct explanation of the observed neutrino oscillation data. 
Even more stringent upper limits on Mt follow, in principle, from the non-observation of 
LFV decays. A smaller Mt implies larger shifts of the DM region. However, the "exact" 
upper limit on Mt from LFV decays depends strongly on tan/3, mo and M1/2, and thus 
can be quantified only when at least some information about these parameters is available. 

On the down side, we need to add a word of caution. We have found that the DM 
calculation suffers from a number of uncertainties, even if we assume the soft masses to be 
perfectly known. The most important SM parameters turn out to be the bottom and the 
top quark mass. The focus point line depends extremely sensitively on the exact value of 
the top mass, the Higgs funnel shows a strong sensitivity on both, nif, and nit. 

It is clear that quite accurate sparticle mass measurements will be necessary, before any 
quantitative conclusions can be taken from the effects we have discussed. Unfortunately, 
such accurate mass measurements might be very difficult to come by for different reasons. 
In the focus point region all scalars will be heavy, leading to small production cross section 
at the LHC. In the co-annihilation line with a nearly degenerate stau and a neutralino, the 
stau decays produce very soft taus, which are hard for the LHC to measure. And the Higgs 
funnel extends, depending on tan/3 and Mt, to very large values of (mo,Mi/2), at least 
partially outside the LHC reach. Nevertheless, DM provides in principle an interesting 
constraint on the (supersymmetric) seesaw explanation of neutrino masses, if seesaw type- 
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II is realized in nature, a fact which to our knowledge has not yet been discussed in the 
literature. 

Finally, in Chapter 7, based on [171], we have studied LFV processes and Dark Matter 
in the framework of a seesaw Type III model and compare it to Type II and in some extent 
to Type I. In case of type-II and type-Ill models we have embedded the SU{2) triplets in 
the corresponding SU{b) representations to maintain gauge coupling unification, e.g. 15- 
plets in case of type-II and 24-plets in case of type-Ill models. For definiteness we have 
assumed mSUGRA boundary conditions for the soft SUSY breaking parameters. The 
additional heavy charged states lead to changes in the beta-functions and, thus, also in the 
running of the SUSY mass parameters. We have calculated the soft masses as a function 
of the seesaw parameters. As discussed in some detail, there are certain combinations 
of soft masses, which are approximately constants over large regions of mSUGRA space. 
These "invariants" contain indirect information about the seesaw scale assuming the type 
of seesaw model. In certain parts of the parameter space, e.g. for low seesaw scales, one 
might even be able to exclude certain seesaw models by combining mass measurements 
at the LHC with the mSUGRA paradigm. We note, that using 2-loop RGBs will be 
crucial to obtain reliable results, he changes in the spectrum leads obviously to changes 
in the phenomenology. We have calculated lepton flavour violating observables, such 
as Br{li — >■ Ij + 7). We find that for fixed (degenerate) seesaw scale these branching 
ratios are in general largest for type-Ill models followed by type-II and type-I. This is 
a consequence of the fact that for a given set of mSUGRA parameters the spectrum in 
type-Ill is lighter than for type-II models which is again lighter than in type-I models. 
However, the difference in the predictions of type-II and type-Ill is somewhat smaller than 
expected from these considerations because in type-II models the flavour violating entries 
are larger compared to the case of type-Ill models. We also investigated the predictions 
for the relic density Q.h? in the type-Ill model and compared them with the other models. 
We find the usual four regions in the mSUGRA parameter space but of course they are 
shifted due to the changes in the spectrum. It has been found that in particular in case of 
the Higgs-funnel the use of 2-loop RGBs is crucial to identify the correct allowed region. 
Last but not least we note, that for low seesaw scales the co-annihilation region vanishes 
for both, the type-II and the type-Ill models, as the required mass difference between the 
lightest neutralino and the stau cannot be obtained anymore. 

In summary we have considered some direct and indirect consequences of models 
leading to massive neutrinos, with a two-way perspective. From current phenomenological 
data we tried to obtain a hint on the actual high scale framework, adopting some previous 
models [46, 75, 86], and have constrained its parameters, using a couple of simplifying 
assumptions, like setting to the identity the Casas-Ibarra i?-matrix [21] - see (1.188) - 
assuming the Tri-Bimaximal mixing ansatz [29] - formula (1.228), the Normal Hierarchy 
for neutrino masses and setting to zero the Dirac and Majorana CP-phases in Chapters 
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5 and 6. Then, with the model defined at GUT scale and within these assumptions, we 
did a study of some of its phenomenological implications at the electroweak scale and 
below for the present baryonic density in the Universe in Chapter 3, for lepton flavour 
violation processes in Chapter 5 and for Dark Matter in Chapters 4 and 6. There are 
still many open questions regarding these matters, like what is the magnitude of CP- 
violation, if any, on the leptonic sector, or if it is seesaw and which type the actual 
mechanism that leads to such tiny masses for neutrinos, or what is the magnitude of 
lepton flavour violation (does it remains within the predictions of the Standard Model 
extended minimally to include neutrino masses or is it there some new Physics), or else 
why there is an asymmetry between matter and anti-matter in the Universe and what 
is the mechanism that generates it... LHC surely will give answers, some partial, some 
deflnitive, as much as there are deflnitive answers in Physics. 
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